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Zusammenfassung 

Das Mittelungsproblem in der Allgemeinen Relativitatstheorie besteht in der Definition 
eines wohldefinierten Mittels iiber TensorgroBen und wir beleuchten dieses Problem von ver- 
schiedenen Seiten. Zunaclist gelien wir auf die kosmologische Riickreaktion ein, die dadurcli 
verursacht wird, dass der gemittelte Einstein- Tensor nicht identisch zu dem Einstein- Tensor 
der gemittelten Metrik ist. Es gibt Vermutungen, nach denen diese Riickreaktion die 
Erklarung fiir Dunkle Energie sein soil. Wir zeigen numerisch, dass im Buchert Formalis- 
mus die Korrekturen von (quasi) linearen Storungen nur von der GroBenordnung 10~^ sind 
und die Eigenschaften von Dunkler Materie aufweisen. AnschlieBend beschaftigcn wir uns 
mit der Formulierung eines allgemein kovarianten Mittelungsprozesses, der die Metrik in 
Vielbeine zerlegt und diese mit Hilfe eines relativistischen Wegner- Wilson Operators an 
einen gemeinsamen Punkt parallelverschiebt, wo sie anschlieBend gemittelt werden. Fiir 
die Festlegung des entsprechenden Vielbeinfeldes wird der Lagrange- Formalismus verwen- 
det. Die Funktionsweise des Mittelungsprozesses wird an speziellen Beispielen in zwei 
und drei Raumdimensionen verdeutlicht. Dazu werden partielle Differentialgleichungen 
numerisch mit dem Simulationspaket Gascoigne gelost. 



Metric Renormalization in General Relativity 

Abstract 

The averaging problem in general relativity concerns the difficulty of defining meaningful 
averages of tensor quantities and we consider various aspects of the problem. We first 
address cosmological backreaction which arises because the averaged Einstein tensor is 
not the same as the Einstein tensor of the averaged metric. It has been suggested that 
backreaction might account for the dark energy. We show numerically in the Buchert 
formalism that the corrections from (quasi) linear perturbations are only of the order of 
10~^ and act as a dark matter. We then focus on constructing averaged metrics and 
present a generally covariant averaging process which decomposes the metric into Vielbeins 
and parallel transports them with a relativistic Wegner-Wilson operator to a single point 
where they can then be averaged. The Vielbeins are chosen in a Lagrangian formalism. 
The functionality of the process is demonstrated in specific examples in two and three 
space dimensions. This involves the numerical solution of partial differential equations by 
the aid of the simulation toolkit Gascoigne. 
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CHAPTER 1 



INTRODUCTION 



The averaging problem is one of the greatest unresolved problems of general relativity. It 
concerns the question of how to define an average over tensor quantities in a physically 
sensible manner when they are defined on separate spacetime points. Such an average can 
only be easily defined for scalar quantities and most of the proposed averaging techniques 
thus far have failed to average higher-order tensors in a generally covariant way. Key to 
the problem is the construction of an averaged metric that describes the corresponding 
smoothed manifold. Such a metric would allow for the classification of manifolds, but an 
immediate focus are applications in astrophysics and cosmology. In particular, knowledge 
of the averaged metric gives us knowledge of the averaged causal structure of a space- 
time. In a cosmological setting, it would give us the paths of hght propagation through 
the averaged spacetime and as all our observations are necessarily of light this must be 
clearly understood. More fundamentally, the ability to take an averaged metric will enable 
us to construct a physically meaningful metric for a system of bodies. While obviously 
key, there is more to general relativity than the metric. Other important tensors are the 
stress-energy tensor and the Einstein tensor linking the dynamics of spacetime to the en- 
ergy and momentum content of the universe. Neither of these can in general be averaged 
in a physically-meaningful way, but the average dynamics of a domain will be governed by 
Einstein's equation averaged across that domain. 

This thesis focuses on these different aspects of the averaging problem, concentrating first 
on the averaging of the Einstein tensor in a simple, cosmological setting. In chapter |2] we 
give an introduction to cosmological issues that are connected to the "dark energy" prob- 
lem and specifically the attempts to address this with cosmological backreaction. Modern 
cosmology is built on the Friedmann-Lemaitre-Robertson- Walker (FLRW) metric, which 
is maximally-symmetric on hypersurfaces of constant time. The extreme isotropy of the 
Cosmic Microwave Background and the "cosmological principle" together imply that the 
universe is both homogeneous and isotropic to a high degree of precision. This makes 
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the assumption of an FLRW background reasonable. We have by now overwhelming evi- 
dence that when the observations of distant supernovae (of type la) combined with those 
of CMB and LSS are interpreted in this model, the universe must be undergoing a phase 
of accelerated expansion. This was first detected in 1998 [S21 El] and has been refined by 
different studies since, examples being [83lll06j . Combined with observations of the cosmic 
microwave background [921 193] and the large-scale structure of the universe [231 1108[ l22l 173] 
this has lead to the ACDM model of the universe which is in good agreement with all con- 
firmed observational data. In this model the universe is described as FLRW on large scales, 
almost flat and composed of only around 4% standard matter and about 20% "cold dark 
matter" which interacts only through gravity. The acceleration can only be generated in 
this model if 76% of the energy density in the universe is made of the so-called dark energy. 

Probably the most important task of modern cosmology is to understand the nature of 
this dark energy, which may also have implications on particle physics. The most com- 
mon source suggested for the "dark energy" is Einstein's cosmological constant A. This is 
usually explained with the energy density of the vacuum polarisation, but estimates show 
that it should then be many orders of magnitude larger. The most common alternatives 
are scalar fields minimally-coupled to gravity [71] I102[ [T9] which in particular cases can 
act in the present universe to mimic a cosmological constant. While successful, it would 
perhaps be preferable to test models which do not require new physics. 

One physical possibility to explain the origin of dark energy is based on the non-linearity 
of the Einstein tensor and the consequence that the averaged Einstein tensor is not the 
same as the Einstein tensor of the averaged metric. The difference between these two 
tensors in a cosmological setting must necessarily yield corrections to Einstein's equation 
which is known as "cosmological backreaction" . This was first recognized in the 60's by 
Shirokov and Fisher [88]. After this the idea remained generally unconsidered until the 
80's [25] when Ellis brought many questions about averaging in relativity and cosmology 
to a wider audience. The discovery of dark energy has however created increased interest 
in this problem and recent efforts aim to explain dark energy as a pure backreaction ef- 
fect. However, most of the works done in this field are based on averaging schemes that 
either directly or indirectly depend on the used coordinate system raising the question of 
whether or not they are physical. The most common formalism used at present is the 
so-called Buchert formalism introduced in 2000 [T3l [Tl|. This separates cosmology into 
spacelike hypersurfaces and uses a three-dimensional average inspired by Newtonian grav- 
ity. Alternatives with improved averaging schemes include a generally-covariant scheme 
introduced by Zalaletdinov [1111 1112] . Whether the backreaction can form the observed 
dark energy is still subject of debate [1031 EHl EH |52l |53l |80] but there are few quantitative 
predictions from realistic cosmological models. 

In chapter [3] we address this issue. It has been recognised for some time that the back- 
reaction from linear perturbations on an FLRW universe must be small [1031 176] . This 
has, however, never been quantified. We set up the Buchert formalism in "Newtonian 
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gauge" to avoid any problems from gauge modes and then use a "Boltzmann integrator" 
to calculate the linear perturbations in particular cosmological models. Doing this for the 
ACDM model yields a backreaction that is far too small to drive the acceleration of the 
universe, and of the order of the rough predictions. For the alternative "Einstein-de Sitter" 
cosmology without dark energy the backreaction is larger but still negligible. The effective 
fluid also mimics a dark matter instead of a dark energy. The analysis can be extended to 
smaller scales using the "halo model" of structure formation. This is a phenomenological 
model based on numerical simulations and observations of large-scale structure. Although 
the approximation breaks down at small scales it provides an estimate for the corrections 
from "quasilinear" perturbations. We find that these are of the same order as the linear 
effect. The corrections also still act as a dark matter and not a dark energy. 

Chapter |4] forms the core of this thesis because it contains the formulation of a generally 
covariant averaging process which can be used to construct an averaged metric within the 
framework of general relativity. The chapter starts with an introduction to the constituents 
needed for this process and their modes of operation. First we present a Wegner- Wilson 
line operator reformulated for general relativity, which is used to parallel transport tensor 
quantities from one spacetime point to another along the geodesic connecting them. The 
metric can only be averaged if it is first separated across a tetrad field. We present a vari- 
ational formalism to select a particular tetrad from among the several options determined 
by the metric. This is done with a proposed Lagrangian based on the covariant derivative 
of the tetrad. We refer to the field selected as the "maximally smooth" tetrad field. We 
end by presenting the actual averaging procedure itself along with an explanation of its 
functionality. 

Any reasonable averaging process must leave a space of constant curvature invariant. These 
spaces are also the clearest aid to visualizing the process and have an immediate cosmolog- 
ical relevence as they form the basis of the FLRW metric. We consider the three possible 
constant curvature spaces in two dimensions in chapter [5j After employing stereographic 
projection onto a plane, the different constituents of the averaging process are computed 
in terms of the coordinates of the projection plane. It is then shown that these spaces are 
invariant under averaging, as should be expected. 

In chapter [6] we turn to a more physically-interesting case. We take a two-sphere whose 
metric is slightly perturbed and calculate the constituents of the averaging process for a 
general perturbation function. This involves the solution of a geodesic equation that is 
parametrized by a parameter different from the arc length to account for the changes due 
to the perturbation. The Lagrangian formulation for the maximally smooth dyad field 
leads to a differential equation of Neumann type, which cannot be solved analytically. In 
principle the solution to this equation, which determines the dyads, can be found numeri- 
cally. 

The analysis of chapter [6] is extended to the three-sphere in chapter [7} It is shown that 
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the averaging process leaves the smooth three-sphere invariant as should be expected. We 
then consider the case of a linearly-perturbed three-sphere which is directly relevent to 
cosmology and to closed linearly-perturbed FLRW models in particular. Varying the La- 
grangian to find the maximally smooth triad results in differential equations analogous to 
those in chapter [6] but significantly more complicated. The equations are intractable both 
analytically and numerically. However, in more suitable coordinates these equations would 
be better controlled and they act as a proof of concept for the averaging process applied 
to immediately relevent models. 

We solve the differential equation for the maximally smooth dyad field of a perturbed 
two-sphere in chapter [8} with the aid of the simulation toolkit Gascoigne. The particular 
model chosen is that of a gaussian perturbation. It is shown that the averaging process 
does not yield the expected result. 

In chapter [9] we analyse the shortcomings of the averaging method and present a method to 
overcoming the problems with the formulation of an improved Lagrangian. The proposed 
Lagrangian specifies a tetrad that is rotated in a manner taking into account the curvature 
of the manifold characterised by the Ricci scalar. The strength of the average is therefore 
directly linked to the fiuctuations in the manifold. 



In chapter 10 we summarize the investigations and the open issues and give a brief outlook 
on future work. 



CHAPTER 2 I 

THE BACKREACTION PROBLEM 



2.1 The Standard Model of Cosmology 

2.1.1 The Friedmann-LemaTtre-Robertson-Walker Model 

The cosmological principle states that the universe is homogeneous and isotropic. On 
account of this, Friedmann-Lemaitre-Robertson- Walker (FLRW) model has become the 
basis of the standard "big bang" model of cosmology. It is based on an exact solution to 
Einstein's equation (A. 17), 

Gf^u = Rf,u -^gfiuR = SnGTf,,, + Ag^,,,, (2.1) 
with a metric that is homogeneous and isotropic on constant time hypersurfaces, 

ds^ = df - a\t) ( + r^de^ + sin^ OdA . (2.2) 
\ 1 — kr^ J 

By a suitable choice of units for r the constant k can be chosen to have the value +1, or -1, 
indicating constant time hypersurfaces of positive, zero or negative curvature, respectively. 
The only dynamical variable in this metric is the time-dependent scale factor a{t). The 
universe is assumed to be filled with a perfect fluid with the energy-momentum tensor 

T^^ = (p + p)u''u'' - pg^". (2.3) 

Here, p denotes the pressure, p the energy density, and the velocity of the fluid in the 
comoving frame. The terms "pressure" and "density" must be interpreted in a very general 
manner and correspond to the usual sense only in case of a perfect fluid. A negative value 
for the pressure, for example, is not unusual in cosmology. 
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Inserting the energy-momentum tensor (2.3) and the FLRW metric (2.2) into Einstein's 



equations (2.1) we find the Friedmann equations, 



a 
a 



k SttG a 

3 3 

47rG , , A 
'^(p + 3p) + -. 



(2.4) 
(2.5) 



The latter is sometimes referred to as the Raychaudhuri equation and demonstrates the 
fundamentally attractive nature of gravity, and the repulsive effect of a positive cosmo- 
logical constant A. We can also see that, if we interpret the cosmological constant as a 
uniformly distributed fiuid, its energy density and pressure are defined as 



A 



Pa 



SttG 



and 



A 



PA = 



8nG 



(2.6) 



While a universe containing matter and a cosmological constant fits the observations, its 
tiny energy scale has lead many authors to set A = and instead employ a time-varying 
fiuid known as "dark energy" to mimic a cosmological constant. The density and pressure 
of the dark energy are still often denoted by phit) and PA(t). The most popular model is 
known as quintessence |102l [7T| 119] and consists of a minimally-coupled scalar field (f with 
a potential V{ip). The energy density and pressure for such a field are pa = + V{ip) 
and Pa = ^ip^ — V{ip). For a sufficiently fiat potential, <^ V{'^) we find pa — —PA, 
mimicking a cosmological constant. 



momentum conservation (A. 18) 



The Raychaudhuri equation (2.5) can be derived from equation (2.4) and the energy- 



p + 3-(p + p) = 0. 
a 



(2.7) 



Therefore, equations (2.4) and (2.7) form an equivalent, commonly-used pair of dynamical 



equations to the Friedmann equations. We furthermore introduce the Hubble rate. 

Hit) ^ ^, 



(2.^ 



With this definition the usual set of dynamical equations for the universe can be expressed 
in the following way. 



k 



i-kG a 



3 3 
p + 3H{p + p) = 0. 

Furthermore, the Hubble rate allows us to define the critical density 

3H^ 



(2.9) 



Pcit) = 



SttG' 



(2.10) 
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for which the universe is exactly flat A; = 0. A larger energy density would imply a closed 
universe with k = 1, and a smaller value would imply an open universe with k = —1. 
It is customary to express the matter density Pm(t) (including baryonic matter and dark 
matter), the radiation density Prad(^) (photons and neutrinos) and the effective density of 
dark energy pA{t) as fractions of the critical density: 

n.w.^. (2,11) 



Sometimes, even the geometry term fc/a^ in the Friedmann equation (2.4) is assigned an 
energy density contribution, which is defined as 

Then, we find immediately that FLRW universes obey 

ftm + ^r^d + ^A + ^k = l. (2.13) 

The set of equations so far is not closed. To solve for a{t) we must further specify an 
equation of state, 

p{t)=w{t)p{t), (2.14) 

for the dominant component of the energy density. For example, if the universe is fiat 
{k = 0) and radiation dominated {w = 1/3) we find a oc t^/^. In a matter dominated fiat 
universe, on the other hand, we can neglect the pressure, such that w = 0. The scale factor 
then behaves as a oc t^/^. If the equation of state is w < —1/3, the universe accelerates, as 
can be seen from the Raychaudhuri equation. For a cosmological constant, w = —1 and 
we find a oc c^*. 



As the universe expands, the other galaxies are moving away from us and the light they 
emit gets redshifted on the way. This can be used to determine their distance from our 
galaxy. Assuming that the laws of atomic and molecular physics haven't changed since the 
emission of the light of wavelength Acmit, the redshift z in the spectral lines is defined from 
the observed wavelength Aobs as 

1 + ;^=^. (2.15) 

^emit 

In terms of the scale factor a{t) the redshift is given by 



l + z = 




(2.16) 



where ao is value of the scale factor at present t = t^. Redshift surveys such as the Sloan 
Digital Sky Survey (SDSS) tlQSj and the 2dF Galaxy Redshift Survey (2dFGRS) [23] 
measure the redshifts of hundreds of thousands of galaxies and produce maps of the galaxy 
distribution to determine the large-scale structure (LSS) of the universe. The 2dFGRS 
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measured the spectra for 245,591 objects and probed the structure in the local universe 



out to z ^ 0.3 (see Figure 2.1). The SDSS is ongoing [2j and currently covers 585,719 
galaxies up to 2; ~ 1 with a mean at 2; ~ 0.3 and 103,647 quasars up to 2; ~ 3 with a mean 
at z ~ 1.5. For a better comparison to observations, cosmological parameters are usually 
expressed as functions of redshift z instead of the time coordinate t. In the local universe 
the expansion, commonly called the Hubble Flow, is smooth, and for small distances d and 
small radial velocities v of the galaxies the Hubble law 

V = Hod (2.17) 

holds. This can be used to measure the Hubble constant Hq, the present value of the 
Hubble rate. Combining a number of different observation methods, the Hubble Telescope 
Space Key Project [30] found 

Ho = 100 h km s~^Mpc"^ with h = 0.72 ± 0.08, (2.18) 

giving a value of to = = 13.7 Gyr for the age of a flat universe with the fractional 
energy densities = 0.7 and = 0.3 at the present epoch. A precise measurement of 
Ho is important for the determination of many cosmological parameters including in turn 
the current values of the fractional energy density contributions, Q^, fl'^^^^ and The 
measurement of the current deceleration paramter, 

is more involved and requires the comparison of the apparent magnitudes of high-redshift 
supernovae la with those of lower redshift. After calibration the supernovae of type la have 
a nearly uniform maximum luminosity and can therefore be used as "standard candles" 
for distance measurements. 

Another important quantity is conformal time, the distance photons could have traveled 
since t = 0: 

dt' r da' 1 



Jo a{t') Jo a' a'H{a') ' ^^'^^^ 

Just like t, z and a, the conformal time r serves as a reference parameter when discussing 
the evolution of the universe. Furthermore, particles which are separated by distances 
larger than r have never been in causal contact with each other. Therefore, it is also 
called the comoving horizon, or more precisely the comoving particle horizon. For a fiat 
matter-dominated universe with a oc 

^2/3 

it is identical to the comoving Hubble radius or 

Hubble length, 

Rh = (2.21) 
aH 

Particles separated further than this are currently not in causal contact. Since the expan- 
sion of the real universe has changed over the time, the comoving horizon and the comoving 
Hubble length are actually not identical, but we can estimate the particle horizon to be of 
the order of the Hubble length today. 
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Figure 2.1: Map of the galaxy distribution from the 2dF Galaxy Redshift Survey [23] 



2.1.2 Cosmological Perturbation Theory 



As one can see from Figure 2.1, the observed matter distribution in the universe is nei- 
ther homogeneous nor isotropic on small scales. It is dominated by giant bubble-like voids 
separated by sheets and filaments of galaxies, with superclusters appearing as occasional 
relatively dense nodes. Quantitative estimates employing the fractal dimension of the data 
suggest that on scales above ~ 70Mpc//i the large scale structure becomes homogeneous. 
See for example [98j for a recent analysis of the SDSS data, in broad agreement with older 
studies [H I107j . However, it has been suggested that a large cold spot in the WMAP data 
is generated by the largest observed void with a radius of ~ 140 Mpc [H]. Regardless of 
this, we can say that structure beneath a scale of 70Mpc is largely inhomogeneous. This 
structure must have grown from small initial fluctuations via gravitational instability. 



Since the universe has been expanding from a much smaller state, the very early universe 
was made up of a hot plasma. This plasma was opaque since the photons were interacting 
constantly with the plasma through Thomson scattering. As the universe expanded the 
plasma cooled down adiabatically until the electrons and protons recombined to hydro- 
gen atoms. The photons scattered off the neutral atoms and continued their path freely 
through the now transparent universe. These photons from the surface of last scatter- 
ing (2;= 1100) can be observed today as cosmic microwave background (CMB). Due to 
the expansion of the universe their wavelength has been stretched and they now form a 
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black-body spectrum of 2.726 K. The primeval plasma, tightly-coupled to the photons by 
Thomson scattering, underwent acoustic oscillations driven by gravity and radiation pres- 
sure which are imprinted on the CMB as temperature anisotropies in the form of acoustic 
peaks in the angular power spectrum. These are measured to high precision by CMB 
experiments such as the Wilkinson Microwave Anisotropy Probe (WMAP) [921 [93], [9] and 
ACBAR [81j and contain much information about cosmology at the time of recombination. 
In particular, they tell us about the corresponding fluctuations in the matter density which 
were the seeds of the large scale structure we observe today. The inflationary scenario links 
these seeds in turn to quantum fluctuations in the very early universe and provides a rel- 
atively natural justification for the amplitude and near scale-independent power spectrum 
the observations imply for the primordial perturbations. However, the density fluctuations 
remain small enough to safely assume that the universe was very nearly of FLRW type 
in the early universe. In the early universe, therefore, we can apply linear perturbation 
theory (see IM] for a review) to find the dynamical equations for the inhomogeneities. 

The full metric describing the mildly inhomogeneous universe is split into the background 
metric plus perturbations, 

9,^ = 9^'i! + Sg,u. (2.22) 
The line element for the background FLRW metric in Cartesian coordinates is 

ds^ = gf}{x)dx^'dx'' = a\T){dT^ - jijdx'dx^), (2.23) 

where dr = a^^dt is the conformal time and 

7., = (l + \k{^' + y' + ^')) ' • (2-24) 

The metric perturbations may be decomposed into scalar, vector and tensor perturbations 
according to the way they transform under coordinate transformations on the hypersurfaces 
of constant conformal time. The decomposition theorem states that this decomposition is 
preserved by Einstein's equation to first order, and that each of these types of perturbations 
evolves independently. The most general form of the line element for the background and 
perturbations is 

ds'^ = a2(r) ((1 + 2<I))dr2 - 2{V,B + Si)dx'dT 

- ((1 - 2^)7ij + 2V{DjE + {VjFi + ViFj) + Kj) dx'dx^) . ^'^''^^^ 

The covariant derivative on the three-dimensional background hypersurface with respect to 
some coordinate i is denoted by Vi. Since vector perturbations {Si, Fi) decay kinematically 
in an expanding universe and tensor perturbations {hij) grow at a slower rate than scalar 
perturbations (0, ip, B, E) and represent gravitational waves, which do not couple to energy 
density and pressure inhomogeneities, we concentrate on scalar perturbations. Therefore, 
we start from the most general form of the scalar metric perturbations 

^»-aV)( o.,„. "''Jn'^ J . (2.26) 
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There is a freedom in the choice of coordinate system with which to study the space- 
time. The most common gauges are the synchronous gauge, where one sets = 0, I? = 
and the conformal Newtonian, or longitudinal, gauge with = 0, S = 0. The advan- 
tage of synchronous gauge is that it closely resembles Minkowski space with purely spatial 
perturbations, because we have chosen the hypersurfaces of freely falling observers. More- 
over, there is an unambigious time parameter which eases numerical implementation. The 
drawback of synchronous gauge is that one freedom remains unfixed and this introduces 
unphysical artifiacts, or gauge modes, which must be dealt with carefully. (See [TSl E] for 
more discussion on this.) A more convenient choice is the Newtonian gauge, where one 
chooses hypersurfaces of isotropic expansion rate perturbation. If we restrict our attention 
to spacetimes without anisotropic stress, which by the Einstein equations implies (j) = ip, 
the line element 

ds^ = a'^{r]){{l + 2(t))drj^ - (1 - 2(j))-fijdx'dx^) , (2.27) 

is the same as the Newtonian approximation to GR, whence may be interpreted as a 
generalization of the Newtonian potential. 

The best way to get around the gauge problem is to express all quantities in a gauge- 
invariant way. This approach was pioneered by Bardeen [6]. Consider infinitesimal coor- 
dinate transformations 

^ = + (2.28) 

with the four- vector C,^ = Then one can show that only and the function ^, 

which is defined as the solution to 

V'V,^ = V,e, (2.29) 

preserve the scalar nature of the metric fluctuations and therefore contribute to their 
transformation. The most general gauge transformation for scalar perturbations in terms 
of these functions is 

T^f = T + ^^{x) and x' ^ x' = x' + YWj^{x). (2.30) 

Denoting the derivative of a function / with respect to the conformal time r by /', the 
change in the scalar perturbation functions 0, ip, B and E under the gauge transformation 
can be expressed as 

4> = <j>--e-e. ^ = ^+-e°, B = B+e-^\ e = e-^. (2.31) 

a a 

Although these functions are not gauge- invariant, certain combinations of them are, and 
the simplest gauge-invariant combinations are the Bardeen potentials 

^ = cf) +-{{B - E')a)' and ^--{B-E'). (2.32) 



In Newtonian gauge (2.27) the Bardeen potentials $ and \l/ coincide with the scalar per- 
turbations and ^, respectively. 



12 



2. The Backreaction Problem 



Depending on the choice of problem one wishes to study, the corresponding energy mo- 
mentum tensor can now be decomposed in a similar way. At very high energies it is no 
longer reasonable to believe that a hydrodynamical description of matter will be valid. 
The matter will then be described in terms of fields as is for example done for inflationary 
cosmology. With a different motivation the same formalism is used in the quintessence 
model of dark energy. The theory of a scalar field minimally coupled to gravity is given by 
the action 

^-D'^'^Do,'^ - ^(v?) ) y/^d^x, (2.33) 



S 



where Z)^ denotes the covariant derivative (A. 11) and V{ip) is the potential of the scalar 



field. The corresponding energy-momentum tensor is 



(2.34) 



To be consistent with the metric (2.26) the scalar field (p{x) must also be approximately 



homogeneous and can be decomposed into 

(p{x) = Tp{t) + 5(p{x) 



(2.35) 



where an overbar represents the homogeneous component. The energy- momentum tensor 
can thus also be decomposed into background and perturbed parts, 



= t:+ 5t^, 



(2.36) 



where STf^ is linear in matter and metric perturbations 6ip and SgaB- Inserting this and 



the metric (|2.26|) into (|2.35|) gives the background energy-momentum tensor 

1 



2a2 
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Vijp) = T. = 0, T, 



and the first-order perturbation 



(2.37) 



(2.38) 



Here, the comma with the space index means differentiation with respect to the correspond- 
ing coordinate and = dV/dip. This energy-momentum tensor can also be expressed in 



a completely gauge-invariant form if one wishes. From (2.37) we can see that if the field 
is slowly- rolling so that ip^ <C a?V{Tp), = —p^ and the scalar field mimics a cosmolog- 
ical constant. In the opposite hmit = p^ and the scalar field acts as a "super-stiff" fluid. 



For the hydrodynamical description of matter we can linearise the density by p = p(l -|- 5) 
for a dimensionless density perturbation 6, and write the barotropic pressure as p = wp. 
Then one can find the linearised energy-momentum tensor 

= p(l + S),T^ = p(l + w)v„ T] = p{w + clS)Si (2.39) 
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with the background sound speed = dp/ dp. 

Both the perturbed metric and energy-momentum tensor are then inserted into Einstein's 
equation and the resulting equations are combined in such a way, that they form a gauge- 
invariant set of equations for the Bardeen potentials. These equations can be split into 
equations for the various different forms of energy. The conservation equations for fluids 
interacting only through gravity, such as neutrinos and cold dark matter, can be then sep- 
arated from the photon/baryon system since in the concordance model the only important 
interaction is Thomson scattering of photons by electrons. One then solves the Boltzmann 
equation for the photons coupled to a perfect baryon fluid. (See for example [55] for details.) 

The perturbations are best described in terms of harmonic modes, characterized by the 
magnitude of the wavevector k = y/Wkl. In fiat space, these take the form of Fourier 
modes. 



A(x) = / A(k)e-^-^. (2.40) 

In linear perturbation theory, the different Fourier modes evolve uncorrelated from one- 
another. This implies that their initial statistical distribution is conserved, and one may 
link their power spectrum at late times directly to the initial power spectrum. A given 
mode initially has the power spectrum 

{A*{k)A{k')) = j^V4k)6{k - k'), (2.41) 

where V^{k) is the primordial power spectrum of the metric fluctuations and the Dirac delta 
function enforces statistical homogeneity. According to the theory of inflation the power 
spectrum is nearly a scale-invariant Harrison-Zel'dovich-Peebles spectrum and follows a 
power law 

V^{k) oc (2.42) 
where is called the spectral index. For adiabatic fluctuations ~ 1. 

Numerical CMB packages such cmbeasy [23] , CAMB [51] or cmbfast [87] solve the Einstein 
equations, along with Boltzmann equations for the photons and neutrinos, and perfect 
fluid equations for the baryons and CDM. (For this reason they are often referred to as 
"Boltzmann codes".) The solutions to these equations then tell us how the amplitude of 
each wavemode changes between the earliest times and recombination and can be used to 



wrap (2.41) onto the CMB by integrating 

oc y V^{k) \/\Ti{kf d{\n{k)) (2.43) 

where A^-^ is the photon "brightness function" expanded across the spherical harmonics and 
/ is the multipole. Observations of the CMB combined with studies of the large-scale struc- 
ture and Supernova type la data remain consistent with an initial, Gaussian distribution of 
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Figure 2.2: The observed matter power spectrum from the SDSS Luminous Red Galaxy 
survey (figure taken from [73]). The sohd fine is the ACDM prediction. 



adiabatic perturbations with a spectral index very close to Harrison- Zel'dovich-Peebles [93] . 



CAMB, cmbfast and cmbeasy also evolve the matter perturbations onward to the present 
day. The quantity is known as the matter power spectrum P{k) and quantifies the 



nature of the clustering of baryons and CDM (see Figure 2.2). The solid line in the figure 
is the prediction from the linear concordance ACDM model, while the data points are from 
the SDSS Luminous Red Galaxy survey [73]. For large scales and a Harrison- Zel'dovich 
primordial power spectrum, P{k) cx k, while on small scales P{k) oc /c~^. The turnover 
refiects the time at which the energy densities matter and radiation were approximately 
equivalent and can be seen to be at 2; ~ 4500. It is important to notice that modes outside 
of the Hubble horizon {k < 1/aH) are gauge-dependent. While in the synchronous gauge 
they retain a power spectrum P{k) oc k, in the Newtonian gauge they grow as one tends 
to larger scales. Within the horizon the two gauges agree, as they should for observable 
quantities. 



Imprinted on the matter power spectrum are ripples entirely analogous to the oscillations 
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in the CMB. These are the imprint of the oscillations the baryons underwent prior to de- 
coupling from the photons and their wavelength provides a "standard ruler" with which to 
probe the expansion of the universe since the formation of the CMB (see for instance [72] 
and its references). Above k^]^ ?s 0.1 Mpc//i, the matter power spectrum is dominated by 
nonlinear effects and difficult to evaluate. This is clearly evidenced in Figure 2^ with the 
data points increasingly diverging from the linear prediction on smaller scales. It may be 
determined for specific models from A^-body simulations |9lj, evaluated from exact inho- 
mogeneous solutions to Einstein's equations |17], or evaluated phenomenologically using 
the halo model of galaxy clustering as in, for example, the Halofit code 



2.1.3 The ACDM model 

From observations of galactic rotation curves, velocity dispersions in galaxy clusters, and 
mass tracing via gravitational lensing, we know that there must be more mass in the uni- 
verse than can be observed. Since this energy component of the universe is not interacting 
with photons it has to be non-baryonic and "dark", but it has to interact gravitationally 
and form structures. This implies that it moves non-relativistically and is "cold" . There- 
fore, this energy component has been given the name "cold dark matter" (CDM). 

The direct evidence for the existence of an additional dark energy component comes from 
distance measurements of type la supernovae [TH [821 |83l 1106] . which indicate that the 
expansion of the universe is currently accelerating. This fact can only be explained in 
terms of a FLRW model if it contains a dark energy component with a negative pressure. 
Recently, Melchiorri et al used combined data sets to constrain the redshift of transition 
from deceleration to acceleration, finding a transition at Zacc = 0.76 ± 0.10 and equality 
between the dark energy and matter at Zeq = 0.40 ± 0.08 [6T] . 

The power-law ACDM or "concordance" model is the simplest known model that is in 
agreement with all astronomical observations, including the CMB, the cosmic expansion 
history determined from supernovae type la data, and the matter power spectrum, mea- 
sured for example by Lyman a Forest, galaxy clustering and motions, gravitational lensing, 
cluster studies, and 21 cm topography. The model assumes the universe to be spatially flat, 
homogeneous and isotropic on large scales. It has a cosmological constant and is composed 
of ordinary matter, radiation, and dark matter. The primordial fluctuations in this model 
are adiabatic, nearly scale-invariant Gaussian random fluctuations. CMB anisotropy mea- 
surements combined with the LSS data and numerical studies (e.g. [91]) indicate that 76% 
of the energy density in the universe is from dark energy, 20% is from dark matter, and 
only 4% is baryonic matter [85] |97] . 

One issue with the ACDM model concerns the CDM and structure formation on small 
scales. The model predicts a significantly higher number of dwarf galaxies in our local 
group than are found [62l |13]. This is called the "missing satellite" or "missing dwarf" 
problem. Although improvements in the observation technology diminish this discrepancy 
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[90] . the problem still remains unsolved. There has also been a long-standing discrepency 
between CMB predictions and LSS observations of as, the normalisation of the matter 
power spectrum, with the CMB tended to suggest lower values than that measured by 
LSS probes. These are now converging (see for example [HI]) but some tension remains. 
Furthermore, the observed void of radius ~ 140 Mpc [HI] is far outside the current expec- 
tations of the ACDM model. 

Other than the problems on small scales the ACDM model has been phenomenologically 
very successful in fitting a wide range of cosmological data. However, it remains purely 
phenomenological and does not provide any explanation of the nature of dark energy and 
dark matter. If dark energy is explained by a cosmological constant a severe fine-tuning 
problem arises. Estimating the Casimir force of vacuum fluctuations in the universe gives 
a value of the order of Mf, with the Planck mass Mp, but the cosmological constant is of 
the order of A 10~^^^Mp. Why the cosmological constant is so extremely tiny, but not 
quite zero, is an unresolved problem known as the "cosmological constant problem". If 
we interpret the cosmological constant as a fluid with a dynamical equation of state, we 
have to explain why this component was negligible for most of the expansion history of the 
universe and became dominant only recently. This problem is known as the "coincidence 
problem" . 

It is important to notice that models without a cosmological constant are consistent with 
the WMAP data alone, provided Hq < 45 km s~^Mpc~^ and the primordial power spec- 
trum is slightly modified [931 IHE]. The flat, matter-dominated case is the so-called Einstein- 
DeSitter (EdS) model. Combined with large-scale structure observations (particularly the 
baryon acoustic oscillations), supernova data and measurements of local dynamics these 
models seem unhkely, but are nevertheless useful toy models to test the data and challenge 
the ACDM model. 



2.2 The Averaging Problem 

Considering that the expansion of the universe has changed from deceleration to accelera- 
tion in the recent past, when the formation of structure has grown increasingly non- linear 
on small scales, it seems sensible to suspect that these effects may be connected. If so, this 
would be an elegant way to solve the coincidence problem. Solving Einstein's equation on 
large scales for a general inhomogeneous universe raises the "averaging problem", which 
we will now discuss in detail. 



The first important indication that the standard model of cosmology is based on an in- 
correct dynamical equation came from Shirokov and Fisher in 1963 ([88j, and reprinted as 
|89j). They realized that in the standard approach an averaged energy-momentum ten- 
sor, corresponding to a continuous matter distribution, is used, while the left-hand side of 
Einstein's equation (2.1) remains unaveraged. Instead, the solution is assumed to be the 
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Einstein tensor for an averaged metric, 

G^A^gpx)) = SttG (T^,) + A {g,,) . (2.44) 

This equation is clearly wrong, since from the nonlinear nature of the Einstein tensor we 
know 

{G,A9px)) ^ G,,{{g,^)). (2.45) 

Shirokov and Fisher suggested introducing a correction term SnGT-j^^ into Einstein's equa- 
tion, 

G,.{{gpx)) = SttG {T,,) + SnGT^^, + A {g,,) , (2.46) 

which depends on the fluctuations of the microscopic gravitational field and is of purely 
geometrical origin. They called it a "polarization term" and remarked that it might be 
interpreted formally as extra terms in the energy and pressure of the cosmic fluid. 

The authors proposed the following averaging process 

where g^y{x) is the exact value of the metric tensor at the point x, and the domain of 
averaging D, containing the points x, is assumed to be large compared to the fluctuations 
and small compared to the scale at which systematic variations of the mass or metric dis- 
tribution are possible. 

Shirokov and Fisher applied this averaging process to the Newtonian approximation and 
estimated the effect of small fluctuations. Their main result was that the "polarization 
term" mimics repulsive gravitation and prevents the Big Bang singularity, although they 
stated that this result was preliminary, since they extended their investigations beyond the 
applicability of perturbation theory. 

The averaging problem remained generally unconsidered until Ellis emphasized its impor- 
tance in great detail in 1984 [25]. He discussed different scales of inhomogeneity and argued, 
that a suitable averaging process has to reproduce the same metric on the largest scales 
both when applied only once or when applied subsequently to the intermediate scales. 
Such a process has to be able to smooth out high-curvature phenomena such as cusps, 
caustics and singularities in a consistent way. Furthermore, the averaged Einstein equa- 
tion has to describe the same dynamical behaviour of the manifold on all scales when one 
additionally includes the correction term. Ellis also raised the "fitting problem" , which he 
discussed in more detail in a subsequent paper with Stoeger |26j. This concerns the prob- 
lem of finding the FLRW model which is the best fit to the inhomogeneous real universe, 
in the same sense that an idealized perfect sphere is fitted to the real earth. This analogy 
is taken further by posing the question which coordinates in the manifolds of different 
averaging scale must be identified with one another, in order to be able to quantify the 
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deviation from the best fit FLRW model and hence to indicate how good the fit actually is. 
Noonan ^\ defined the macroscopic average {Q') of a microscopic quantity Q' by 

j ^/-9{x')d^x' 

Applied to the Newtonian weak-field, slow-motion approximation, the author found that 
the contribution by the small-scale gravitational fields to the macroscopic density and 
energy-momentum tensor are the Newtonian gravitational energy density and energy- 
momentum tensor, respectively. However, g{x') is assumed to be the determinant of the 
macroscopic metric and it is not explained how that metric is obtained. 

Following pioneering work by Brill and Hartle [12], Isaacson defined an effective energy- 
momentum tensor for gravitational waves [39] and showed its gauge-invariance after av- 
eraging for high-frequency waves. The author realized that tensors defined on different 
spacetime points cannot be compared directly, and therefore simple volume averaging de- 
stroys the tensor character of the averaged quantities. Therefore, he used the "bivector 
of geodesic parallel displacement", denoted by V^^ {x,x'). This bivector transforms as a 
vector with respect to coordinate transformations at either x or x' and, assuming that x 
and x' are sufficiently close together to ensure the existence of a unique geodesic of the 
metric gfj,i, between them, given the vector Au' at x' then = Vfj,'^ A^i is the unique vector 
at X which can be obtained by parallel-transporting Ayi from x' to x along the geodesic. 
The average of a tensor T^p is then defined as 

{T^p{x)) = j V/{x,x')V/ix,x')n,yix')f{x,x')d^x', (2.49) 

where f{x, x') is a weighting function which falls smoothly to zero when x and x' differ by 
the averaging radius, and which is furthermore normalized to unity, 

j4 / 



f{x,x')d^x' = 1. (2.50) 

Although Isaacson's procedure looks quite promising at first sight, Stoeger, Helmi and 
Torres [95] realized that it cannot be used to average the metric: the parallel-transported 
metric at x' is identical to the metric at x since the covariant derivative of the metric 
vanishes everywhere, 

V^"' {x,x')Vp^' {x,x')gu'x'{x') = gpp{x). (2.51) 

The metric is hence invariant under Isaacson's averaging procedure, {gp,p) = g^p- Stoeger, 
Helmi and Torres proposed instead another averaging process, 

(g(x)) = !Qi^+ ^'U-9{^ + ^') d'^' (2.52) 
/ \/-9{x + x') d'^x' 
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The authors apphed this process to weak-field and perturbed FLRW models and demon- 
strate approximately tensor character of the averaged quantities to linear order. Unfor- 
tunately, this is not necessarily the case in general and the result is in any event only 
approximate. 

Many of the early approaches to the Averaging Problem employed spatial averaging. 
This was applied to perturbed FLRW models (e.g. pT] [32l W\\ [95]). after implement- 
ing spacetime-slicing (e.g. |3T]) or in the context of perturbed Newtonian cosmologies 
(e.g. [IZ])- A review can be found in |17]. The most promising approach is the Buchert 
formalism and we will therefore discuss it in detail. 

2.3 The Buchert Formalism 

Although a simple form of averaging cannot be used to construct an averaged metric or 
to average over other tensor quantities in a generally covariant way, it can still be used 
to find certain averages of scalars. Scalar fields are defined to be invariant under general 
coordinate transformations and their values on different spacetime points can therefore be 
compared directly. To make use of this fact, we need to decompose Einstein's equation into 
a set of dynamical equations for scalar quantities. This is usually done in the framework 
of the 3+1 formalism of general relativity, which we briefly review here. In line with most 
of the literature in this fleld, we use a metric of signature (— , +, +, +) for the rest of this 
section. 

2.3.1 Einstein's Equation in 3+1 Form 

The 3+1 formalism, sometimes also called the ADM formalism, is a standard method in 
general relativity and an essential ingredient of numerical relativity. Details about this 
formalism can be found for example in [361 I109i llOlt [13]. Provided that the spacetime 
(7Vl,g) is globally hyperbolic, we may foliate M. with a family of spacelike hypersurfaces 
S, such that each hypersurface is a level surface of a smooth scalar fleld t, which we will 
later identify with the time coordinate. Then the entire future and past history of the 
universe can be evaluated from the conditions on one of the "Cauchy surfaces" E. The 
normal unit vector n to the surface S being timelike must satisfy 

n-n = -l. (2.53) 

It is necessarily colinear to the metric dual Vt of the gradient 1-form dt and the propor- 
tionality factor is called the lapse function a 

n = -aVt, (2.54) 

where a = (— Vt ■ Vt) ^^'^ = — ((dt, Vt))""*^^^. To ensure that n is future oriented, a 
must be larger than zero and in particular never vanishes for a regular foliation. Since 



20 



2. The Backreaction Problem 



{dt, an) = 1 the hypersurfaces are Lie dragged by the normal evolution vector an. 

Now introduce coordinate systems on each hypersurface E that vary smoothly between 
neighbouring hypersurfaces, such that = [t, ) form well-behaved coordinate 

systems on M.. The 1-form dt is dual to the time vector dt. Hence, just like the normal 
evolution vector, the time vector Lie drags the hypersurfaces. In general these vectors do 
not coincide. The difference defines the shift vector (3 

dt = an + (3. (2.55) 

The shift vector is purely spatial since {dt,f3) = {dt,dt) — (dt, an) = 0. In terms of 
the natural b 
expressed as 



the natural basis of the introduced coordinates x'^ the normal unit vector can then be 



= -(l,-p'). (2.56) 
a 

With the projection operator into S 

hf,u = gf,u + n^n^ (2.57) 
the components of the 3-metric hij become 

hij = g^^uKh'y (2.58) 
From here one may find that the line element in these coordinates is 

ds^ = (-a^ + PiP')dt^ + 2Pidtdx' + hijdx'dx^, (2.59) 

where Pi = hijP^ . 

The projection of the gradient of the unit normal vector defines the extrinsic curvature 
tensor 

Kij = -KV^Hj = -ViUj. (2.60) 

Using Frobenius' theorem we find that the normal unit vector is hypersurface orthogonal 
if and only if 

ny^NvUx] = 0, (2.61) 

where square brackets denote antisymmetrization on the enclosed indices. Therefore, we 
can express the extrinsic curvature tensor as the Lie derivative (A. 13) 

K,, = -"-LnK, (2.62) 

and the evolution of the 3-metric hij is given by the Lie derivative along the normal evo- 
lution vector an. 
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We can use the freedom in the choice of coordinates to set the shift vector = 0. We 
further identify the normal unit vector with the fluid four- velocity u'^, thereby identifying 
the hypersurfaces as the rest-space of the fluid. The line element in these coordinates is 
given by 

rfs2 ^ _c^2^t2 _^ hijdx'dx^, (2.63) 
and the four-velocity of the observers becomes 

n^=Q,oV n^ = (-a,0). (2.64) 



Using (2.55) and (2.62) the evolution equation then becomes 



aCnhij = Cdthij - ^ishij = diKj = ahij = -2aKij, (2.65) 

where an overdot denotes a total derivative with respect to the proper time given by 
dr = adt, T>i is the covariant derivative on the 3-surface. 

Now we perform a 3+1 decomposition of the stress-energy tensor 

T^u = pn^riy +ph^u p = T^^n^rf , phij = Tij (2.66) 

and project the Einstein equations onto the hypersurface and along its normal. This results 
in the Hamilton constraint equation 

TZ + K^- K}Ki = IQnGp, (2.67) 
where TZ is the Ricci scalar on S and K = Kl; the momentum constraint equation, 

Vj {K'^ - h'^K) = 0, (2.68) 
and, with the "acceleration" , defined as 

= n^V^n^" = h'', = h'^^a^, (2.69) 

which describes the deviations from a geodesic flow, the evolution equation for the extrinsic 
curvature tensor can be expressed as 

kij = TZij + KKij - AT^GKj (p-p)- (Vja' + a'aj). (2.70) 

Now we can introduce the expansion tensor Gjj = —Kij and decompose it in terms of 
kinematical quantities and their scalar invariants, 

®- - = V.-, = + a., + lh,,e. (2.71) 

These are the trace-free antisymmetric "rotation tensor", 

uJij = V[inj], (2.72) 
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the trace-free symmetric "shear tensor'' 



(2.73) 



where round brackets denote symmetrization on the enclosed indices, and the "rate of 
expansion" , which is the trace of the expansion tensor, 



(2.74) 



Using once more Frobenius' theorem (2.61), we find that the manifold representing the 
universe is globally hyperbolic if and only if it is irrotational, and hence uoij = 0. 



With the "rate of shear" a, defined by 



a 



we write down two out of the three scalar invariants of the expansion tensor 

-K = e and 
The four-divergence of the acceleration field is 



- (K^ - K%K\) = -e^ - a' 
2 V ^ ?, 



A = V^a^" = Via' + aia\ 



Inserting this and (2.67) into the trace of (2.70) gives the Raychaudhuri equation, 

1 



9 



-9' - 2a' - 47rG(p + 3p) - A. 



(2.75) 



(2.76) 



(2.77) 



(2.78) 



2.3.2 The Buchert Equations 

Having decomposed Einstein's equations into a system of equations for scalar quantities, 
we now want to find their averaged form. Following Buchert ITi] , we define the volume 
of a domain "D C S of the hypersurface by 



V 



h(ryi. 



(2.79) 



V 



Then we define the average of a scalar field A in this domain as 



Ay/h£x. 



(2.80) 



V 



Since dt and (i^x commute, we can furthermore introduce the effective Hubble rate by 
^dtax) dtV 1 d 



3H: 



V 



V 



Vdt 



— 1 r 1 ■ ■ ■ 
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Using (2.76) and (2.65) we thus have 

^dtav 



3H' 



V 



(aK) = {a9) 



9). 



(2.81) 



Hence, the effective Hubble rate in a domain is given by the average of the scaled expansion 
9 = a9. We can also quickly find the commutator between time and space derivatives for 
any scalar field A, 

'^"^^ -i-{A) = (A9)-{A)(9). (2.82) 



dt 



dt 



The averaged Raychaudhuri equation (2.78) for the scaled densities p = a p and p = a p 
is 



+ 47rG (p + 3p) = Qv + Vv, 



(2.83) 



and with the scaled spatial Ricci scalar TZ = a^TZ the averaged Hamilton constraint (2.67) 
becomes 



Q 



V 



n 



(2.84) 



Here, we have introduced the domain dependent "backreaction terms": the "kinematical 
backreaction" , 

Qv^l(i9-(9))A-2{a'), (2.85) 



with the scaled shear scalar a = aa, and the "dynamical backreaction" , 

Vv = (a) + (a9^ 



(2.86) 



with the scaled acceleration divergence A = a A. To ensure that (2.84) is an integral of 
(2.83), we additionally need the integrability condition 

dtQv+6Hr,Qv+dt {n)+2H-D I1i)+AH-dPt,-IQitG (dt {p)+3Ht, {p+p)) = 0. (2.87) 



Having now an averaged system of Einstein's equation in scalar form, we can fit it to a flat 
FLRW model as follows. Define the effective energy density and pressure as 



pes = (p) 



Q 



V 



{n 



Pes = (P) 



Q 



V 



167rG 167rG 



V 



167rG 487rG 127rG 



Then the averaged equations (2.83) and (2.84) can be expressed as 

d Qj 

3—+471G (Peff + 3peff) = 0, 



6H. 



V 



IQnGp, 



'eS 



0, 



(2.88) 
(2.89) 

(2.90) 
(2.91) 
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and the integrability condition (2.87) has exactly the form of a conservation equation law 



dtPcS + 3Hv (pcflf + Pcff) = 0. 



(2.92) 



This set of equations is known as "Buchert equations". They show that, despite the 
non-commutativity of time and space derivatives, the averaged variables obey the same 
equations as the variables in the flat FLRW model. Therefore, in the limit of vanishing do- 



main of averaging the equations (2.90), (2.91) and (2.92) smoothly reduce to the equations 



(2.4), (2.5) and (2.7), respectively. However, as pointed out by Buchert (e.g. in ^T5\ 116]). 
there are correction terms to the fiat FLRW variables that counteract gravity. As we can 
see from the definition of the effective energy density (2.88), the first term in the kine- 



matical backreaction (2.85) is positive and hence reduces the value of the energy density. 



This stems from the fact that an average correlates the local fluctuations, which then act 
in the sense of a global "kinematical pressure". Therefore, one cannot exclude large-scale 
effects from averaging inhomogeneities with the argument that the perturbations are of 
small amplitude and act gravitationally. 



The Buchert equations reduce the solution of the averaging problem for scalars to finding 
an effective equation of state that relates the effective densities, 

PeS 

Wes = —. (2.93) 

If, however, we consider the backreaction terms independently as an effective fluid, it has 
an equation of state 



I (n) - 4:Vv - 3Qv 



n) + Q 



(2.94) 



V 



We can then see from this result that should these terms dominate the matter contribution 
with Qti < —Vv and hence w-d < — |, the modifications act to accelerate the averaged 
scale factor. 



The statement that the kinematical backreaction term Qx> in Buchert 's equations can 
play the role of a cosmological constant has been confirmed by many authors (see e.g. 
[HI HHl [531 HZl HH], for a review, see [I5]). A promising approach has been suggested by 
Wiltshire |105[ I104j . He introduced two scales, one that belongs to "finite infinity" re- 
gions, and one that describes the scale factor of negatively curved voids. Then the Buchert 
formalism is applied and the average spherically symmetric geometry is reconstructed in 
terms of a spatially averaged scale factor and time. According to the author, this "fractal 
bubble" model might simulateously resolve key anomalies relating to primordial lithium 
abundances, CMB ellipticity, the Hubble bubble feature and the expansion age. The key to 
the latter is the fact that within this model an observer in a galaxy measures a significantly 
older expansion age of the universe than that of a ACDM model, giving a volume average 
age of the universe of 18.6Gyr. 



2-4 The Lemaitre-Tolmann-Bondi and Swiss Cheese Models 
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2.4 The LemaTtre-Tolmann-Bondi and Swiss Cheese 
Models 

The Lemaitre-Tolmann-Bondi (LTB) solution is the most general spherically symmetric 
solution to Einstein's equations for dust, 

Gf,u = MGpn^n^. (2.95) 

Because of the spherical symmetry the line element is expressed in spherical coordinates 
plus time, 

ds" = -de + ^'^^'''}^ dr^ + Rit, rfide^ + sin^ ed<p''), (2.96) 
1 + E[r) 

where the functions R{t, r) and E{r) are related to each other and to the energy density 
p{r,t) as follows 

B(«,r)2 = 8,rG-^^ + E(r), 

^^^'""^ R{t,ryR'{t,r)' 

where dots and primes denote derivatives with respect to t and r, respectively. E{r) can 
be interpreted as the total energy per unit mass and M(r) as the mass within the sphere 
of comoving radial coordinate r. 



This model is a useful toy for studying backreaction effects, because the symmetry allows 
quantitative studies without the use of approximations. Therefore, it is probably the most 
studied model in this context. In the earliest work [20] the author used the LTB model 
as a toy model to obtain a reasonable fit to supernovae luminosity densities and showed 
that the acceleration implied by the supernovae data can be explained by a large scale 
inhomogeneity without the need for a cosmo logical constant. A later study ^76j showed 
that backreaction slows down the expansion if measured in terms of the proper time, but 
speeds it up if measured in terms of the energy density or the scale factor. The authors of 
[UU] used a solution with both a region with positive spatial curvature and a region with 
negative spatial curvature. They found that after the region of positive spatial curvature 
begins to re-collapse, the averaged universe starts accelerated expansion. This implies a 
strong coupling between averaged scalar curvature and kinematical backreaction. It is im- 
portant to notice, that such an effect cannot be studied in the LTB model with vanishing 
scalar curvature, since for this model Qx) = [70] . 

However, although there are many studies about the averaged luminosity and angular di- 
ameter distances (see for example [21 [23 ESI ES]) in this model, it remains merely a toy 
model which and unsuitable to describe the real universe, since it violates, for example, 
the cosmological principle. On the other hand, it is a suitable choice to model the observed 
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voids in the universe and is as such implemented in the swiss-cheese modeL 

The swiss-cheese model is a more realistic model which consists of LTB patches embedded 
in a fiat FLRW background containing only matter. In (STf l56] the authors did not analyze 
the averaged domain dynamics in this model but instead the propagation of photons. 
They found that the light-cone average of the density as a function of redshift is affected 
by inhomogeneities in a way that the phenomenological homogeneous model (which is 
identical to the EdS background) behaves as if it has a dark-energy component. The effect 
arises because, as the universe evolves, a photon spends increasingly more time in the 
large voids than in the thin high-density structures. The authors conclude further that 
within their toy model the voids must have a present size of ~ 250 Mpc to be able to 
mimic the ACDM model. In contrast, photon propagation in a swiss-cheese model has 
been independently studied and very small overall effects were found if the observer 
sits outside the void, while they are large if he sits inside. 



2.5 Super-Hubble Fluctuations 

In inflationary cosmology matter is treated as a scalar field in the context of cosmolog- 



ical perturbation theory (which was introduced in section 2.1.2). It was suggested that 
backreactions from perturbations of wave-lengths larger than the Hubble radius, generated 
during inflation, can explain the accelerated expansion of the universe without the need 
for a cosmo logical constant. In [5] and [35] the authors computed the luminosity distance- 
redshift relation in a perturbed flat matter-dominated universe, taking into account the 
presence of cosmological inhomogeneities up to second order in perturbation theory in the 
adiabatic case. They found that the time evolution of the super-Hubble modes produce a 
large variance of the deceleration parameter, which could mimick dark energy. 

A number of studies (e.g. [291 EHl EH [79]) have criticised this result. In [351 ES] it was 
shown that the corrections only amount to a renormalization of local spatial curvature 
and hence that their magnitude is tightly constrained by observations. The authors of [58] 
pointed out that this argument only excludes the approach of [3 |15] , since in that work 
only the leading gradient terms in the energy-momentum tensor are considered. 

The effects of long-wavelength scalar metric perturbations with the complete effective 
energy-momentum tensor were studied in ^E\. Their starting point was the metric with- 



out anisotropic stress in Newtonian gauge (2.27). Then, both the Einstein and energy 



momentum tensor are expanded up to second order in metric and matter perturbations 
6(f. The linear equations are assumed to be satisfied, and the remnants are spatially av- 
eraged, providing the equation for a new background metric which takes into account the 
backreaction of linear fluctuations computed up to second order 

G^, = 87iG{T^, + T^,), (2.98) 



2.6 Macroscopic Gravity 
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where t^u consists of terms quadratic in metric and matter fluctuations and is called the 
effective energy- momentum tensor. Its explicit form was first calculated in [63] : 



roo = ^ {l2H (00) - 3 ((0)2) + 9a-' ((V^)^)) + ((50)^> 
+a-2 {(VS^f) + ^V'ilp) ( V> + 2V'm {<PS^) , 



(2.99) 



and 



Tij a 



'S,, (^i24H' + 16H) (02) +2AH (00) + ((0)^) +4 (00) - ^a~' (( V0)2)^ 



44^'(02) + ((5vi)2)-a-2((VM'>-#('500)-^V^"(^)(V> + 21^'(^) {<PSv) 

(2.100) 

The role of backreaction of super-Hubble modes in those inflationary models in which 
inflation ends through the reheating dynamics of if was considered in [5H]. This work 
found that then the effective energy-momentum tensor acts as a tracker during the period 
of radiation domination, but redshifts less rapidly than matter in the matter era. Using 
standard values for the preheating temperature and the amplitude of the inflation following 
preheating , they found that this mechanism leads to a possible explanation of dark energy. 
It should be commented that the issue of super-Hubble fluctuations is complicated by 
gauge-dependence. From the matter power spectrum it is quick to see that adiabatic 
fluctuations on such scales are highly-gauge dependent as in the Newtonian gauge the 
P{k) grows for decreasing k in contrast to the synchronous gauge where P{k) remains 
well-behaved. It has been shown [33l [3l] that super-Hubble isocurvature modes are gauge- 
independent but it is still fair to say that the backreaction of super-Hubble fluctuations on 
the evolution of the universe is a controversial topic. 



2.6 Macroscopic Gravity 

Macroscopic gravity is a non-perturbative, geometrical approach to resolve the averaging 
problem by its reformulation as the problem of the macroscopic description of classical 
gravitation. This theory was suggested by Zalaletdinov |110|, lllll 1112] . 

According to this approach, the average value of a tensor field (5^ "(a;) on the manifold Ai 
over a region S C at a supporting point x G S as 

(Qti^)) = ^ I iy-'TAx\x)...V/{x\x)...Q^;-{x')dn'. (2.101) 

Here Vs denotes the volume of the averaging region S and dVL denotes the invariant volume 
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element. Vu^{x', x) is the averaging bivector which is required to satisfy the two conditions, 

(i) limV/{x',x) = 6/, (2.102) 

(ii) Vp'^ {x' ,x) = 613'^ for vanishing curvature. (2.103) 

The first condition ensures the existence of the inverse bivector, V„^'(V^i)^^' = 6 J, and 
the correct hmit for Vj] ^ 0, the second condition guarantees a correct fiat space hmit. 

Since averaging and covariant derivation of a tensor field do not commute in a curved space, 
a second bivector, the coordination bivector W^"(x',x), is introduced. This bivector Lie- 
drags the averaging region from x' to x along integral lines, such that a comparison of 



the averages at different points is possible. It also obeys the two properties (2.102) and 



(2.103), but is also required to satisfy 

(iii) w/.y = 0, (2.104) 

(iv) W[/,p] + Wia^'y)^pf' = 0. (2.105) 

The first condition ensures volume preservation while Lie-dragging, and the second condi- 
tion ensures analyticity of the averaged tensor field. It is reasonable to choose the averaging 
bivector to coincide with the ccordination one, V/j'^ = W/?'^ , and such a choice turns out 
to simplify the formalism greatly. 

Next, the bilocal extension of the Christoffel symbols is defined as 

= + W/;.W/'). (2.106) 



When averaged by (2.101), this serves as afiine connection for the averaged manifold J^°/3-y. 



From this the curvature tensor M"^^^ is calculated in the usual way. It is called the in- 
duction tensor. Contracted it gives M^^,. 

Then the connection correlation tensor is introduced, 

Z^'p-y^ua = Z°'l3["i^vu\ = (-^"/3[7^^i/(7]) ~^ I3[y^^ ua\, (2.107) 

and the quantity 

Z"'^vi3 = 2^"^[p^|i,|/3], (2.108) 

where the index between bars || does not participate in the antisymmetrization, and the 
quantity Q^^^ = Q^^pu defined by 

Q"/375 = 2Z''/3[-y°'|^|5] (2.109) 

give the induction tensor Z"" p^^p — \5pQfj^y. Eventually, the averaged Einstein equations can 
be expressed as the field equations of macroscopic gravity: 



2. 7 Conclusions 
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where g^'^ is the averaged metric. 

If all correlation functions vanish and the macroscopic spacetime is highly symmetric, this 
equation reduces to the usual Einstein equations 

M., - \g^,r''M,, = -kt!^;'-\ (2.111) 

It has been shown |112] that, in the special case of the FLRW metric with an energy- 
momentum tensor of a perfect fluid, this macroscopic gravity equation takes the form of 
the Friedmann equations with additional terms pgrav and Pg^av that obey the equation of 
state Pgrav = — |pgrav, therefore mimicking dark energy. 

In additionally, with reasonable cosmological assumptions Buchert's formalism can be re- 
alized as a consistent limit of the macroscopic gravity formalism with identical corrections 
to the Friedmann equations [69] . 

2.7 Conclusions 

We have presented in this chapter an overview of standard cosmology, some of the problems 
it faces and discussed in some detail one particular approach at solving these, concerning 
the averaging problem. The different approaches to explain the observed acceleration of 
the universe as a consequence of deviations from exact FLRW symmetry without invoking 
the presence of a cosmological constant or dark energy have attracted much attention in 
recent years. However, there is also much criticism and controversy about this ansatz. As 
pointed out by Ishibashi and Wald in [10], an alternative model must be compatible with 
observations. They showed on a specific choice of time-slicing that a quantity representing 
the scale factor may accelerate even in Minkowski spacetime without any physical cause. 
Furthermore, they calculated Buchert's kinematical backreaction for two disconnected dust 
filled FLRW models, obtaining the requirement for an average acceleration although all ob- 
servers see only deceleration. They concluded that acceleration, as defined by the Buchert 
formalism, can easily arise as a gauge artifact produced by the choice of slicing or the choice 
of averaging domain. This shows quite clearly the need for a generally covariant approach 
to the averaging problem which is not subjected to any gauge ambiguities. Macroscopic 
gravity is a systematic attempt to set up such an approach, but it has not been widely 
accepted. 



CHAPTER 3 



.COSMOLOGICAL BACKREACTION FROM 

PERTURBATIONS 



3.1 Introduction 

There are two main ways in which averaging can be seen to be a vital, if imphcit, com- 
ponent in cosmology. On the one hand, cosmology is founded on the FLRW metric which 
is assumed to be the large-scale average of the true, inhomogeneous metric. This state- 
ment cannot be proven without a proper generally covariant averaging scheme capable 
of smoothing inhomogeneous metrics. While it is certainly logical that the metric of the 
universe on the largest scales should resemble an FLRW metric (and recent studies [SU [60] 
have considered this problem from an observational point of view) this is not a solid proof. 

The other way in which averaging in cosmology has a vital impact concerns the non- 
linearity of the Einstein tensor. Even should we possess a covariantly averaged metric 
on large scales we cannot directly employ this in the Einstein equation, since the Einstein 
tensor of the average metric is not equivalent to the average of the inhomogeneous Einstein 
tensor. The corrections to the dynamics are known as the "cosmological backreaction" . 
In the remainder of the thesis we will develop a generally covariant averaging procedure 
for metrics in general relativity, but in this chapter we consider a simple test-case to eval- 
uate quantitatively the deviations from a standard FLRW evolution arising from linear 
and mildly nonlinear perturbations. Such a calculation complements the recent studies 
by |100] , in which the authors reconstruct the impact of backreaction effects from the ob- 
servational data, |l2] where the authors evaluate the size of the effective density of the 
backreaction as a function of redshift in a structured Robertson- Walker model and [80] 
which considers a statistical "peak model" for large scale structure in an Einstein-de Sitter 
universe. While the current literature generally concerns exact inhomogeneous rather than 
perturbative models these models are not constructed to accurately model the universe. 
In contrast, the deviations introduced by linear perturbations are expected to be small 
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when averaged over a domain approaching the Hubble volume in size but the model is 
viable phenomenologically and the perturbations can be calculated numerically employing 
popular Boltzmann codes such as cmbeasy or cmbfast. Employing the Buchert formalism 
we will calculate the deviations in both the concordance, ACDM model and in a toy EdS 
universe with a low Hubble rate. We model mildly nonlinear ( "quasilinear" ) scales using 
a modified Halofit [91] code. 

This chapter is based on work performed in collaboration with Iain Brown and Georg 
Robbers [8] and was recently published in JCAP. In this chapter an overdot represents a 
derivative with respect to coordinate time t rather than proper time r, and in line with 
much of the literature we again employ a signature (— , +, +, +). 



3.2 The Buchert Formalism in Newtonian Gauge 

While much of the previous study into backreaction in perturbed FLRW spaces has been 
undertaken in synchronous gauge (see for example this gauge contains unphysical 

modes and the metric perturbation can grow to be relatively large. While this is not nec- 
essarily an issue, Newtonian gauge is unambiguous and the variable remains small and 
well-defined across almost all scales. Moreover, it is easily incorporated into the cmbeasy 
Boltzmann code and should serve as a complementary probe to the previous studies, al- 
though it should be remembered that we cannot directly compare calculations performed 
in different gauges. It has been appreciated for some time |1U31 177] that the impact on 
large-scale evolution from perturbations is not expected to be large, something clearly seen 
in Newtonian gauge; since the perturbations themselves are consistently small across al- 
most all scales the impact from backreaction would naively be expected to be at most of 
the order of 10~^. 



As we intend to evaluate the perturbations with a multi-fiuid Boltzmann code in Newto- 
nian gauge some modifications to the formalism presented in sections 2.3.1 and 2.3.2 is 
necessary. In particular, we do not identify the normal vector to a hypersurface with 
the 4- velocity of any one fluid and will instead connect the lapse function to the Newtonian 
potential. As and are not equal, we cannot employ the acceleration as in equation 



(2.69) and our Raychaudhuri equation will look slightly different. Moreover, the split of 
the stress-energy tensor with respect to the foliation does not coincide with the definitions 
used in the Boltzmann codes and this will introduce a further modification in both the 
Friedmann and Raychaudhuri equations. 

The Einstein equations separated with respect to this foliation can be written 

n + K^- K^Ki = IQnGp + 2A, Dj [K'^ - K) = SirGj', 
^kij = TZij - 2K^Kaj + KKij - SnGSij + ATcGhij {S - p) - Ahij - ^'^i^a 



3.2 The Buchert Formalism in Newtonian Gauge 
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where Sij = Tij and p = T^yu'^n^ . We also have the normal vector and extrinsic curvature, 
which we write as 

= [^,0^ , = (-«,0), K,, = -^h,j, (3.1) 

recalling that hij = dthij. We again define an average Hubble rate by 

3^ = ^ = - (aK) . (3.2) 
a-r, V 

The average of the Hamiltonian constraint can then be written as the average Friedmann 
equation, 

(^)^?|l(„V> + f(«V^S. + '^.) (3.3) 
with the kinematical backreaction 

Qv = {a'{K'-K;Ki))-l{aKf (3.4) 

and 

TZv = (a^n) (3.5) 

the correction arising from the spatial curvature. The trace of the evolution equation for 
the extrinsic curvature is 

ak - aK'-aKi + aV'Via = AnGa^p + S) - a^A. (3.6) 

Averaging this, employing the commutation relation and using the averaged Friedmann 
equation (3.3) gives us the average Raychaudhuri equation, 

''^ ^""^ ' \p + S)) + ^ {a') + l{Qv + Vr^) (3.7) 



aj) 3 3 3 

with 

Vjy = (aD'D.a) - (aK) , (3.8) 

the dynamical backreaction. 

Finally we have the integrability condition. Writing Hx> = (Ti.) where Ti. = —aK is a local 
Hubble rate, and using local matter continuity we find that 

6 + 2 = 12 (a \ ^^P+ IT ] - 

+167rG {{a^S) (H) - (a^SH)) + 6A {{a^U) - (a^) (H)) - AHvVv (3.9) 

It is easy to see that in dust-filled synchronous gauge models, with a = 1 and 5 = 0, 
the source on the right-hand side vanishes. For non-dust models, there is still an extra 
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source term dependant on the local pressure and the local and averaged Hubble rates. In 
[52| |53] the authors make much use of the integrability condition, employing it iteratively to 
recover the backreaction at higher-orders in perturbation theory. In Newtonian gauge the 
source term is somewhat complicated. Since we will be evaluating the backreactions from 
perturbations calculated in a Boltzmann code, we will not make use of the integrability 
condition. 

3.3 The Connection with Linear Perturbation Theory 

The Buchert equations are exact for any inhomogeneous model; for further progress we 
must specify this model. An increasingly common approach is to employ a model that 
averages out to be Robert son- Walker on large scales but is a relatively realistic approxima- 
tion to the local universe on smaller scales. We choose instead to use a linearly-perturbed 
FLRW model, employing Newtonian gauge. We consider a universe filled with a cosmo- 
logical constant and pressureless dust. This is a reasonable approximation to the current 
universe on the largest scales. Should the backreaction turn out to be significant the cosmo- 
logical constant can be reduced or set to zero, while should the backreaction be insufficient 
to account for observations it can still remain phenomenologically viable. Working to first- 
order in the gravitational potential (p, we consider only scalar perturbations and work with 
the line-element 

rfs^ = -(1 + 24>)dt^ + a2(t)(l - 24>)6ijdx'dx^ (3.10) 

where the scale factor a is distinct from the averaged scale factor ax>. We neglect the 
tensor perturbations as sub-dominant to the scalars although a full approach should nat- 
urally take these into account. When performing our averages across scales approaching 
the Hubble horizon, we consider the quadratic order in linear perturbations and neglect 
the averages of pure first- and second-order perturbations. Were we to consider averaging 
across relatively small scales then these neglected terms should also be taken into account 
for a full evaluation. Likewise we can neglect the second-order vector and tensor pertur- 
bations that scalars inevitably source, but again these should be taken into account for a 
full treatment on small scales. These issues would require a generally-covariant averaging 
procedure. Since we are considering purely scalar perturbations, all vector quantities can 
be written as the gradient of a scalar quantity; in particular, f * = —d'^ijj for some velocity 
potential ip. 

By unambiguously identifying the 3-1-1 and Newtonian gauge coordinates and expanding 
quantities up to the second-order in perturbations (see also [Mj), we can quickly see that 

a^ = l + 2(j), a^l + 0-^0^ a-^ ^ 1 - (f) + "^(p^ , ^ 1 - 20 + 40^. (3.11) 

We also have 

hij = a'^{t){l-2(j))5ij, /i'^' ^ a-2(t)(l + 20 + 402)(5'^'. (3.12) 



3.3 The Connection with Linear Perturbation Theory 
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We can now proceed to evaluate the geometric quantities we need in the averaged equations. 
The extrinsic curvature is 

aK] = - - 0(1 + 2</))^ 5]. (3.13) 

Using the Ricci scalar on the spatial hypersurface the curvature correction is 

Uv = {a'n) = 4 < (2 + 3( V0)' + 120V V) > (3-14) 
d 

where V = {d/dx,d/dy,d/dz) as in standard vector calculus. 
From the extrinsic curvature we can immediately see that 

2,— = - (aK) = 3 /- - 0(1 + 20) 
av \a 

implying 

^ = --(0(1 + 20)) (3.15) 

giving the relationship between the physical averaged scale factor and the scale factor 
employed in the perturbative approximation. If = then these coincide; this occurs in 
an Einstein-de Sitter universe, or when one considers a domain sufficiently small that its 
time variation can be neglected but sufficiently large that linear perturbation theory may 
be employed. We can now rapidly find the kinematical backreaction, 

Qv^6{{ ( ^-0(l + 20)y^-(^^-0(l + 20)y'] =6f(<^')-(<^y )• (3-16) 
Using the covariant derivative on the 3-surface and 

dX^ 30' -3-0(1 -20), (3.17) 

the dynamical backreaction is 

Vv = 3^ ^0(1 - 20)^ - 3 <J02^ + ^ (VV + 20VV - 2(V0)2> . (3.18) 
This gives us the geometric information that we need. 

Consider now the fluids. We separated the stress-energy tensor with respect to the normal 
vector 

n''= (1-0+ -02,0) . (3.19) 
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However, the stress-energy tensor employed in the Boltzmann code is defined with respect 
to the 4-velocity and for the case of pressureless dust is typically taken to be 

Tf^v = £u^u^ (3.20) 

where e is the energy density as defined in this reference frame, which will be linearised to 
e = p{l + S) for an average p. The density p and FLRW density e thus do not coincide 
unless the velocity is orthogonal to the normal, which is not in general the case. The 
4-velocity, to second-order in first-order perturbations, can be seen to be 

u^^(^l-(l>+^ {a\^ + 30^) , v'^ (3.21) 

with Vi — 5ijV^ . Then the density is 



(3.22) 



Scaled by the lapse function, the quantity that enters the Buchert equations is then 

{a^p) = p (1 + 5 + 20 + + 2(1)5) (3.23) 
and so we define a "density correction" 

r^, = ^p(5 + 2(/) + 205 + aV>. (3.24) 
The Friedmann equations can then be written as 

— =^Pm + i^ + ^F, -^ = - p^+ +AR (3.25) 
av J o 3 Ox) 3 3 

with corrections 

AF^Tjy-^ {Qt, + Uv) , AR^-^To + ^ (Q^ + Vv) ■ (3.26) 
The effective pressure and density of the corrections are then 

^Peff ^Tv-^{Qv + nv), 167rG'pefr = ^TZv - Qv - ^Vv (3-27) 
and so the effective equation of state is 

^'^ = -o\^ , ^ • (3-28) 



We take the domain to be large enough to allow us to neglect the averages of first order 
quantities on the background. Specifically, we take the domain size to be of the order of 



3-4 Ergodic Averaging 
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the comoving Hubble scale, across which first-order averages can consistently be neglected. 
While this implies that our analysis can only be taken to apply on the very largest scales, 
taking the domain to be so large also implies that we can invoke the ergodic principle to 
convert our spatial averages into ensemble averages; we exploit this in the next section. 



On such large scales, the modifications become 

TZv = -^((V0)' + 40V'0>, (3.29) 

Vv = (3.30) 

Qv = -6^(00)-3(02) + ^(</)V2</)-(V</))2>, (3.31) 

Tv = ^p(205 + aV>. (3.32) 



We can see that, depending on the signs of (0V^0) and (05), for an Einstein de-Sitter 
universe with = (as considered in [96]), one might get either an enhancement or 
reduction of both the effective Hubble rate and of the effective acceleration. 

3.4 Ergodic Averaging 

The simplest attack on the problem is to employ a Boltzmann code such as cmbeasy [21] 
or cmbfast [87] . The use of linear theory automatically implies that we are working on rel- 
atively large scales; we can render the system developed in the last section more tractable 
by taking our domain to approach the Hubble volume itself - large enough, at least, that 
we can employ the ergodic theorem and convert the spatial averages into averages across 
a statistical ensemble. (Similar approaches were also employed in |103l 177].) 

Consider first the general average (A()f:, x)i?(t, x)) where {A, 5} G {0, 5, f,0}. Then 

(A(t,x)5(t,x)) = jj^ (3.33) 

At linear order and in the absence of decoherent sources, different wavemodes are decoupled 
from one another and the evolution equations depend only on the magnitude of k and so 
we can write 

A{t,k) = a{k)A{t,k) (3.34) 
for any quantity. The initial correlation between two arbitrary quantities is 

(a(k)/5*(k')) = —P4k)5ik - k'). (3.35) 

Inserting this into the correlation and integrating once yields 

{A{t,^)B*{t,^))^ V4k)A{t,k)B*{t,k) — . (3.36) 
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The averages involving gradients are only a little more complicated; for ((V0)^) we have 



(W)') = / / (-zk)(zk')(<^(t,k)<^*(t,k' 



))e- 



-i(k-k')-x_ 



(27r)3 (27r)3 



dk 

k^v^{k)\ci>{tM^ 



'k ^ 

which is positive-definite as it should be. The other curvature term is 

dk 

k'V4k)\<j>{t,k)f^. 




(3.37) 
(3.38) 

(3.39) 
(3.40) 



This term is then negative definite. 



The kinematical, dynamical, curvature and density corrections (3.29 3.32), then become 

18 /" . r,_ . . .2 dk 



Qv 
Vv 



k'V^{k)\<p(t,k)[ 



k 



6 / V^{k) ^{t,k) 



2 dk 



-3 lV^{k)[ <P{t,k)^ +"^[cP{t,k)<P*{t,k) + <P*{t,k)<P{t,k) 



3a2 



+.^i0(u)r 



k ' 



Tr 



V 



^-^-pnevma\t)\v{t,k)ff 



V^k) {(Pit, k)5*it, k) + 0*(t, k)5{t, k)) 



dk 
~k 



(3.41) 
(3.42) 
(3.43) 
(3.44) 
(3.45) 
(3.46) 



For a more realistic approach where we separate the baryons from the cold dark matter 
one can readily see that the contribution from the density can be expanded as 

= '-f fk'v,{k)j:-pMt,k)f'^ 



V^k) [<j>{t,k)Y,PXit,k) +<l>*{t,k)J2pX{t,k)] -■ (3.47) 

a a / 

Since Boltzmann codes tend to be written in conformal time we must these into conformal 
time expressions using d/dt = {l/a)d/dri\ denoting conformal time quantities with a tilde 
we thus have 

(3.48) 



0' a la' 
a ' a a a 



in the expressions above. 
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3.5 Corrections to the FLRW Picture from Linear 
Perturbations 

We incorporate the above formalism into the cmbeasy Boltzmann code and run consis- 
tency checks with the cmbfast code, modified to output conformal Newtonian quantities. 
We will evaluate the terms across all post-recombination redshifts and employ an infra-red 
cut-off at the comoving Hubble scale, k^i^ = l/rj, avoiding the unphysical gauge-dependent 
super-horizon contributions. The small-scale limit A^max of our domain, determined from 
the stability of the integration with respect to changing fcmax is about fcmax ~ 30Mpc~^ 
and we integrate to /cmax = lOOMpc"^. Naturally, we do not claim that our results at such 
small scales are complete, merely that we are evaluating the contribution of such large-scale 
modes on these scales. 



At a linear level, there are scaling relations that hold to a high level of accuracy; it is 
obvious from Poisson's equation that 

oc S/P; (3.49) 

from Euler's equation v + {a/a)v oc V0 we can also predict that 

\v\ oc S/k. (3.50) 

From here, one may immediately state that on smaller scales where <^ and for 
an approximately scale-invariant primordial power spectrum, 

Qvik) oc ^, Tr^ik) oc P^ik) oc TZvik) oc ^ (3.51) 

where Qv{k) = 6P^(A;)|0p is the integrand of Qv, with similar definitions for Toik), 7lv{k) 
and Vvik). Qd is thus generally subdominant to the other corrections except on very large 
scales. This also implies that each correction is approximately of the form 

AD = aj ISf'^dk (3.52) 

for some constant a, reminiscent of the approximation in equation (32) of |103] . 



For our models we take a low-Hubble constant Einstein de-Sitter model with ilj, = 0.05, 
Qrn = 1 and h = 0.41 and a WMAPIII ACDM concordance model. In both cases we con- 



sider only adiabatic initial conditions. The left panel in Figure 3.1 shows the (integrands 
of the) four different correction terms at redshifts of z = 10 and z = for the EdS case and 
the right panel the same for the ACDM case, and in Figure 3.2 we present the correction 
terms as a function of z for the EdS and ACDM cases; the premultiplication by a^{z) acts 
as a volume normalisation. The subdominance of Qv{k) and proportionality between TZ-p, 
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T-p and Vv is very clear, with TZxi the strongest correction. 

From the proportionahty of the other corrections (which holds up to relatively large scales) 



we can write Tx, = tTZ-p and Vv = pT^-v and express the equation of state (3.28) as 



'^T^ ~ -\ i\ — % ) • (3-53) 



3 V 1 - 6t 

which is simple to evaluate numerically by selecting some pivot scale on which to evaluate 
the ratios p{z) and t{z). Doing so at the current epoch and a pivot scale of A; = O.OlMpc"^, 
we find p ^ 2/9 and t ~ 1/20. Both of these are also only slowly evolving. This gives an 



estimate of the average equation of state as wj} ~ —1/19. In the left panel of Figure 3.3 



we plot the evolution of wt> for both the model EdS and ACDM cases evaluated directly 



from equation (3.28). In both cases the effective equation of state from linear perturba- 
tions remains around the order wv ~ —1/19 and the correction terms to the usual FLRW 
as a whole thus act as a form of non-standard dark matter and not as a dark energy. 
These results compare reasonably well with the estimate in [103] where the author found 
wr, ^ -1/270 

To quantify the impact on the Friedmann and Raychaudhuri equations, define 

~ — 3a^ — ' ~ — 3a^ — 

as the standard contribution from baryons and CDM. We can then consider AF/Fm and 
AR/Rm as a well-defined measure of the impact of the correction terms in both EdS and 



ACDM models. By equation (3.27), AF/Fm for an EdS model is just pes/Pm- For the sake 
of clarity we choose to focus on deviations from the standard behaviour but, naturally, we 
can re-express these quantities in terms of the effective energy density and pressure. 



For the EdS case we can see in Figure |3.3| that the total impact from linear perturbations 
tends to ~ 4 X 10~^ on the Friedmann equation and ~ —3.2 x 10~^ on the Raychaudhuri 
equation; the corrections thus act with a positive effective density and with insufficiently 
negative pressure to accelerate the universe and instead act to decelerate it at a negligible 



level. The dashed curves in Figure 3.3 are for the ACDM case; they tend to impacts of 
~ 1.3 X 10~^ on the Friedmann equation and ~ —1 x 10~^ on the Raychaudhuri equation. 
The behaviour, as might be expected, is qualitatively similar to that in the EdS case and 
the impact is significantly less (by a factor of roughly 10/3 at the current epoch). 

This compares reasonably well with the approximation in |103j which in our notation is 

^ = ^(8vrGpJ f^dk. (3.55) 



^They also found wxi ~ —1/15 for a clustering cosmon field but we have not considered such a component 
here. 
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Figure 3.1: The correction terms Tj), Vv, Tlv and Qv as a function of A; at 2; = 10 (black) 
and z = (brown) for (left) the sample Einstein de-Sitter model and (right) 
the WMAPIII concordance model. 



Evaluating this approximation for the EdS universe, we find that this underestimates the 
correction terms with respect to our more detailed study by a factor of about 30%. In the 
same paper, the author found a "cosmic virial theorem" , 

Pm = -PeS (3.56) 

where is the correction to the matter pressure arising from gravitational interactions. 
Crudely modelling the matter pressure as ~ (l/3)p^ (f^) we find that underesti- 
mates —pes, again by ~30%, but the proportionality between them holds remarkably well 
for a wide range of redshifts. 



For the ACDM case we can consider an alternative normalisation that directly quantifies 
the impact on a universe with a cosmo logical constant, 

AF _ AF AR _ AR 

av/av Fm + A/3' \av/av\ |-Rm - A/3| 

The benefit of doing so is that it demonstrates the declining contribution of the corrections 
with respect to the cosmo logical constant; however, it also introduces a singularity at a 
redshift of z ~ 0.8 when the reference FLRW universe undergoes a transition from decel- 
eration to acceleration. The impact on both the Friedmann and Raychaudhuri equations 
tends to ~ 10^^ and the maximum contribution is at a redshift of 2; ~ 1.3. 
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Figure 3.2: The evolution of the corrections for (left) the sample EdS model and (right) 
the WMAPIII concordance model. 




Figure 3.3: Left: The effective equation of state wd as a function of redshift for a sample 
Einstein de-Sitter and WMAP concordance models. Right: The impact of 
the corrections onto the Friedmann equation (black) and the Raychaudhuri 
equation (brown), normalised in the ACDM case to the matter content (dashed) 
and to the standard equations (solid). 
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3.6 Corrections to the FLRW Picture from Quasilinear 
Perturbations 

In the previous section we evaluated the impact on the large-scale evolution of the universe 
from linear perturbations, demonstrating that it is small, as should be expected, but also 
that it acts as a dark matter and not a dark energy. In this section we employ the halo 
model to estimate the impact of perturbations on smaller, quasilinear scales which are 
much less understood. 



The publicly-available Halofit code [HI] converts a linear CDM power spectrum into the 
CDM power spectrum from the halo model. We employ a modified Halofit code that 
instead takes the linear matter power spectrum and estimates the nonlinear matter power 
spectrum. The square-root of this power spectrum can then be used to estimate the non- 
linear density contrast 5nl(^) and we can then recover the other relevant quantities from 
the scaling relations v oc 6/k, (p oc If we define 

m = ^ (3.57) 

we can employ this as a scaling factor to recover estimates for the quasilinear behaviour 
of S{k), v{k), (f){k) and (f){k). We can then employ the formalism we developed for the 
linear case to estimate the impact in the quasilinear case, retaining the same domain size. 
We should immediately note two problems with this technique. Firstly, while the scaling 
relation for the velocity is extremely precise for linear perturbations it is not for non-linear 
perturbations and so we have automatically introduced a source of error. Perhaps more 
importantly, we are employing a formalism developed for linear perturbations in which 
different wavemodes decouple from one another, allowing us to separate the statistical 
problem into transfer functions and a primordial power spectrum. This does not hold for 
a nonlinear problem. For both of these reasons, our results are only intended to be taken 
as good approximations on quasilinear scales until the velocity virialises, at which point a 
more detailed study is necessary. 



In Figure 3^ we present approximations for the corrections to the ACDM model at redshift 
z = 0. As before, the kinematical backreaction is strictly negligible and as we have merely 
scaled the previous corrections by the same quantity, the effective equation of state remains 
unchanged. The total impacts on the Friedmann and Raychaudhuri equations for both are 



shown in Figure |3.4[ We see that although the halo model provides a boost in power 
on smaller scales, and even though TZ-p and T-p in particular include factors of k"^ in the 
integral that increase the contribution from smaller scales, the impact is not significantly 
greater than from linear scales. More quantitively, for an EdS-universe the impact at the 
current epoch on the Friedmann equation is of order 5.6 x 10~^ and on the Raychaudhuri 
equation is of order —4.7 x 10~^. Normalised to the matter content, the impact at the 
current epoch from ACDM perturbations is of order 1.6 x 10~^ on the Friedmann equation 
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Figure 3.4: Left: The halo model approximations to the quasilinear corrections for the 
WMAPIII concordance model at z = 0. Right: The halo model approximations 
to the quasilinear corrections to the Friedmann (positive) and Raychaudhuri 
(negative) equations; the linear prediction is in brown. 



and of order —1.4 x 10 ^ on the Raychaudhuri equation. Normalised to the full evolution 
these become 4.4 x 10"^ and —2.7 x 10~^, respectively. 

Our results compare well with the observational estimate from Vanderveld et. al. |1U0] 
and the structured FLRW estimate of Khosravi et. al. [12] . It is also in broad agreement 
with the larger-scale estimates of Li and Schwarz [521 [53] although their method allows 
them to evaluate the backreaction on much smaller scales at which they recover signif- 
icantly greater impacts, up to about 10~^ at scales of 50Mpc. This should perhaps be 
treated with some wariness, however, since they work to fourth-order perturbation theory, 
at which level vector and tensor perturbations should be considered in a consistent ap- 
proach. Recently, Rasanen employed a statistical peak structure model at relatively small 
scales in an EdS universe [80] and recovered significantly larger modifications than our 
large-scale perturbative approach. Although the backreaction remained at the 10~^ level 
he found the curvature contribution to be significant. As with Li and Schwarz, Rasanen's 
results are not contradictory to ours given the different scales considered in the two cases. 
Moreover he considered only EdS universes and approximate transfer functions. 

3.7 Discussion and Conclusions 

We have presented quantitive estimates of the corrections in Newtonian gauge to a standard 
perturbative FLRW model from an explicit averaging procedure, which can be separated 
into four distinct quantities: the kinematic backreaction, which remains strictly negligible 
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across the scales of interest, the dynamic backreaction, the impact of spatial curvature 
and a correction to the FLRW energy density. In line with expectation, the impact from 
linear perturbations is insignificant and of the order of 10~^ for both ACDM and an Ein- 
stein de-Sitter model. Specifically, the effective energy density arising at the background 
level from inhomogeneities is peff ~ (4 x 10~^)p^ and pefr ~ (1-3 x 10~^)p^ for the EdS 
and ACDM models respectively. Moreover, we have presented estimates arising from halo 
model corrections on quasilinear scales which are not much larger than those from the 
linear perturbations, the largest contribution arising from the quasilinear modes in an 
EdS universe remaining below 10"'^; we find pcS ~ (5.6 x 10~^)p^ for the EdS model 
and peflf ~ (1-6 x 10~^)p^ for ACDM. This is in broad agreement with recent calculations 
[1001 H2l [53l [80] . Additionally, the total effective equation of state arising from the dif- 
ferent modifications does not act as a dark energy and instead as a dark matter with a 
slowly-varying equation of state w-d ~ —1/19 for z G (0, 10). The impact is to decelerate 
the universe and, contrary to other studies, the corrections impact positively on the Fried- 
mann equation implying a positive effective pressure. 

This does not, however, necessarily imply that there are no observational consequences 
arising from inhomogeneities, particularly on small scales, as our analysis is limited to 
the impacts in very large volumes and breaks down when mode-coupling and virialisation 
become significant. In particular, one can readily imagine situations in which the local 
Hubble rate is significantly larger than the global average (the so-called "Hubble bubble"). 
The impact of local inhomogeneities directly on luminosity distances has also long been 
well studied and remains an active area of interest. One avenue for further research would 
thus be to employ the same formalism and consider more accurate models on smaller scales. 
This remains very much an active field and one possibility we wish to consider is to employ 
the fully non-linear approaches to cosmology of Vernizzi and Langlois and Enqvist et. al. 

Ha ill EDI Eg. 

However, useful as this work will be, it will be necessarily limited until we have a usable, 
generally-covariant averaging procedure. On nonlinear scales one cannot neglect vector 
and tensor perturbations, which will be of a similar amplitude to the scalar perturbations. 
Without a generally-covariant averaging procedure, then, we cannot perform an average 
on relatively small scales without missing potentially significant contributions from the 
vector and tensor modes. Moreover, the Buchert averaging procedure relies on a vanish- 
ing vorticity, which holds true down to approximately galactic scales but breaks down if 
one wishes to consider, for example, the contribution from neutron stars. Perhaps more 
importantly, it relies on taking scalar projections of the Einstein equations and averag- 
ing these. This naturally loses much information about the system one is considering. It 
would be preferable to be able to directly average the Einstein and stress-energy tensors 
and so recover the large-scale behaviour directly. We will not address this point further 
in this thesis and instead refer to Zalaletdinov [112j for encouraging recent work in this area. 



To summarise, in this chapter we have demonstrated the existence of the cosmological 
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backreaction arising from Newtonian gauge perturbations to a FLRW universe on the 
largest scales. In the domains we have chosen the averages of the perturbations themselves 
can be taken to vanish and leading contribution comes from averages of the squares of 
the linear perturbations. We found that the total modifications to the FLRW picture 
from such modes remains of the order of lO^^p^, both for the concordance model and 
for a low Hubble rate Einstein-de Sitter model, and that their equation of state is only 
marginally negative and certainly not sufficient to drive an accelerated expansion. We have 
also extended the analysis to smaller scales using the Halofit code and determined that 
the modifications to the FLRW behaviour remain negligible on smaller scales. In doing 
so, we have demonstrated that one can quantitatively study the backreaction in realistic 
cosmological models. 



CHAPTER 4 



A GENERALLY COVARIANT AVERAGING PROCESS 



In this chapter we present a generally covariant averaging process in the framework of gen- 
eral relativity. It involves the decomposition of the metric into tetrads via a Lagrangian 
formulation and then parallel transport of the tetrads with a reformulated, general relativis- 
tic Wegner-Wilson line operator. In the first section we give an overview over the various 
constituents that we need to formulate the process in the subsequent section. Sections of 
this chapter are based on work undertaken for my diploma thesis [7j. 



4.1 Constituents of the Averaging Process 
4.1.1 Wegner-Wilson Lines in General Relativity 

In QCD, the Wegner-Wilson line operator is used to transport quantities from one point 
to another at finite separation in a gauge-invariant manner. To build a covariant averag- 
ing process, we need the general relativistic analogue of this operator, one which parallel 
transports tensor quantities from one spacetime point to another in a generally covariant 
way. To define this operator, which we will call the "connector" , we follow the calculations 
in [HS] and reformulate them in the framework of general relativity. 



Let r^(x) denote four 4x4 matrices with the components 



(4.1) 



where T^j^{x) are the Christoffel symbols (A. 10). We define the connector along a spacetime 



curve CxQxi, which connects the points Xq and Xi (Figure 4.1 a)), in the following way: 



V{xi,xo;Cxoxi) = Pexp I - / 



(4.2) 
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Figure 4.1: Different curves in spacetime: a) connecting two points, b) connecting three 
points and c) oriented in opposite directions. 

V denotes the path-ordering of the integral. To show what this means, we choose an 
exphcit parametrisation for the curve, 

= {{zni^" = ^^r). To < r < n, z>^{to) = z'^{n) = a;^}, (4.3) 

such that we can express the connector as 

\/(rO =Pexp £d/-^ TMr))) ■ (4.4) 

For Ti = Tq the connector reduces to 

Vin = To) = 1. (4.5) 

Now we divide the interval (ri — tq) into small pieces, denoted by Api = pi — pj„i for 
z = 1, . . . , A^. For very fine granulation the exponential in the connector can be expanded 
in a power series, 

Vin) = jhn ((l - Ap^^^^TMPn))) . . . (l - Ap,^^TMPi)))) ■ (4-6) 
Path-ordering means that matrices with higher values of p are placed on the left, such that 

To = po < pi < . . . < pn = n. 
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Due to the path-ordering, the connector has the property that for two adjoining curves 
Cxoxi,CxiX2 (Figure 41] b)) the connectors are multiphed 



V{X2, Xq; Cxoxi + C^ixa) = V{X2, Xi, Cx^x^ W^Xi, Xq] Cxox^)- (4.7) 

Furthermore, let Cx^^xo be the same curve as Cx^xi, oriented in the opposite direction (Figure 
4.1| c)). Then the connector obeys 

V{xi,xo;CxoxjV{xo,xi;Cx^xo) = 1- (4-8) 

This can be seen by sphtting and expanding the curve Cx^xo i^i the same way as Cxgxi- 
Since 2;''(r) has then the reverse sign, we find 

V{ro) = jim ((l + Ap,^^TMPi))) ■ ■ ■ + Ap^^^^T Mpn))) ) , (4.9) 



such that in the product (4.8) all Api terms cancel 



4.1.2 Transformation Law for the Connector 



The connector (4.4) is the solution to the differential equation 



dr dr ' 

with r = Ti and the initial condition, from ( 4.5[ ), 



r'MT))V\iT). (4.10) 



V\{to) = 5^. (4.11) 

Under a coordinate transformation z'{t) = f{z{T)), contravariant and covariant vectors 
transform as 



A'^ = —A^ ^ U^^A^ and A'^ = —A^ ^ {U'^ ,A^. (4.12) 



Therefore, the transformation laws for dz'' /dr and d/dz'' are 

and the Christoffel symbols transform as 

rU^'ir)) = U^,iU~y,{U^')\T';Mr)) + U',^^^^ (4.14) 
Furthermore, we find from d{U^ ^{U^^)^ s) / dr = that 

dr^ ' ' dz'f dr dz'^ ^ ■ ^ 
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Consider now the following quantity: 

V^air) = U%{T)V^,{r){U-\ro))\. (4.16) 
Differentiation of this with respect to r gives 

= if -^iU-'Y^V\ - U^,-^T^^,{U-')\U'^V\{U-\r,))\ (4.17) 



dr ''^ 



This equation is identical to that for the connector (4.10) in new coordinates, and be- 
cause of the existence and uniqueness theorem, which states that there exists only one 
solution to a first order differential equation which satisfies a given initial condition, we 
can identify V^ai^) with V'^aiT)- Thus, ( |4.16 ) is the transformation law for the connector. 

This transformation law shows that the first index of the connector transforms as a con- 
travariant vector at z{t), while the second index transforms as a covariant vector at z{to). 
In this sense the connector connects different spacetime points. Hence, the indices of 
V^ai^) are lowered and raised by applying the metric on the corresponding point, 

VJ{t) = g^M^W'u{T)g''^{z{To)). (4.18) 

Differentiation with respect to r gives 

dVa^ir) d 



{go.Mr))V^Ar)9'''{z{ro))) 



dr dr 

= '^-^^'^^^V\{r)g^^{z{^^^^ 

dr dr 

= ^^"^j^(^)^ ^M^(^(r))^,,(^(r))\/^(r)^-^(^(ro)) (4.19) 

-9c.Mr))^^T>i^{z{T))V'>^{T)g'^^{z{T,)) 
_ dz^ir) [ dgai,iz{T)) 



dr \ dz^ 



Inserting the explicit expression for the Christoffel symbols (A. 10) yields the differential 
equation 

ar dr 
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and thus (4.18) is the solution to this differential equation with the initial condition 

Vc^^iro) = 6J. (4.21) 
We refer to this as the "hatted connector" . 



4.1.3 Effect on Tensor Quantities 



We define the application of the connector (4.4) to a contravariant vector A"{z{to)) at the 
point z{to) as 

A^izir)) ^ y^(r)A"(z(ro)). (4.22) 
Differentiation with respect to r gives 



dz^'ir] 



r'Mr))V\{T)A-{z{To)) 



dr '^"^ 



(4.23) 



where we used (4.10). Therefore, the covariant derivative of A^{z{t)) along the curve z{t) 
vanishes, 

^^!i£M) . l^Mr)) + "^KMrMr)) = 0. (4.24) 



With (4.5) we find that A^{z{tq)) = A^{z{to)) and conclude that the application of the 



connector to a contravariant vector A"{z{to)) at z{to) yields the unique parallel trans- 
ported vector yl^(z(r)) at z{t). 



For a covariant vector A.y{z{T)) at z{t) we find 

A,{z{r)) = g,Mr))A^i4r)) = g,Mr))V^a{r)g''\z{To))As{z{To)) 
= %\T)As{z{ro)), 



(4.25) 



where we used (4.18). Differentiation with respect to r and employing (4.20) gives 



^i,(r) = ^ (v,\r)As{z{ro)) 



'^''^^^-rJz{T))Vj{r)As{z{r,)) 



dr 



(4.26) 



Hence the covariant derivative of A^{z{to)) along the curve z{t) vanishes, 

DA^(t) d ~ . . dz'^ir)^,, 

— ^ = j^AM-^)) - -^T^,M^))A,^z{t)) = 0, 



(4.27) 
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and with (4.21) we find Aj{tq) = A^{to) and thus the apphcation of the hatted connector 



(4.18) to a covariant vector As{z{to)) at z{to) yields the unique parallel transported vector 



A^{t) at z{t). 

In a similar way, one can show that general tensors are parallel transported by applying 
the connector to each contravariant and the hatted connector to each covariant index. We 
will now turn to two special cases, namely, the tangent vector and the metric. 

A geodesic is defined as a curve z{t) whose tangent vector is always self-parallel. Therefore, 
such a curve has to obey 

^(r) - VUr/^M. (4.28) 
Differentiation with respect to r gives 

Thus, we find the following differential equation for the curve, identical to the geodesic 



equation (A.9): 



^M-.e(.M)|^MfM^o. (4.30) 

The parallel transported metric q^u^'t) obeys the differential equation 

^9,Mr)) = ^ {%-{T)V/{T)g^,{z{ro)) 

dz^ir) 



dr 



^iMr))V-{T)V/{T)g^Mro)) (4-31) 

+ V,-{r)^^rl{z{r))V/{r)g^Mro)) 

dz^(r) d 
= {^iM^)y9pMT))+TiMr))g,M'^))) = ^9,u{zir)) 

with the initial condition 

^M^(i"o) = 9tM^iTo). (4.32) 

From the existence and uniqueness theorem we conclude that ^^j,(z(r)) = g^^{z{r)). This is 



consistent with the fact that the covariant derivative of the metric always vanishes (A. 12). 
As a consequence, we find an important relation between Xq; C^qx^)) and Xq; C^^x^] 
From 

5\ = g^^{z{r))g,,{z{T)) = g''^{r)%'^{r)V/{r)g^p{T,) = V\{t)V.^{t) (4.33) 



and (4.8), we conclude that in general 

V{xi,xq]Cxox^) = {V'^){xo,xi;Cx,xo)- (4-34) 
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4.1.4 The Tetrad Formalism 

As we saw in the last subsection, the parallel transported metric g^^{z{T)) is identical to 
9fiu{z{r)). This is the reason why a covariant averaging process such as the one suggested 



by Isaacson (2.49) cannot be used to construct an averaged metric. To overcome this prob- 



lem, we decompose the metric into a set of vectors with the help of the tetrad formalism. 

Consider a set of four continuous, pointwise orthonormal vector fields Ea^{x) with a = 
0, 1, 2, 3. An orthonormal set must be linearly independent, and so Ea^{x) forms a vector 
space basis for the tangent Minkowski space at each point x of the manifold. Three of the 
vectors of this basis are necessarily spacelike and one is timelike. We shall always so label 
the vectors that Eq^{x) is timelike. Then the condition of orthonormality may be written 
as 

Ea^{x)Eb''ix)g^^ix) = riab, (4.35) 

where r]ab = rj""^ = diag(l, —1, —1, —1) is the Minkowski metric. We can now interpret the 
index a as a Lorentz vector index which is raised and lowered by the Minkowski metric. 
Then we find 

E\{x)Eb''{x) = 51 (4.36) 

Provided the manifold is time and spacetime orientable, we choose Eq^{x) to be future- 
pointing and Ei^[x) to be a right-handed triad of vectors. In this case Ea^{x) is called a 
right-handed orthochronous Minkowski tetrad. 

The components of a given tensor field T^'^ p„{x) in the orthonormal basis are 

T"^,(x) = E\{x)E\{x)Ef{x)E/{x)T^%,{x). (4.37) 
If the Lorentz vector fields are subject to a spacetime-dependent Lorentz transformation, 

E\{x)=k\{x)E\{x), (4.38) 

the metric transforms as 

g^,u{x) = E^f,{x)E'^u{x)vcd{x) 

= A\{x)E%{x)A\E\{x)vUx) 
= E%{x)E\{x)vab{x) 
= 9f,v{x). 



The metric is therefore gauge-invariant under local Lorentz transformations and only a 
certain class of local tetrads is determined by the metric. 

Finally, we remark that under general coordinate transformations the tetrads still transform 
as coordinate vectors, 

E\{x) = g^E%ix). (4.40) 
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Figure 4.2: Parallel transport of all tetrads in a spacelike geodesic R region to the center. 

4.2 The Averaging Process 

To perform any sort of averaging we first need to define the domain of averaging. Let xq 
be a particular point of the manifold. Provided the domain is a sufficiently small neigh- 
borhood around Xq, it contains no conjugate points and hence there is a unique geodesic 
between each point of the domain and the center xq. We take the geodesic congruence 
generated by all spacelike vectors of the tangent space at xq. We then parametrize the 
geodesies by their length s and cut them off at the geodesic distance s = r = const, from 
the origin. We call the resulting region the spacelike "geodesic region" R of radius r around 
Xq. This region is bounded by the congruence of null geodesies. 

The metric inside the geodesic region is then decomposed into a right-handed orthochronous 
Minkowski tetrad, 

g^,{x) = E\{x)E\{x)r,ab. (4.41) 

Since the result of the averaging process will depend on the choice of E"-^{x), we need to 
choose a particular gauge. This is done by emplying a Lagrangian formulation. We define 
the "maximally smooth tetrad field" (ii^Ms)"/t(^) be the tetrad field that minimizes 



JR 



(4.42) 



4-3 Summary 
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within the geodesic region R with the Lagrangian 

Cms = {D,E%{x)) {D,E\{x)) g^^^{x)gP''{x)r^ab. (4.43) 
This Lagrangian describes the variation of the covariant derivative of the tetrads. It should 



be noted that (4.42) has to be minimized within a certain class of tetrads E^-^ix) that obey 
(4.41), namely, the group of proper, orthochronous Lorentz transformations. Furthermore, 
since we want to explicitly break gauge-invariance, the covariant derivative does not 
include the spin connection. 



Recall the expression for the hatted connector Vy^{xQ, x'\ C^'x^^-, defined by (4.34) and (4.2), 

i>(xo,x';C,,,J = Pexp {-I dz^ T^{z)] . (4.44) 
The hatted connector allows for the parallel transport of the maximally smooth tetrads 



(EMsTtiix) from all points x' in the geodesic region R to the center xq (see Figure 4.2), 
where we perform an average defined by 

{{EmsTu{xo)) = [ f{xo, x'- C,,,ji>/(xo, x'- C,'xo){Ems)%{x')V-9{x') d^x' . (4.45) 



R 



Here, /(xq, x'; Cx'xq) is a positive weighting function along the geodesies Cx'xq, which smoothly 
falls off to zero at the geodesic distance r and is furthermore normalised to unity over the 
geodesic region R, 

f{xo, x'; Cx'xo)V-9{x') d^x' = 1- (4.46) 



R 



This procedure has to be repeated for all points x of the manifold. We can then recompose 
the averaged metric from the averaged tetrads by 

{g,,{x)) = {{Eus)%{x)) {{Ems)K{x)) r^^b- (4.47) 



4.3 Summary 

In this chapter we have defined a covariant averaging process which can be used to smooth 
metrics within the framework of general relativity. We presented a Lagrangian approach 
to determining a maximally smooth tetrad field used to decompose the metric at each 
point in the averaging domain. We defined a general relativistic analogue of the Wegner- 
Wilson operator to parallel transport all tetrads in the averaging domain to the centre 
and we defined an average at this point. The averaged tetrads can then be employed to 
rebuild the averaged metric at the centre point. Repeating the procedure for all points of 
the manifold finally yields the averaged metric. To demonstrate the functionality of this 
process and to visuahze its effects, we apply it to various toy models in the rest of this 
thesis. 



CHAPTER 5 



AVERAGING CONSTANT CURVATURE SPACES IN TWO 

DIMENSIONS 



As we have seen in section 2J^, standard cosmology is based on the assumption of a homo- 
geneous and isotropic universe. The only spacetime metric consistent with this assumption 
is the FLRW metric (2.2), which is homogeneous and isotropic on constant time hyper- 
surfaces. This provides us with a strong physical motivation for studying such geometries. 
A more immediate motivation concerns the averaging process itself. Surfaces of constant 
curvature should be invariant under any reasonable averaging process and testing this 
should be an immediate aim. To study the geometry of constant curvature spaces, it is 
easiest to consider the case with two spatial dimensions. In this case, the space of constant 
curvature is either a positively curved two-sphere {S'^), a flat plane (H^), or a negatively 
curved hyperbolic plane [H'^)- In the next section we will apply the averaging process to 
the two-sphere (S*^) and investigate its effects in detail. The analysis of the flat plane (H^) 
and the hyperbolic plane {H'^) will be the subject of the subsequent sections. 



5.1 Averaging the Sphere 

The two-sphere of radius a can be embedded in a three-dimensional Euclidean space with 
coordinates (^^, ^^). This enables us to visualize it and to easily read off the line element. 
In spherical coordinates 

= a sin 9 cos (p, 

= a sin 9 sin 0, (5. 1) 

= a cos 9, 

where < 9 < n and < (f> < 2n, the line element is given by ds^ = a'^d9^ + a? sin^ 6*^0^. 
There is an obvious similarity between this line element and the FLRW line element (2.2). 
We could equally have used coordinates (r, 9) in which the line element takes the form 
ds^ = a^(l — r^)^^(ir^ + a^r'^d9'^. However, both line elements lead to metrics whose 
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Figure 5.1: Stereographic projection of the two-sphere into the (a;^a;^)-plane 

corresponding geodesies have a comphcated parameter representation. To circumvent this 
problem, we take the hne element ds'^ = a^dO^ + sin^ Odcl)^ and furthermore apply a 
stereographic projection into the tangent plane of the south pole. The geodesies through 
the origin of the projection plane are then given by straight lines. This simplifies the 
calculation of the connector and hence the calculation of the averaged metric considerably. 

5.1.1 The Metric from Stereographic Projection 

The stereographic projection of the two-sphere into the tangent plane of the south pole is 



illustrated in Figure |5.1[ The corresponding coordinate transformation is given by 

(5.2) 



a = (/), 

r = 2atan(7r/2 - 6/2). 
In terms of the new coordinates, the line element becomes 

16a^ 



ds' 



, ^ (ir^ -|- sin^ (vr — 2 arctan ( — )) dcx^ . (5.3) 

Using the addition theorems this can be simplified to the following form: 
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This reduction of the hne element to a simple form comes at a cost: The angle a is 
undetermined at r = 0, causing severe problems. Therefore, we change to Cartesian 
coordinates on the projection plane, 



= r cos a. 
x"^ = r sina. 



(5.5) 



If we additionally introduce the function 

L\x) = 4a' + {xy + {x''y, (5.6) 
we can express the line element in the following compact form: 

ds' = ^ {{dx^f + {dx^f) . (5.7) 

The associated metric and its inverse are 

where the indices i, j G {1,2}. 



5.1.2 The Geodesies and the Connector 



The two-dimensional analogue of the four-dimensional tetrad is the "dyad". To find the 
connector in the (x^x2)-coordinates, which parallel transports dyads from any point of the 
geodesic region to the origin, we first need to determine the geodesic congruence generated 
by all tangent vectors of the tangent space at the origin. The partial derivatives of the 



metric (5.8) 

yield the Christoffel symbols 

pi 
^ jk 



dgij 
dx^ 



Axk 



L{x) 



i9ij 



X 



L(x) 



{xjSi + XkS] - x'6jk) . 



(5.9) 



(5.10) 



Inserting these into the geodesic equation (A. 9) gives the differential equation for the 
components z^t) of the geodesies. 



d'^z' 
~d^ 



dz^ dz^ 



dz^ dz^ 



^ dr dr ~'~ ^'''^ dr dr 



0. 



(5.11) 



For the initial conditions z*(0) = and dz'^/dT{0) = const, the solution to this equation is 

r \ dz"" 



Z^T 



2a tan 



2a/ dr 



(0). 



(5.12) 
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With 



L(.W)-4a^co=-(I;)((|; 



(0)) H-^^(O) 



we find for tlie arc lengtli 



Tlius, if we clioose r to be tlie arc lengtli, r = s, we find 



-1/2 



dr 



(5.13) 



(5.14) 



(5.15) 



This property makes it possible to introduce an angle 7 with < 7 < 27r, such that the 
initial conditions for the geodesies are given by 



dz^ 
dr 



(0) = cos 7 and 



dz^ 
dr 



(0) = sin 7. 



(5.16) 



Recall the general form of the connector (4.2) 



y(xi,xo;Ca;oxi) = Pexp 



dz' Ti(z) 



(5.17) 



We parametrize the connector by r and insert the Christoffel symbols (5.10) and the 



geodesic (5.12) into this expression. Considering the parallel transport from the origin 



parametrized by r = to the point parametrized by r gives 
V{T,0;Cor) =Vexp 



Pexp (^-21,1 (cos 



Applying the inverse of this, 



V{0, T- C,o) = V-\t, 0; Co.) = cos^ j 1, 



(5.18) 
(5.19) 
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to a contravariant vector at z{t) yields the parallel transported vector at ^(0) = 0. The 
parallel transport of contravariant vectors is achieved through the application of the hatted 
connector (4.18), 

V^/(0,r;C.o) =^?,A.(0)\/^(0,r;C.o)/(r) = (y^)/(r, 0; Co.). (5.20) 

The hatted connector appears in the averaging process and according to the calculation 
above its explicit form in this case is given by 

\//(0,r;C.o) = cos-2(^)5i. (5.21) 

5.1.3 The Maximally Smooth Dyad Field 

In this section we want to find the maximally smooth dyad field. We start by writing the 
metric in terms of dyads, 

g,,{x) = E\{x)E\{x) Sab. (5.22) 

In two spatial dimensions the indices a,b, . . . are raised and lowered by 6ab and the indices 
... by gij. Let a reference field be 

The maximally smooth dyad field has to be a rotation of this dyad field by a position 
dependent angle 0(x^,x^): 

rpa _ Tja pfe ^,-^.1^ rra _ ( COS0(x\x^) sin0(x\x^)\ (^0A\ 

E,-UbE, with U b-[_^.^^., 2^ ^os6(x\x') ■ (^-2^) 



The partial derivative of the reference field is 
Hence, the covariant derivative can be expressed as 



(5.25) 



Ai?", = d.E'^j - T%E\ = ^ {x,5^ - x'g,,) E\, (5.26) 

The partial derivative of the rotation matrix Uab is 

d,U\ = e\U\{d4), (5.27) 
where e""}, is the totally antisymmetric tensor density, 

e%=(_°i J), (5.28) 
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and so we find for tlie covariant derivative of the maximally smooth dyad field: 



DiE'^j = {d,U\)E'^ + U^D.E'j = e\U%{d4)E'', + —U\ (a;,(5f - x'^Qi,) E\. (5.29) 



We will furthermore need the following identities: 



(5.30) 

(5.31) 
(5.32) 



We now find the maximally smooth dyad field by varying the Lagrangian (4.43). For the 
chosen reference field, this Lagrangian is 

£ms = (A^^) {DkE'i) g'^g^'dab 



Jab 



= e\U%E',e'fU^d{di^){dk(P)E''ig'^g^'5ab + U\{D,E',)e'' jU^ d{dk(p)E''ig'^g^'5ab 
+ U\{D,,E^)e\U\{d4W,g'^g^'5ab + U\U\{D,E'j){Dj,E''i)g'^g^'5ab. 

(5.33) 

This can be simplified by applying the identities (5.30 5.32) and inserting the expression 
(5.26) for the covariant derivative: 



£ms = 6,aE^,E\g^'g'\d4){du(l>) + 2e,, — (x.^f - x^g,,)E'^E\g^'g'\dk^ 



(5.34) 



With (5.22) and (5.23) and this expression gives 



Introducing the vector 



the Lagrangian becomes 



u'' = e,dE'^E\—{g 



mk^l _ rf.'rriglk'j 



(5.35) 



(5.36) 



Cms = 2g'\di<P){dk<P) + u\dk^) + —^gmnx'^x^. 
If we now vary the corresponding generally covariant action 



(5.37) 



S 



'R 



R 



d'x^ Cms = I d'x^ { 2(9,0) (9,0) (7*'= + n''-(9,0) + -^^7™„x"x" ) (5.38 
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with respect to 0, we find 



5S = / <fx^ {Ag'\di^){dM) + 2u\dk 



R 



u 



(5.39) 



Jr 



The fact that the integration area R is finite leads to two conditions that ought to be 



fulfilled by the dyad field. The first term in (5.39) can be converted to a line integral 
according to Gauss' law, 



dR 



ds ^nk{Ag'\d,(l)) + 2u'^) 50 = 0. 



For arbitrary 50 this leads to the boundary condition 



on OR. With the general identities 



90 1 , 



1 



(5.40) 



(5.41) 



g'^dkgij = -dug 
9 
1 



(5.42) 



the covariant derivative Dk of any vector can be expressed as 

DkA' = dkA'^ + lg'"'idigkm)A' = dkA' + ^idig)A' = {^y'dki^ A'^). (5.43) 



Hence, for arbitrary 50 the second term in (5.39) leads to the differential equation on R 

(5.44) 



9'' {dA^ 



-Dku\ 



In this case the covariant derivative of u vanishes, DkU = 0, as can be seen from (5.36). 



The geodesies through the origin (5.12) are parametrized by r and 7, where 7 is the initial 



direction of the geodesic, as can be seen from (5.16). We will call these coordinates the 



"geodesic coordinates" . They are parametrized as follows: 



z^{T,'-f) = 2a tan ( — ) cos 



7 



z (r, 7) = 2a tan y— | sm 7 



(5.45) 
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where < 7 < 27r and < r < 00. The contour hne dK of the geodesic region R is defined 
as all points of geodesic distance s = r from the origin. It therefore has the parameter 
representation 

a^{^) = 2a tan ( — ) cos 7, 
\2a/ 

0^(7) = 2a tan (^^j sin7. 
The coordinate invariant definition of the normal vector to dR is 

da'' 



^kl 



(97 



and the normal unit vector is 

nk = {9'^{r,l)viVjY^^'^ Vk- 
Explicit calculation shows that UkU^ = and we find the differential equation 



in the area R 



and 



dn 



on dR. 



(5.46) 



(5.47) 



(5.48) 



(5.49) 



The solution to this equation is = 0. We have thus seen that the dyad field (5.23) was 
already the maximally smooth dyad field, and so 



2a 
L(x) 



(5.50) 



5.1.4 The Result of the Averaging Process 

In the last sections we have determined all the constituents that we need for the averaging 



process (4.45). The general form of the average is 



{{Eusrj{x))= / f{x,x';C,>^)^\x,x';C,,,){Eusri{x')^/9{^ d^x'. 



(5.51) 



Since we calculated the connector (5.21) in terms of the affine parameter r it is use- 



ful to carry out the integral along the geodesies through the origin (5.45), which were 



parametrized by 7 and r. Expressing the maximally smooth dyad field (5.50) and the 



connector (5.21) in terms of these coordinates, we see that 

(EMs)^(7,r) =cos2 



-) 




2a) 




) 


{- 


\2a 





(5.52) 



and insert these into the averaging process (5.51). In terms of the coordinates (5.45) the 
square root of the determinant of the metric is 
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JQ 



r r2^ X (5-54) 

/(^(7, r)) v/^(^ drd^] 5% = 5%. 



Expressed completely in terms of 7 and r, the averaged dyad field at the origin becomes: 
{{EusTAm)= f r'^/(^(7,r))K/(^(7,r))(£;Ms)^(7,r)V'^(^t/rrf7 

') 

~'o / 

The term in brackets equals one since the averaging function f{z{'~f,T)) is taken to be 
normalized to unity over all space. Hence, the averaging process returns the averaged 
metric at the origin 

(^7,,(0,0)) = ((^Ms)^(0,0)) ((^Ms)^(0,0)>5,, = S,,. (5.55) 

To get the full averaged metric, we need to repeat this procedure for each point of the 
sphere. Doing so, however, leaves the results for each point expressed in different coordi- 
nates. To express them all in the same coordinate system we need to know the transfor- 
mation law which we now turn to 

5.1.5 The Coordinate Transformation 

Let the coordinates of the sphere be x = (x, y, z) and those of the stereographic projection 
plane at the south pole X5 = (0,0,— a) be (x^x^). We want the coordinate transfor- 
mation to the coordinates ix^^x^) of the stereographic projection plane at another point 



= [xq, Ho, zq) = (a sin 6q cos 0o, a sin Oq sin o- cos 6^0) of the sphere (see Figure 5.2 ). This 
coordinate transformation can be decomposed into three steps: the inverse stereographic 
projection of (x^x^) onto the sphere x = {x,y,z), the rotation of xq to the south pole 
X5 of the sphere and the associated transformation to the coordinates x = {x,y,T), and 
finally the stereographic projection into the tangent plane of the new south pole (x^,x^). 

The inverse stereographic projection from (x^, x^) to the coordinates of sphere x = (x, y, z) 



is a composition of the coordinate transformations (5.1), (5.2), and (5.5): 



4a^x^ 4a^x^ / 8a^ 



The rotation D of the sphere consists of three rotation matrices D^, Dy and D~^. First Dz 
rotates the sphere around the z-axis until the point xq intersects with the (x2;)-plane then 
Dy rotates the sphere around the |/-axis until the point D~^xo is aligned with the south 
pole, and finally D^^ rotates the sphere back to the original orientation. The rotation 
D = DzDyD~^ is of the following form: 

sin 00 0\ 

cos 00 . (5.57) 
1/ 
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The stereographic projection of the rotated sphere with the coordinates x = (x, F) into 
the tangent plane of the new south pole with the coordinates is given by 

~i 2a5; ~ 2ay 

= -, a;2 = ^. 5.58 

[a — z) (a — z) 

The explicit transformation from the coordinates to the coordinates (x^,5;^) and 

vice versa are given in appendix |Bj 

5.1.6 The Mobius Transformation 

If we identify the [x^x^)- and (x^x^)-projection planes with complex planes C, we can 



express the coordinate transformation (B.4) as a Mobius transformation. The general 
form of a Mobius transformation is 

~ "" + ^ (5.59) 



cz + d 



where a, 6, c, c? G C and ad — be ^ 0. If we use extended complex planes C U {oo}, we can 
furthermore include the north pole of the sphere in the transformation. In that case the 
point z = —d/c is mapped to the complex infinity z = oo and the point z = oo is mapped 
to = a/c. 

One method of finding the Mobius transformation that describes the coordinate transfor- 



mation (B.4) is to insert the complex numbers a = ai + ia2, b = bi + ib2, c = Ci + ic2, and 



d = di + id2 and to compute 

~_~ _ (fli + ia2){x'^ + ix'^) + {bi + ^62) 

^ ^ {ci + ic2){x^ + ix"^) + {di + id2) 

{aix^ — a2x'^ + bi) + i{a2X^ +aix^ + b2) {cix^ — C2x'^ + di) — 2(02x^ + 01x^ + ^2) 



(ciX^ — C2X^ + (il) + i(c2X^ + 01X^+6/2) (ciX^— C2X^ + (ii) — i{c2X^ + cix'^ + d2) 

explicitly. Then a, b, c and d follow by comparison of coefficients. 

A more elegant approach originates from Mobius transformations preserving the cross-ratio 

Z — Zi Z2 — Z^ Z — Z\ Z2 — Z'i 



Z2 — Zi Z — Zz Z2 — 23 



(5.61) 



From here one can compute the Mobius transformation if at least three points and their 
images zi — > Si, Z2 — > 22, and — >• S3 are provided. 

To construct the Mobius transformation in this way, we list here some particular points 
of the (x^x^)-plane and their associated images in the (x^x^)-plane. First of all, there is 
the origin of the (x^x^)-projection plane and its image in the (x^x^)-plane, which we name 
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xq- Secondly, there is the origin of the (x^x^)-plane and its preimage in the (x^x^)-plane, 
which we call xo- Finally, there are the two fixed points xpi and xf2 of the transformation, 
which are the stereographic projections of the antipodal intersection points of the sphere 
with the rotations axis into both projection planes: 



(0,0)^^0 = 
, '2asin6'o COS0O 2a sin 6^0 sin 0o \ , . 

1 — cos fc'o 1 — COS Oq I 



2a cos 6q cos (f)Q 2a sin cos (f)Q 



1 — cos 6*0 1 — cos ^0 



xpi = (2a sin 00, — 2acos0o) < — > xfi = (2asin0o, —2a cos ( 
xf2 = (—2a sin 00, 2a cos 0o) < — > xf2 = (—2a sin 00, 2a cos ( 



We choose the two fixed points zi = zi = {x^^ + ixf^a) ^^cl Z2 = Z2 = {x\2 + ^^^^2)^ 
as well as the origin of the (x^x^)-plane ^3 = and its image 23 = (2asin6'ocos0o/(l + 
cos6'o) + z2asin6'osin0o/(l + cos 6*0)). This leads to the following representation of the 
Mobius transformation: 

~_ 2: - 2acot(6'o/2)(cos0o + zsin0o) 
^ ~ (l/2a) cot(6'o/2)(cos0o - '^sin0o)^ + 1' 

5.1.7 The Metric Transformation 

Having derived the coordinate transformation from the stereographic projection plane of 
the south pole to an arbitrary projection plane and back, we now turn our attention to the 
way the metric behaves under such transformations. In general, the metric transforms as 
a covariant tensor field of rank two, 

9ki{x\x') = ^^^ g,,{x\x'). (5.63) 



For the coordinate transformation (B.4) the Jacobi matrix elements are given by (B.7). 



Recalling the metric in terms of the (x^x^)-coordinates (5.8), which were the coordinates 



of the stereographic projection plane around the south pole, 

2a\' 



9vi^\x')={-j (5.64) 



we can now use (5.63) to compute the transformed metric in the projection plane specified 
by 9o and 0o: 



(9a;* dx^ f2a\'^^ (2a 



6,, = (5.65) 



This result is not unexpected, since the plain two-sphere is invariant under spatial rotations. 
Hence, applying the averaging process to the origin of the transformed coordinate system 
is analogous to averaging in the original coordinate system and yields the averaged metric 



(5.66) 
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For the transformation of the averaged metric back to the original coordinate system 
we need to assume that it transforms as a tensor, which is justifiable considering the 
general covariance of the averaging process. Keeping in mind that in terms of the x^x'^- 
coordinates the origin of the transformed coordinate system (5;^,x^) = (0,0) is given by 
(a;^,x^) = {x\,Xq), the transformation law can be expressed as 



dx^ 



dx"- 



dx^ 



(5^m(0,0)) 



(x'l,x2)=(xj,xg) 



(5.67) 



The required Jacobi matrix elements are obtained from the back transformation (B.6) 
by partial derivation. The explicit results are given in (B.5). Inserting them along with 
the averaged metric at the origin of the transformed coordinate system (5.66) into the 
back transformation formula (5.67) returns the averaged metric at the point {x\,x1) = 
(2a sin 6*0 cos0o/(l ~ cos 6^0), —2a sin Oq sin 0o/(l ~ cos6'o)), 

64a^ 



(L^ — (8a2 — L^) cos 6q + 4a sin 9o{x^ cos (f)o + x"^ sin 0o))' 
^(1 - cos9of6ij. 



(x^ ,x'2)={xl,xl) 



Sij 



Doing the same for all points of the sphere leads to the assembled averaged metric in the 
(x^x^)-plane, 

{g,,{x\x')) = ]{l-cose)%. (5.68) 



This result can be best interpreted when expressing 6 in terms of x^ and x^. From (5.1) 
we get 



6 = arctan 



r 



(5.69) 



and from the stereographic projection formula (5.56) we know that 



z = a \ 1 — 



(5.70) 



Combining these coordinate transformations and inserting them into (5.68) leads to the 
complete averaged metric in the stereographic projection plane of the south pole. 



\9iji-^ ) )) 



2a 

T 



Sij- 



(5.71) 



Comparing this to the original unaveraged metric (5.8) leads to the conclusion that the av- 
eraging process preserves the plain two-sphere. Considering that this is the two-dimensional 
space of constant positive curvature, this is in complete agreement with our requirements 
on a reasonable averaging process which would be expected to leave smooth manifolds 
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invariant. 



An obvious question to ask is whether this result also holds for the space of constant 
negative curvature, namely the hyperbolic plane [H"^), and the space of constant zero 
curvature, the fiat plane (H^). The investigation of this will be the topic of the next 
sections. 



5.2 Averaging the Hyperbolic Plane 

5.2.1 The Stereographic Projection of the Hyperbolic Plane 

In contrast to the two-sphere, the two-dimensional space of constant negative curvature 
cannot be embedded in a three-dimensional Euclidean space. While this space can be 
embedded in a five-dimensional Euclidean space this does not help its visualization. On 
the other hand, it can be embedded in a three-dimensional Minkowski space with the 
line element ds"^ = dx^ + dy^ — dz"^, where it becomes a hyperboloid of two sheets. The 
parameter representation of the embedded hyperbolic plane is given by 

X = a sinh 6 cos 0, 

y = a sinh 9 sin 0, (5-72) 
z = a cosh 6^, 

where a =const. The associated line element has the form 

ds^ = a^de^ + sinh^ ed(f)\ (5.73) 
We will apply the analogue of the stereographic projection we used in the case of the two- 



sphere, illustrated in Figure 5.1 Concentrating on the upper sheet of the hyperboloid, we 



project it from the maximum point of the lower sheet onto the plane which intersects the 



minimum point of the upper sheet. This is illustrated in Figure 5^ The corresponding 
coordinate transformation is 

a = 0, 

'0\ (5.74) 



2a tanh 



and, in terms of the new coordinates, the line element is 

''^'^(S;^**' + (5,75) 

Since the angle a is undetermined for r = 0, we again convert to Cartesian coordinates. 



= r cos a, 
= r sin a. 



(5.76) 
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Figure 5.3: Stereographic projection of the hyperbolic plane into the (x^x^)-plane 



which yields the line element 



(4a2 - (xl)2 - (x2)2)2 



The metric and its inverse are then the same as for the two-sphere (5.8) 

4 

9i3 



^ I Sij and 



2a 



except that here is defined as 



4a' 



The Christoffel symbols associated with the metric (5.78) are then those for the 
(5.10) apart from the new definition of and a minus sign, 



and the geodesic equation is hence of the form 



dz^ dz^ 



dr"^ ' dr dr "'''^ dr dr 



2z'gjk- 



dz^ dz^ 



0. 



(5.77) 

(5.78) 

(5.79) 
sphere 

(5.80) 
(5.81) 
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For the initial conditions z*(0) = and dz'^/dr^O) = const, the solution to this equation is 
given by 

^Xr)=2atanhf;^)^(0). (5.82) 
\2a/ ar 

Again, we find that if r is equal the arc length, the initial conditions have to fulfill the 
identity 

'Sm)\[^iO)pi, (5.83) 

and we introduce an angle 7 such that the initial conditions for the geodesies are given by 

— (0)=cos7 and — (0)=sin7. (5.84) 
CLT dr 

Thus, the geodesic coordinates, which we need in order to carry out the integral in the 
averaging process (5.51), are given by 

= 2a tanh ( — ) cos 7 

(5.85) 



2^(r, 7) = 2a tanh cos 7 

2;^(r, 7) = 2a tanh sin 7, 



where < 7 < 27r and < r < 00. 



5.2.2 The Result of the Averaging Process 



To apply the averaging process (5.51) we still need to determine the connector (5.17) and 
the maximally smooth dyad field. In terms of the geodesic coordinates (5.85) the hatted 
connector for the projected hyperbolic plane is 

T 



' 1 ' ) >-'Tuy — cosh 

\2a 



V^/(0,r;C,, 



(5.86) 



and 



dn 



on dK 



(5.87) 



To calculate the maximally smooth dyad field we make the same ansatz as for the two- 
sphere (5.23) except with the new definition of L^. Varying the action gives again the 
differential equation 

{didk(j)) 5'^ = in the area R 
where Or is the geodesic region of radius r, 

^^(7) = 2a tanh 
^^(7) = 2a tanh 



2a 

r 

2a 



cos 7, 
sin 7. 



(5.88) 



The equation (5.87) is again solved by 0(x^,x^) = 0. In terms of the geodesic coordinates 

(5.89) 



(5.85) the maximally smooth dyad field is thus given by 

E^(7,r) = cosh2(^^) 5\ 
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Inserting this and the hatted connector (5.86) into the averaging process (5.51), and ex 



pressing the other quantities in terms of 7 and r as well, the averaged dyad field at the 
origin becomes: 



JO 



27r 



(E",(0,0))=/ / /(^(7,r))t//(;^(7,r))E^(7,r)V^(^drd7 = 5J, (5.90) 



where we used again the fact that the averaging function is normalized to unity over all 
space. Hence, the averaging process returns the following averaged metric at the origin: 

(^7,,(0,0)) = (E^(0,0))(E'',(0,0)>5,fe = 5,,-. (5.91) 

As before we need the coordinate transformation to the stereographic projection plane of 
any other point of the embedded hyperbolic plane and back, as well as the metric trans- 
formation law. 

Rotations that map the hyperbolic plane into itself are different from the rotations that 
map the two-sphere into itself. Having embedded the hyperbolic plane into a three- 
dimensional Minkowski space, we have to account for the minus sign in the line ele- 
ment ds^ = dx^ + dy^ — dz^. Therefore, rotations around the z-axis are the same as 
in the case of the Euclidean space, but rotations around the x-axis and y-axis have to 
be replaced by hyperbolic rotations. The rotation matrix D, which rotates the point 
xo = (a;o,Z/o,-2o) = (asinh6'ocos0o,osinh^osin0o,«cosh6'o) to the minimum point of the 
upper sheet of the embedded hyperbolic plane (0, 0, a), is again constructed from the three 
rotation matrices DzDyD~^. The first matrix, Dz, rotates the hyperbolic plane around 
the 2;-axis until the point Xq intersects with the (x2)-plane. The second, Dy, rotates the 
hyperbolic plane around the ?/-axis until the point D~^xo aligns with the minimum point, 
and the third, D~^, rotates the hyperbolic plane back to the original orientation. Hence, 
the rotation D = D^DyD^^ is 

0\ /cosh 6^0 sinh^^oX / cos0o sin0o 0\ 
1 -sin0o COS0O . (5.92) 




sinh^o cosh^o/ \ 



In this case, the inverse stereographic projection from (x^,x^) to the coordinates of the 
hyperbohc plane x = (x, y, z) is given by 



4a^x^ 4a^x^ / Sa^ 



and the stereographic projection of the rotated hyperbolic plane with the coordinates 
{x,y,'z) into the tangent plane of the new minimum point with the coordinates (x^,x^) is 

~i 2ax _2 2ay 
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The complete coordinate transformation from the coordinates to the coordinates 

(5;^,x^) and vice versa is given in appendix [Cj 



In terms of the coordinates {x\x ) the origin of the (x x )-plane is 



(Xq, Xq) 



2atanh ( j cos0O)2atanh (-^ ) smcpQ 



(5.95) 



Since the rotation axis associated with D is in the xy-plane and therefore does not intersect 
with the embedded hyperbohc plane, the fixed points of the coordinate transformations 



(C.4) and (C.6) are not as easy to find as in the case of the sphere. We calculate the 



Mobius transformation by computing (5.60) and comparing the coefficients. In this way 
we find the the Mobius transformation 



z 



z — 2atanh(^o/2)(cos0o + i sin0o) 
{I /2a) tanh(6'o/2)(— cos 4>o + i sin^o)^; + 1 



(5.96) 



For the transformation of the metric we use the Jacobi matrix elements of the coordinate 



transformation (C.6), which are given in (C.7). Inserting them along with the metric (5.78) 



into the transformation law, we find the metric in terms of the coordinates (x^,x^). 



gki{x^,x^) 



dx^ dx^ ( 2a 



dx^ dx^ V L 



2a 
L 



Jkl- 



(5.97) 



As for the two-sphere the metric is invariant under the coordinate transformation and we 
already know the result of the averaging process in the rotated coordinate system: 



{gki{0,0)) = Ski. 



(5.9^ 



For the inverse transformation of the averaged metric into the original coordinates (x , x ) 



we use the Jacobi matrix elements of the coordinate transformation (C.4), which are given 



in (C.5). Inserting them along with the averaged metric at the origin of the transformed 
coordinate system into (5.67) returns the averaged metric at the point (xq,Xq), which is 



given by (5.95) 



[giji^Q, Xq)) 



64a^ 



(L^ + (8a^ — L^) cosh 6*0—40 sinh 6q{x^ cos (po+x"^ sin 0o))" 
^(l + cosh^o)%. 



(xi,x2)=(xl,a:2) 



Therefore, the complete averaged metric in the (x^x^)-plane is 

1 



[gij{x\x'^)) 



4 



1 + cosh0)%. 



(5.99) 
(5.100) 
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For the correct interpretation of this result we express 9 again in terms of and by 
first transforming it to x,y and z, which follows from (5.72) , 



9 = arctanh 



(5.101) 



and then applying the z-component of the stereographic projection 



z = a 



1 . 



(5.102) 



The result is 



[gij{x\x^)) 



2a 

T 



(5.103) 



This equals the unaveraged metric (5.78 ). We conclude that the averaging process preserves 
not only the plain two-sphere but also the hyperbolic plane. 



5.3 Averaging the Flat Plane 

For the sake of completeness we present the case of the flat plane (H^) explicitly, although 
it follows from the two previous cases for a — > cxd. Due to its simplicity the fiat plane does 
not need to be embedded in any three-dimensional space. The line element is 

ds^ = {dx^f + {dx^f, (5.104) 

and the metric and its inverse are 

gi,=5ij and g'' = 5'^, (5.105) 
where i,j = 1,2. All Christoffel symbols vanish and therefore the geodesic equation is 

dh' 



dr^ 



0. (5.106) 



For the initial conditions z\0) = and dz'^/dr^O) = the solution to the geodesic equation 
is the bundle of straight lines through the origin. 



dr 

The parameter r is equals the arc length if 

dz^ \ ^ / dz'^ ^ ^ 



z-|t) = ?1(0)t. (5.107) 
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and we again introduce an angle 7 such that the initial conditions for the geodesies are 
given by 

— (0)=cos7 and — (0)=sin7. (5.109) 
dr dr 

The geodesic coordinates are given by 



z^{t, 7) = r cos 7 
2;^(r, 7) = rsin7, 



(5.110) 



where < 7 < 27r and < r < 00. In terms of these coordinates the maximally smooth 
dyad field and the hatted connector are 



E\{j,t)=5\ and V/iz{^,T)) = 5]. 
Inserting them into the averaging process results in the averaged dyad field, 

(E%.(0,0))=5%-. 
Therefore, the averaging process returns the metric 

{gij{0,0)) = (E'^,,(0,0))(E^(0,0)>5,f. = 5,, 



(5.111) 



(5.112) 



(5.113) 



at the origin. The coordinate transformation is a pure translation from the plane with 
the coordinates (a;^,a;^) to the plane with the coordinates (5;^,x^), and the reference point 
(xJ,Xq) in the untransformed coordinates becomes the origin (0,0) of the transformed 
coordinates. Therefore, the coordinate transformations are 



X 



X 



X ~\~ Xq^ 
2 _|_ 2 



X 



X 



0; 



X 



(5.114) 

This can be expressed as a Mobius transformation for and z = xl + ix^, 

(5.115) 



1- Z±(Xn + ixVj 



0-z + l 



The Jacobi matrix elements for the two transformations (5.114) are given by 
dx^ 



dx^ 
dx^ 

dx^ 



1, 
0, 



dx 

dx'^ 
dx^ 

dx'^ 



0, 
1, 



dx^ 


= 1, 


dx^ 


dx^ 




dx"^ 


dx"^ 


= 0, 


dx"^ 






dx^ 


dx'^ 



0, 

1. 



(5.116) 



From the metric transformation law (5.63) we find the metric in the x^x^-coordinates, 

gki(x'^,x'^) 



dx^ dx^ 



dx'' dx 



-^gki{x^,x^) = 5ku 



(5.117) 
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and from the back transformation (5.67) we get the averaged metric at the point (xJ,Xo): 



dx'- 



(a;l,a;2)=(a;l,a;2) 



(5.118) 



The complete averaged metric in the coordinates {x^^x^) equals the unaveraged metric 

{gij{x\x^)) = 5,j. (5.119) 



(5.105), 



In summary, we have shown in this chapter that the averaging process preserves all two- 
dimensional spaces of constant curvature. The obvious question to ask now is what effect 
the averaging process would have on a space of contant curvature, say a two-sphere, if we 
added a small perturbation. This will the topic of the next chapter. 



CHAPTER 6 



THE PERTURBED TWO-SPHERE 



In the previous chapter we tested our averaging formulation on the conceptually trivial 
examples of two-dimensional spaces of constant curvature and verified that it leaves such 
spaces invariant. In this chapter we extend our discussion of the averaging process in two 
dimensions to the more interesting case of a sphere whose metric is slightly perturbed and 
derive the expressions for a general perturbation function. A specific example is discussed 
in chapter |8| 



6.1 The Perturbed Metric 

To begin, we recall that embedding the unperturbed sphere into a three-dimensional Eu- 
clidean space enabled us to specify the stereographic projection from the Euchdean coor- 
dinates of the sphere x = {x, y, z) to the coordinates of the projection plane (x^, x^). The 
associated coordinate transformation (5.56) was 

X = Aa^x^/L^, 

y = Aa^x^L^, (6.1) 
z = a{l-8a^/L'^), 

where a denotes the radius of the sphere and is 

L^ = {x'y + {xy + 4a^. (6.2) 

We now apply a small perturbation defined by 



Xp = {1 + r]f)x = (1 + rjf)Aa^x'^/L 



2 

1 

2/p = (1 + r]f)y = (1 + r//)4aV/L2, (6.3) 
Zp = {l + r]f)z = (1 + vfH^ - Sa'/L^), 
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with a space-dependent perturbation function / = and a small parameter 77 > 0. 

Up to first order in rj the hne element takes the form 

,2 ,2 ,2 ,2 / 16a^ 32a^/\ . , /16a'' 32a''/\ . , 2\2 

ds' = dxl + dyl + d^J = / _ + 77-^ J (da;^) + f — + ^^^1 (^a;') , (6.4) 

and hence the metric and its inverse are given by 

{gp),, = (1 + 2r^/) 5., and {gpf = (1 - 2r//) (6-5) 

where the indices take the values i,j = 1, 2. 



6.2 The Perturbed Geodesies 



In the case of the unperturbed two sphere the solution to the geodesic equation for the 



initial conditions z'^{0) = and dzydT{0) = const. (5.12) was given by 



.•M = 2atan(^)^(0). 



(6.6) 



Since the imposed perturbation accounts for a change in the arc length of the geodesic, 



we are forced to use the general geodesic equation (A.7), where the parametrization r 



of the geodesic is not necessarily the arc length s. This equation is invariant under a 
reparametrization of the form r' = r + t]w{t) with a continuously different iable function 
w{t). That allows us to make the following ansatz for the parametrization of the perturbed 
geodesic: 

dr^ 

z},{T) = z\r)+r]viT)e,, — iT). (6.7) 

Here v{t) is a twice continuous differentiable function which describes the longitudinal de- 
viation of the geodesic caused by the perturbation. This parametrization is also visualized 



in Figure 6.1 The arc length of the perturbed geodesic up to first order in rj is given by 



(6.8) 



with the inverse function 



r{s) = s-7^ / f{s')ds'. 
'0 



Inserting the arc length (6.8) and the derivative of (6.7), 

f(r)^^(.)^-.|^(.)e^f(.) + ..(.)e^^(.). 



(6.9) 



(6.10) 
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Figure 6.1: Parametrization of the perturbed geodesic. 



into the equation for the perturbed geodesic (A. 7) yields the differential equation of a 
driven harmonic oscillator for the function v(t). 



d'^v{T) v{t) h{T) 



dr' 



cos 



Here, the driving force h{T) has the dependence on the perturbation function /(r) 



Mr) 



df 



34 
dr 



(0) 



df 



(zl,x2) = (2l(T),z2(r)) 



04 

dr 



(0). 



(6.11) 



(6.12) 



The initial conditions are chosen such that the perturbed and unperturbed geodesies and 
their derivatives coincide at the origin, 



zUO) = z\0) 



and 



dzl 
dr 



(0) 



dz' 
dr 



(0). 



(6.13) 



Therefore we are looking for solutions to the differential equation (6.11) for the initial 
conditions 

dij 

v{0) = and ^(0) = 0. (6.14) 

dr 



The general solution of the differential equation (6.11) is composed of the solution of the 



homogeneous equation and a particular solution of the inhomogeneous equation, which can 
be found by variation of constants. It is given by 



f (r) = A cos ^- j + 5 sin ^ 



a 

a (T 
- cos - 
2 \a 



sm(^ 



—^dr' + - sm - 

cos^r^) 2 \a 



h{r') 



C0S(^)C0S2(I^; 



(6.15) 
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with A,B = const. For the initial conditions f(0) = and dv/drifi) = the constants A 
and B vanish. Hence, the solution is 



a /r\ r 1 /i(r' 

v{t) = — cos - / -, 

^ ' 2 \aJ Jo sin(^)cos2( 



^dr' + -sin -) / , \,L dr'. (6.16) 

^) 2 \a) J, cos^)cos2(i^) 



6.3 The Maximally Smooth Dyad Field 



Here we perform a calculation analogous to that in section [5.1.3 In this case, we will 
derive the differential equation needed to determine the maximally smooth dyad field for 
a general metric and then specialise to the case of the perturbed sphere. 



We decompose the metric into dyads according to 



(6.17) 



We choose a reference dyad field E^i{x) and express the class of tetrads that obey (6.17) 
as a rotation by a position dependent angle 0(x^,x^), 



E\{x) = U\{x)E\{x) where U\{x) 



cos (p{Xj smcp^xj 
— sm.(hix) cos (hix] 



The partial derivative of the rotation matrix Uab is again 

d,U\{x) = e\U\{x){d4{x)), 
The generally covariant derivative of the dyad field then becomes 



(6.18) 



(6.19) 



D,E% = {d,U\)E', + U\{D,E',) = e\UMd4)E', + U\{d,E' ,) - U\r^,E\. (6.20) 



The Lagrangian for the maximally smooth dyad field (5.33) is 



(6.21) 



Cms = {D,E^j){DkE\)g'^g^'5ab = e\U\E'^e' jU^ d{d4){duct))E''ig'''g^'5ab 

+ 2U\{DiE',)t' fUf d{dk4>)E\g'''g''6ab + U\U\{D,E'^){DkE\)g"'g^'5ab. 

This can be simplified further by applying the identities (5.30 15. 3^ ) and (6.17), which leaves 

Cms = 2{d4){du<P)g''' + 2{dkct>){D,E'^)&ie,dg'^g^' + (A^^)(/^fcM) W/- (6.22) 

We define the vector field 

m'^ = {D,E',)ESe,dg'''g^' (6.23) 

and set the variation of S* = / d'^x^ Cms with respect to equal to zero. As shown in 
section 5.1.3 , the result is a differential equation on i?, 

1 



(6.24) 
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with boundary conditions of Neumann type on dR 

d(f) 1 



dn 



(6.25) 



To investigate whether the solution 0(x^,x^) to this differential equation depends on the 
choice of reference dyad field E^i{x), we choose a different reference field E'^i{x) and find 
the differential equation 



9 



ik 



did) 



(6.26) 
(6.27) 



with the boundary condition 

2' 

Since the reference field E'^i{x) can be expressed as a rotation of the reference field E^i[x) 
by a position dependent angle a(x^,x^), 



E^{x) = U% {a{x))E\{x) 



(6.28) 



the vector field becomes 



u 



{D,E'^^)E'ieUgpyH9pf 

{e%UMd,a)E'j + f/=b(A^^)) U'^fEfieUgpY^gpf 
2{gpy\d.a) + {D,&,)&ie,f{gpy\gpf 



(6.29) 



When inserting this in the differential equation (6.26) and the boundary condition (6.27) 



we can add the a-dependent term to the left-hand side, 

1 



and 



r -9A(0 + «; 



d{(P + a) 



(6.30) 



(6.31) 



dn 2 

The maximally smooth dyad field is a rotation of the reference field E'^i{x) by the solution 
0(x^,x^) to this differential equation. Since E'^i{x) is a rotation of the dyad E^i{x) by the 
angle a(a;^, x^), we can also express the maximally smooth dyad field as a r otati on of E\{x) 
by the angle = (0 + a). The differential equation is then identical to (6.24) and (6.25[) 



and we conclude that the definition of the maximally smooth dyad field is independent of 
the choice of reference dyad field. 



To find the maximally smooth dyad field for the perturbed sphere with the metric (6.5) 



we first need to calculate the vector field (6.23) 



u 



{D,E',)E\e,d{gpy\gp 



(6.32) 
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In the two dimensional case with the general metric 
we can always choose the refernece field 



9ii 912 

921 922 



V911 - (912/922) 




and 



9l2/y/922 
\f922 



For gu,g22 > and g = det{gij) > this dyad field is well-defined and obeys 

E°'iE'^j5ab = gij- 

For the perturbed sphere we therefore choose the reference field to be 



E\ 



{Aa'/L^) + r/(4aV/^') 




and 



E' 







(4aVL2) +r7(4a2//^' 



Then the vector field (6.23) becomes 



u 



2a4 ^Sa^dx^ 



and u 



+ r] 



L^fx^ df 
2a^ ^ Safari 



(6.33) 



(6.34) 



(6.35) 



(6.36) 



(6.37) 



The covariant derivative of this vector field vanishes, DkU^ = 0, and thus the differential 



equation (6.24) reduces to 

(gpY' (OA 



+ 



(dx^y {dx"^) 



0. 



Thus, the differential equation for 0(x^,x^) on R is 

'{x\x^) = A(f){x\x^) = 0. 
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+ 



(dx^y (dx 



2\2 



(6.38) 



(6.39) 



Recalling the initial conditions for the perturbed geodesies v{0) = and dv/drifi) = 0, we 
can introduce an angle 7 with < 7 < 27r in the same way as for the unperturbed sphere 
(pl. 



dzp 



(0) = cos 7 



and 



dzp 



(0) = sin 7. 



(6.40) 



dr dr 
That means in particular that all functions which depend on the initial conditions acquire 
a dependence on 7 which must be included in the argument. We will therefore from now 
on denote v{t) and h{T) by f('r, 7) and h{T,'~f). The geodesic coordinates are then given 
by 

t;(r,7) . 



/ T \ V(T 7) 

2p(r,7) = 2a tan — cosj + i] — sin7 
V2a/ cos^(^) 

Zp{T, 7) = 2a tan (^^^ sin 7 — 77 



(6.41) 



cos 7. 
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The geodesic region R is defined as the region inside the entirety of points at geodesic 



distance r from the origin that form dR. The geodesic distance is the arc length (6.8) 
along the geodesies. In parameter representation the boundary dR is 



/ T \ COS ^ f 

a^(7) = 2atan(-jcos7-r,^^y^ 
a^(7) = 2a tan ^ 



COS7 r t;(r,7) 
V ^TTT / f{sn)ds +r/— :j7^sm7 



2a' 

sin 7 

2a J ^ ^cos^(;^ 



cos2(^) 



(6.42) 



cos 



COS 7. 



An important quantity for the calculation of the Neumann boundary condition is L^. Its 
value at the boundary dR is 



4a2 



4atan( 



cos2(^; 



V- 



2a' 



COS 



2( JL 



^2a! -JO 



(6.43) 



and hence the square root of the determinant of the metric at the boundary is 



f{s\^)ds'. (6.44) 



From this we can compute the generally covariant normal vector Vk to OR, 

I (9a' 

the normal unit vector Uk to dR, 

. . 1/2 

nk = {{gpy^{r,-f)uiUj) z/fc, 



(6.45) 



(6.46) 



and the value of the vector field u , given in (6.37), at the boundary. The Neumann 



boundary condition on dR can then be calculated from (6.24), giving 

/i(r,7) 



dn 



V 



COS' 



(2a) 



+ r] 



2a2 



cos^ 



(2a) 



d^ 



f{s',l)ds' 



(6.47) 



In summary, the maximally smooth dyad field for the perturbed sphere is given by the 



(6.47) on dR, where dR is given by (6.42) . 



solution to the differential equation (6.39) on R with the Neumann boundary condition 



6.4 The Perturbed Connector 



Using that 



cos 



r 
2a 



(6.48) 
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and the perturbed geodesic (6.7) we find the Christoffel symbols 



dzp 
~d7 



(6.49) 



f df ^ df ^ df ^ 

^ox ox-' ax J (^i^^2)=(2i(^_^)^^2(^^^)) 

The connector of parallel transport along the geodesic 2;p(r, 7) from the origin to the point 
with parameter r is then 



y(r,0;Co.) 



Pexp I — / - 

Pexp ( ( - / tan ( — ) dr' 
\\aJo \2aJ 



r',7) r, {zp{T',^))dT' 



, ,\ (1 



1 



+ Tj 

1 



1 r , 1 



2a2 



cos2(|^; 



(6.50) 



dr' + — / tan^ j^— J t;(r',7)rfr' 

fr'\ dv, , , , , r /i(r',7) , A / 1 
+ - / tan — — (r',7)dr' + / \ ' ,^ dr' ^ ' 



\2a J dr 



C0S2(I^) 



-1 



The integrand in (6.50) is of the form Mo(r') + riMi{T') with two matrices Mo(r') and 
Mi(r'). Defining the matrix Uq to be 



Unir) 



(6.51) 



we can expand the general identity 
Pexp ( fdr' {Mo{t') + vM,{r'))] = U^{t) V exp (r, rdr'U^\T')Ivh{r')U,{T') ) (6.52) 



to first order in rj and find 
Vexp(^j^dT'{M,{T')+r^M^{T'))^ = Uo{t) (^1 + r] j^dT'U^\T')M^{T')Uo{r')^ . (6.53) 



For the connector (6.50) the matrix Uq{t) is given by 

Uq{t) =Vexp(^j rfr Mo(r')^ = V exp (^(^- tan (^^^ dr'^ 1 



"Pexp 



—2a In ( cos ( — 
\2a 



(6.54) 
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and therefore the first order expansion of the connector becomes 

'1 



V{t, 0; Co.) = cos-2 (^) (1 - Vfir, 7) + r]/(0, 7)) (^q J 

Using (6.11) and partial integration methods this expression can be further simphfied for 
the boundary conditions ^(0,7) = and c?t>/(ir(0, 7) = 0: 

Vir, 0; Co.) = cos-^ (^) (1 - r//(r, 7) + vf{0, l)) 

+ ,(.'(r,7) + ^tan(^) 



1 
1 



i;(r,7) 



(6.56) 

or' 

4 



For the averaging process we need the connector from the point with parameter r to the 
origin and hence we need the inverse of this expression to first order in 77, which is 

l^(0,r;C.o) = V-'(r,0;Co.) 

(^)(l + 2^/W-2^/(0)))(cos"^ {^)(^-^f('')+^f(^)))(l ? 



cos^(;^)(l+2r//(r)-2r//(0)))(^r^cos-2(^)(^^;'(r) + ^an(^)^;(r)^^^ ° ^ 



.2aJ^ ' "Jy \2aJ\ ^ ' a \2aJ ^ J J\-^ Oy 

In components the resuh can be written as 

V^^(0,r;C.o) = cos2 (^) ((l+r^/(r) -r^/(0)) 5^ - [v'{r) + ^ tan (^) v{r)^ e^^ . 

For the parallel transport of dyads in the averaging process the hatted connector is needed 
and we find it to be of the form 

V/(0, r; Co) = gA^, 7)^^(0, r; C.o)/(r, 7) 

= cos-2 (^)((l+^/(0,7)-r//(r,7)) 5]-r^ {v'{r,i) + ^tan (^) v{t,^)\ e^,") . ^^'^^^ 



6.5 Alternative Computation of the Parallel Transport 

When the vector E°'i is parallel transported along the geodesic 2*(s,7), the angle between 
the vector and the tangent vector of the geodesic is preserved. We can use this to compute 
the parallel transported dyad at the origin in a direct way without the connector. 
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Figure 6.2: Definition of angles in tlie parallel transport of the dyads. 



Recall that we expressed the maxir nally smooth dyad field E'^i in terr as o f the reference 
dyad field E^i, which was given by (6.36). In terms of the arc length (6.9) this reference 
field is 



g".(..,) = ^",oos^(A) (l + ,/(„) + „lta,> (A 



(6.58) 



The maximally smooth dyad field is a rotation of this reference dyad field by a position 
dependent angle (l){x^,x^), 



E^, 



U\E\ 



where 



cos (p sm ( 
- sin d) cos ( 



(6.59) 



Let the angle between E°-i{s, 7) and the tangent vector of the geodesic zp{s, 7) be denoted 
by /3, as shown in Figure [6l2j After the parallel transport the angle between the transported 
dyad i?"j(0,7) at the origin and the geodesic zp{0,'y) is still Since the angle of the 



geodesic at the origin was, according to (6.40), given by 7, we can express the parallel 
transported dyad as a rotation of the reference dyad field E^i by the angle 



0-/3 + 7- 



Since 



7 = arccos 



dzp 
ds 



(0,7) 



and 



(3 = arccos ( -^(s^'l) ] > 



(6.60) 
(6.61) 



we need the perturbed geodesic (|6.7|) in terms of the arc length (6.9), which is 

cos 7 f-r r \ 1 / , sin7 



Zp{s, 7) = 2a tan j cos 7 — 
Zp{s, 7) = 2a tan j sin 7 — 



sin 7 



cos' 



.2a' JO 



f{s\-i)ds' + riv{s,-i)- 
f{s',-f)ds' - r]v{s,-f) 



cos 



2f_s_ 



2a' 
COS 7 



(6.62) 
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From the derivative of this and the identity 



dz p , , 



(6.63) 



we find 

COS0 = cos(</) — /3+7) = cos (/) ( cos 7 cos + sin7sin/3) — sin(/)(sin7cos/3 — sin/3cos7) 



cos ( 



cos ( 



^(0)^(s) + ^(0)^(s) 
ds ds ds ds 



-sm < 



dz p dz p ^ 



dzp , dzp 



cos 



/(0,7) /(5,7) 



ds 
1 sin 



ds 



ds 



ds 



(0) 



sin ri 



2(J_ 

1 dv 



COS 



2(_s_ 



V- 



2a' 



COS 



2( A_ 



COS^f;^-) ds 
^ la' 



(s,7) +r/ 



2a' V2a 

t;(s,7) sin( 



a cos 



^2al 



2a' 



a cos 



2( 

~2a' 



(6.64) 



and 
sin I 



sin(0— /5+7) = sin 0( cos 7 cos (3 + sin 7 sin /3) + cos 0( sin 7 cos /5 — sin/5 cos 7) 
1 /(0,7) /(s,7) 1 sin( 



smi 



cos 



~2a' 



V 



cos2(^) '^cos^ 



s ' 
2a' 



a cos 



^2a' 



.2a^ -'O 



f{s',i)ds' 



+ COS ?7 



cos2(^) (is 



(it; , , v{s,'y) sin 



2a' 



a cos 



2('^-^ 
^2a' 



(6.65) 

The rotation U{4>)/E°-i{s,'y) of the reference dyad field, (6.58), does not yet give the 
correct result. The reason is that the length of the vector E'^i is changed when it is parallel 
transported. We can calculate the rate of change from a quantity which is conserved 
during the parallel transport, and this is the scalar product of the transported vector and 
the tangent vector of the geodesic. 



i?^(s,7)^(.,7) = ^^(0,7)^(0,7). 
ds ds 



(6.66) 



From this we find that to first order in rj the factor of change is given by 

C = l + 2r//(0,7). (6.67) 
Altogether, the parallel transported maximally smooth dyad at the origin is given by 

E'^,iO,^) = CU/i^)E%sn) (6.68) 

with the components 

dv 



E\{0, 7) = ^2^(0) = COS0 (1 + 77/(0)) - sin0 
E\{0,j) = -E\{0,j) = sin 0(1 + r//(0)) + cos0 



ds 
dv 
ds 



[s) + r] 



v[s) 




r s 




tan 




a 




V2a 


v{s) 




r s 




tan 




a 




V2a 



(6.69) 
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This result is identical to the connector result, 

E%(0,7) = i>/(0,s;C,o)^^(s,7), (6.70) 



with the connector (6.57) and the maximally smooth dyad field (6.59) expressed in terms 
of the arc length s. 

6.6 The Averaging Function 

The weighting function f{xQ,x';Cxixo) in the averaging process is required to be positive 
and normalised to unity, 

f{xQ,x';Cx'xQ)\/-9{x') d^x' = 1. (6.71) 



R 

The simplest choice is a step function which drops to zero at the boundary dR. With 
such a choice all tetrads of the geodesic region contribute to the average with equal weight 
and furthermore there are discontinuities at the boundary of the region. It is therefore 
preferable to choose a weighting function that has a maximum at the origin and falls off 
to zero smoothly at the averaging radius r. 

The Beta function is defined as 

B{x,y) = j\^^-'\l-t)^y~^)dt = ^1^. (6.72) 

For positive integers m and n it takes the value 

(n-l)!(m-l)! 
B{m, n) = ^ '-^ — ^. 6.73 

We can extend the interval [0, 1] to [— r, r] if we set 

t=i(l + ^) and (1-^) = ^1-^)- (6.74) 

Motivated by the Beta function, reasonable choices for the averaging function are thus the 
normalised functions 

/„(r) = iV^(l-^)^ (6.75) 
where is the normalisation over the geodesic region. We choose the function /2(t), which 



is shown in Figure 6.3 for the averaging radius r = 5, 



4N 2 



6. 7 The Average Metric 



91 




Figure 6.3: Normalised averaging function /2(t) for the averaging radius r=5 in the flat 
two-dimensional case Qij = 6ij. 

For the perturbed two-sphere we find 

N = (^fj h{r)^^)A^)didT^ = lo' - V lo' If. (6-77) 



where we defined 



and 



/o(r,7) 



2iT rr 



^0 



/2(r) cos^ ( — ) d-idT 



2tt rr 



(6.78) 



(6.79) 



With the normalised averaging function /2(t) we have now all constituents at hand to 
apply the averaging process to the perturbed two-sphere. This will be done in the next 
section. 



6.7 The Average Metric 



The general form of the averaging process, (4.45), is 



(E",(x))= / f2ix,x';C,,,)V/ix,x';C,,,)E%ix')^M^) d'x'. 



(6.80) 



In terms of the geodesic coordinates (6.41) the averaging process becomes 



(^^•(0,7)) 



JO 



/2(r)\//(z(r,7))i?'^.(r,7)V^7p(r,7) d^d 



T. 



(6.81) 
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We insert the connector (6.57) and the maximally smooth dyad field, expressed as the 
rotation of the dyad field (6.36), where the angle of rotation 0(r, 7) is the solution to the 
differential equation (6.39) on the area R with the Neumann boundary condition (6.47) on 
the boundary (6.42), 



cos0(r, 7) sin0(r, 7) 



sin0(r, 7) cos0(r, 7) 



27r fr 



Furthermore, we define 

/c(r,7) 
and the scalar fields 

^(r.oO = cos^(^)( 



5'-' (l + r^/(r,7))cos^(^2a 



Jo 

2tt r-T 



JO 



/2(r) sin(0(r,7))cos'' 
/2(r) cos(0(r,7))cos^ 



2a. 
r 
2a 



d'-^dr, 



T \ / 1 T d 

-tan( — )i;(r,7) + — t;(r,7) 
a Za (XT 



S(r,7) = cos^ (^j (/(0,7) + 2/(r,7)). 
The components {E°'j{0)) of the averaged dyad at the origin are then 



Jo Jo 



(cos(0(r, j))B{t, 7) - sin(0(r, 7))A(r, j)djdr) 



(i?S(0,7)> = - (i^'i(0,7)> = lo' Is - V lo' If Is 



+ V lo' 



2-K rr 



Jo 



(cos(0(r, 7))A(r, 7) + sin(0(r, 7))i?(r, '^)d'^dT) 



(6.82) 



(6.83) 



(6.84) 



(6.85) 



The averaged metric at the origin is again diagonal and the metric elements are recomposed 
by the averaged dyad at the origin. 



and 



(Mii(0,7)) = (M22(0,7)) = [{E\{Q,^))f + {{E\{Q,^))y 



(Mi2(0,7)) = (M2i(0,7))=0. 



(6.86) 



(6.87) 



6.8 The Coordinate Transformation 



Starting from the coordinate transformations which we used for the unperturbed sphere, 
(B.4) and (B.6), we have to make two modifications. First, we notice that for a numerical 
calculation of the averaging process on a grid it is much easier to express the reference 
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point Xq = (xo,?/o,-2o) = (a. sin 6*0 COS 00, a sin 6*0 sin 00, « COS 6*0) in terms of the coordinates 
of the stereographic projection plane. It is then given by 



Xn 



Xfi 



2a sin 6*0 cos 0o 
1 — cos 6*0 

2a sin 6^0 sin 0o 
1 — cos 6^0 



(6.88) 



Second, we have to take the perturbation into account. The coordinate transformation 
is composed of three parts: the stereographic projection from the (x^x^)-plane onto the 
sphere, a rotation of the sphere, and a stereographic projection into the tangent (x^x^)- 
plane of the reference point xq. The perturbation enters the coordinate transformation only 
through the two stereographic projections, but it can be shown that the first order con- 
tributions cancel. We thus use the unperturbed coordinate transformation (B.4) and (B.6). 

To find the perturbed metric in the coordinate system of the reference point (xq,Xq) we 
recall that the transformation law of the metric is given by (5.63): 



dx^ dx^ 



dx'^ dx 



(6.89) 



Inserting the coordinate transformation (B.6) and the perturbed metric (6.5) yields 

16a^ 



{gp)ki{,x^,x'^) 



L4 



il + 2r^f{x\x^)). 



(6.90) 



The change to the coordinate system (x^,x^) affects thus only the perturbation function 
/(x^,x^), which has to be transformed to the new coordinates. 



Applying the back transformation law (5.67) to the averaged metric at the origin of the 
coordinate system (x^,x^) we find the averaged metric at the reference point (xq,Xq) in 
terms of the original (x^x^)-coordinates to be 



{i.9p)ij{.xl,xl)) 



dx^ 



dx^ 



dx^ 



(xl,x2) = (xl,x2) 



1 dxi 



(xl,x2) = {4,x2) 



((fi^p)fe/(0,0)). (6.91) 



metric elements. By inserting (B.4) we find 



Since the averaged metric is diagonal we only need the transformation law for the diagonal 

(6.92) 
(6.93) 

T = r^TT^. (6.94) 



and 
where 



((<7p)n(xJ,x^))=T((5?p)n(0,0)) 

((^7p)22(4,X^)>=r((fi^p)22(0,0)), 

16a^ 



(4a2 + {xlY + 



X, 
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The result for the averaged metric is therefore 

(Mi2(0, 7)) = ((^7p)2i(0, 7)) = (6.95) 

and 

(Wii(0,7)) = (M22(0,7)) = ((i?\(0,7)»'+ ((i?S(0,7)»' 

16a^ 



(4^2 + ^2 + ^2)2 



+ Jq"^ Ic / /2(t) (cos 5 - sin A) d-fdr 

Jo Jo 

+ Jq"^ Is / /2(r) (cos A + sin S) (i7(ir 

Jo Jo 



(6.96) 



In this chapter we investigated the averaging process for a general imposed perturbation 
function /(x\ x^) to the two-sphere. We will use the results of this chapter to compute an 
explicit example of such a function in chapter [8j 



CHAPTER 7 



AVERAGING THE THREE-SPHERE 



In the previous chapters we have considered the smooth and the hnearly-perturbed two- 
sphere in some detail. However, despite the successes in those metrics we should demon- 
strate its validity in a higher dimension. To this aim we present a proof-of-concept for 

a three-sphere, first applying the averaging process to the smooth case and then to the 
linearly perturbed case. Such spaces form the basis of closed FLRW models. 



7.1 The Metric from Stereographic Projection 

As in the three dimensional case we start with a stereographic projection of the three- 
sphere into the tangent plane at the south pole. Let (^^, C^, C^) denote the Cartesian 
coordinates in four dimensions and (x^, x^) the ones of the projection plane. The radius 
of the four sphere is denoted by a. Solving the equation 



/0\ 




\2aJ 



+ A 



e 
e 



for A returns the relations 



2a 






V2a/ 



x"^ 
x^ 



(7.1) 



where the index i G {1,2, 3}. To express ^'^ in the coordinates of the projection plane we 
make use of the fact that the surface of the sphere is given by 



(7.3) 



With the quantity 



= 4a^ + {xy + {x'f + {x') 



(7.4) 
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we can express the coordinate as 

^ 2a{{xr+{xr+{x'r) _ ^^^^ 

Prom here we can calculate the line element in terms of the coordinates of the projection 
plane, 

ds' = {der + ider + {der + id^r = ^ Hdx^r + {dxr + {dx^) . m 

The metric and its inverse are therefore 

^.. = (^)'% and 9''=[^Js'', (7-7) 
with the indices i,j e {1, 2, 3}. 

7.2 The Geodesies and the Connector 

As for the two-sphere the Christoffel symbols of this metric are given by 

- = A 



2 



and the geodesic is 

^ , d^z* , dz^ dz^ dz^ dz^ „ r>\ 

For the initial conditions 2;'(0) = and dzydT{0) = const, the solution to this equation is 
again 

^'(r) = 2atanf;^) ^(0), (7.10) 



■.2a/ dr 

with the additional requirement 



The latter identity ensures furthermore that r is the arc length 



The connector in this metric is 



V{T,Q;Cor) = Vexp (-^^cos'^ (^) ^(0) TMr')) dr'^ = cos' (^) 1, (7.13) 
and hence the hatted connector that appears in the averaging process is 

V^/(0,r;C,o) = cos-2(^) S^. (7.14) 



1.3 The Maximally Smooth Triad Field 
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Figure 7.1: Definition of Euler angles 



7.3 The Maximally Smooth Triad Field 

The rotation group 5*0(3) is the group of homogeneous linear transformations 

x"" = R\x^ 



(7.15) 



that preserve lengths and orientation. A general element of the group can be characterized 
by three angles. A common choice of angles are the Euler angles. Let the Cartesian 
coordinates of the unrotated frame be ( ) and the coordinates of the rotated frame 

be ( '^). Let further the intersection of the (x^a;^)-plane and the (x'^x'^)-plane, 



namely the nodal line, be labelled i^', as shown in Figure 7.1[ Then the Euler angles 
{(f), ip, 6) are defined as follows: (f) is the angle between the x^-axis and K, ip is the angle 
between K and the x'^-axis, and 6 is the angle between the x^-axis and the x'^-axis. Any 
rotation can then be represented as a composition of a rotation about the x'^-axis by 
(0 < < 2tt) and a rotation about the by ^ (0 < < 27r) with a rotation about the 

x^-axis hj ip {0 < ip < 2tt). In terms of the Euler angles the parametrization of a general 
rotation is 



R(0, e, ip) = R{ip) ■ R{e) ■ R(0) 

cosip sinip 0\ /l 

- sin ijj cos ijj I cos 6 sin 6 

1/ \0 -sine cos^, 



cos sin 
— sin cos 
1, 



(7.16) 
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cos'0cos0— sinz/'cos^sm^ cos'0sin0+sm'0cos^cos0 sin ^ sin ■0 
— sin'^cos^— cos-i/^cos^sin^ — sin'0sin0+cos-?/'cos^cos0 cos -0 sin ^ | . (7.17) 
sin 9 sin cj) — cos (p sin 9 cos 9 

Now we express the maximally smooth triad field as a rotation of the reference triad 
field E''i. In terms of the Euler angles this is 

E\ = {R^)\{<P, e, ij) E\ = U\{(P, 9, ij) E\. (7.18) 

To determine the differential equations for the Euler angles we insert this into the La- 
grangian 

£ms = {DiE'^^){DkE\)g'''g^'8a,. (7.19) 

With the substitution 

{D,E''j)E\ g'^g^' = {Q'^'^f (7.20) 

we find 

>CMS = 2(9iV')(9fcV')^''+2(9,e)(9fe^)^''=+2(9,0)(9fc)^''^+4cos^(9,V')(9fe0)^^*^ 
-2 ((9feV')+cose(9fe0)) {{^k^^J)+cos9{^k(|>)) {Q^'f 

-2 {s\ni){dk9) -sin 9 cos ^'(9^0)) (sin ^'(9^^) -sin 9 cos ^'(9^0)) {Q^^f (7.21) 

-2 (cos '0(9fc^) +sin 9 sin ^'(9^0)) (g^^)'=+2 (cos'ijj(dk9) +sin ^ sin V'(9jfc0)) (g^^)'= 

+ (A^'^,)(I?fc^'0^''/<^cd. 
Applying the principle of least action on a finite integration volume R with boundary 
dR yields three partial differential equations with Neumann boundary equations. The 
differential equation for on i? is 

'2{Vg)-'dk (v^ g'HdicI^)) +2cos9{^)-'dk g'Hdii^)) -2sm9g'''{di^){dk9) 
= cos9Dk{Q^'^f - cos 9Dk{Q^y - sm9{Q^^)''{dk9) + sm9{Q^y{dk9) 

- sin 9 cos i^DkiQ'^^f + sin 9 cos i^DkiQ^^f - cos 9 cos i^iQ'^^f {dk9) 

+ cos 9 cos V'(9fe^) {Q^y + sin 9 suii^{Q^^f - sin 9 sin il^iQ^^idki^) ^ '"'^^^ 

+ sin ^ sin i/jD^iQ^^)'' - sin ^ sin ijjDkiQ^Y + cos ^ sin 'il){Q'^^f{dk9) 

- cos 9 sin ilj{Q^'^)''{dk9)+ sin 9 cos ilj{Q^Y{dk^) - sin9 cosil^{Q^'^f{dk'il^) 
with the following boundary condition on dR: 

2^ + 2cos9^ = cos9nk{Q'Y - cos9nk{Q^')' - sin9 cos^nkiQ^"^)' 
on on (7.23) 

+ sin6'cosV'njt((5^^)'' - sin6'sinV^nfe(Q^^)^ - sin6'sin^/'nfe((5^^)''. 
Here, denotes the normal unit vector to dR. The differential equation for ^ on is 
2{^)-'dk g'HdiO)) +2sin9g"'{di^){dk<f>) 

= sim/jDkiQ^^)'' - simljDk{Q^^)'' + cos9cosilj{dk(l)){Q^^)'' + costlj{Q^^)\dki^) 

- cos9cosilj{dk(p)iQ^^)'' - cosilj{Q^^)\dki^) + cos^ljDkiQ^^)'' - cosijDkiQ^^)'' (7.24) 
+ cos 9 sin i^{dk<t>){Q'^^f - sin ij{dk^){Q^^)'' + cos 9 sin'il^{dk<t>){Q^'^f 

+ sin (g^^)*^ + sin 9{Q^''f{dkct>) - sin 9{Q^')\dk(f>) 
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with the boundary condition on OR 

an 

2 = sinijUkiQ'^)'' - sm^rikiQ^')' + cos Vnfc(Q23)'= - cos^nkiQ^^. 
on 

And finally the differential equation for ip on R is 

2{^)-'d, {^g'\d,^P))+2cosei^)-'dk g''id,<l>)) -2smeg^'id,(l))idke) 

- sin ^ sin ^(9^0) (g^^)'^ + sin^sin^/'(9fc0)(Q3^)^ + sin^(Q23)fc(5^^) 

- sin ^/^(Q32)fc(5^^) _ sin 6 cos ^(9fc0)(Q23)fc + gin 6 cos ^/^ (9^0) {Q'^^ f 

with the boundary condition on OR 



(7.25) 



(7.26) 



2-^ + 2 cos ^-^ 
on on 



21\k 



(7.27) 



7.4 The Result of the Averaging Process 

To solve the differential equations for the maximally smooth triad field we first need to 
decompose the metric of the smooth three-sphere (7.7) into triads and choose among the 

E^, = {^] 5%. (7.28) 



options one reference triad field to start from. The easiest choice is 

2 



The vector fields (7.20) then become 



8a 



2a 

T 



IX 



(7.29) 



Furthermore, we need the normal unit vector to the boundary of the geodesic region dR. 
The requirement (7.11) allows us to introduce the angles 71 G [0, vr) and 72 G [0,27r), such 
that 

dz^ , , dz"^ , , dz^ 



-(0) = sin7i cos 72, 



-(0) = sin7i sin 72, 



(0) 



COS71. 



(7.30) 



dT dr dr 

Then we can express the boundary of the geodesic region R of radius r in terms of these 
angles, 



2a tan 



2a tan 



2a 
r 

2a 



a^(7i,72 
a^(7i)72 
"^(71,72) = 2a tan {J^ cos 71. 



sm 7i cos 72, 
sin 7i sin 72, 



(7.31) 



The coordinate invariant definition of the normal vector to dR is 



\/9{r, 71,72) efck 



<97i ^72 



(7.32) 
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and the normal unit vector is 

nk= [g'^{r,li,l2)i'ii'j) i^k- (7.33) 

In this case we find 

ni = cos^ sin7i COS72, n2 = cos^ sin7i sin72, ns = cos^ COS71, (7.34) 

and the differential equations for the maximally smooth triad field reduce to 

A0 + cos - sin V^- = ^ [ - sin^(x29i0 - x^^a^) - cos^ cosV^(x39i^ - x^Sg^) 

+ sin sin'?/'(x^(9i'?/' — x'^d^'ip) + cos^sin'?/'(x^926' — x^d^O) + sin6' cos'?/'(x^(92'?/' — x^d^il))\ 

(7.35) 

for 0, 

2 

A6' — sin^ V?/' ■ V0 = — \co?,ip{x^diip — x^d^ip) — sin'?/'(x^92V' — x'^d^ip) 

— sin 9{x'^di(f) — x^d2(p) + cos 9 cos ip{x"^di(f> — x^d^cp) — cos 9 sin ip^x^d-icj) — a;^9i0)] 

(7.36) 

for 6', and 

2 

+ cos9A(j) - sm9 V9 -Vcp = — \ - cosi){x^di9 - x^d^9) 



L2 

— sin 9 sin (f){x^di(j) — x^d-^cp) + sin 'i/j{x'^d29 — x'^d^9) — sin 6' cos '4'{x^d2(p — x'^d^ 



(7.37) 

for ^Z'. The boundary conditions reduce to 

— + cos 6' — = 0, ^ = 0, — + cos 6' — = 0. (7.38) 
on on on on on 

For this type of boundary condition the solution to the differential equation is (j) = 9 = ip = 
const, up to a global rotation and we can set all the constants to zero. Hence, with the 



reference triad field (7.28) we have already found the maximally smooth triad field. The 



averaging process then yields 

''oo PIT p27r 



{E^j{0))= / / /M(2;(r,7i,72))Q/(2;(r,7i,72))^^(r,7i,72)v/^(^^7^^rfrc/7irf72 



'0 Jo Jo 

r'OO /"TT /*2lT 



= / / /M(2(r, 71,72)) A/fi'(r, 71, 72)c/rc/7ic/72 = Saj, 

Jo JO JO 

(7.39) 

since the averaging function is normalized to unity. 

Therefore, the averaged metric is identical to the original metric, 

{9^,m=S^,=g^M, (7.40) 

and we conclude that the smooth three-sphere is not affected by the averaging process. 
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7.5 Averaging the Perturbed Three-Sphere 

In this section we investigate the effect of the averaging process when the surface of the 
three-sphere is subjected to arbitrary linear fluctuations. While we do not present a closed 
result, leaving such to future study, we demonstrate that our formalism is capable of 
dealing with such spaces. The obvious immediate application lies in cosmology, for which 
the perturbed three-sphere forms the basis of the closed FLRW model 



7.5.1 The Metric from Stereographic Projection 

Following the same procedure as for the two-sphere, we perturb the three-sphere slightly 
and apply the stereographic projection. This yields the coordinate transformation 

e^ = (l + r^/(x))— x\ 

where i G {1,2,3}. The perturbation is encoded in the perturbation fuction /(x) and 
?7 > is a small constant. The line element of the perturbed sphere after the projection is 
to first order in r] given by 

ds' = (dep)' + {dil? + m? + {d^pf = (1 + 2^f) (^) V^')' + ^d^"f + ^d^'f) 

(7.42) 

and therefore the metric and its inverse are 
g,, = (1 + 2r//(x)) and g^^ = (1 - 27^/(x)) <5^^'. (7-43) 

7.5.2 The Geodesies and the Connector 

The Christoffel symbols of the metric ( |7.43 ) are 



T), = ^ {-x,6^, - x,5\ + x^6,,) + V (^5^- + - S-£^6,,^ . (7.44) 

For the parametrization of the perturbed geodesic we need three parameters indicating the 
deviation from the unperturbed geodesic in each of the space dimensions. We therefore 
make the ansatz for the geodesic 

zUr) = z\t) + T] 5'^e,uiv\T) — {T). (7.45) 

An appropriate choice for the direction of the w*(r) is perpendicular to the tangent vector 
of the unperturbed geodesic and hence we demand 
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Making use once more of the requirement (7.11) we find the arc length of the perturbed 
geodesic 

s{r) = T + ri [ fir')dr'. (7.47) 



Since the arc length is changed through the perturbation we have to use the general geodesic 

(7.48) 



equation (A.7). With (7.46) and the identity 

{L{zp)f = Aa^ cos-^ ( 



T 

2a 



which holds to first order in we find 



It? 



+ eij/c^ + cos 



-2 



2a J dr 



df 



(0)1 



-6 



jk 



dx^ 



dz^ 
dr 



(0) = 0. 



For non- vanishing initial values dz^/dT{0) ^ this can be further simplified to 



d V 
dr"^ a^ 



cos 



2a 



dz^ 
dr 



(0) 5^=' 



dx'- 



(7.49) 



(7.50) 



These are three independent driven harmonic oscillator equations for the functions f*(r). 
We require the perturbed and unperturbed geodesies and their first derivatives to coincide 
at the origin, 

4(0)=;2X0) and ^(0) = ^(0), (7.51) 



and therefore we impose the initial conditions 



dr 



dv^ 



dr 



v'ifi) = and —(0) = 0. 

dr 



(7.52) 



Now we switch to geodesic coordinates, namely the parameter along the geodesies, r, and 
the angles which describe the direction of the geodesic in the origin, 71 and 72. The latter 
were introduced in (7.30). In terms of these coordinates the first order expansion of the 
connector becomes 



V{t, 0; Cor) = cos-2 (1 - r]f{T, 71, 72) + r//(0, 71, 72)) 1 



+ 1] cos 



-2 



^) (d^T, 71, 72) + I tan (^) v,iT, 71, 72)") 5'^ A, (7-53) 



where the Aj indicate the matrices 





Ai = -1 , A2 = 0, A3 = 1 . (7.54) 
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These are the generators of the rotation group. They form a Lie algebra with the commu- 
tator relations 

[Ai, Aj] = eijk^k- (7.55) 

The jfc-element of Aj is given by 

i^i)jk = ^ijk, (7.56) 



with which we can express the components of the inverse of the connector (7.53) as 



V\{0,T;C^o)=iV-'fiiT,0;Cor 



■ cos ^ ' 



^^)(l-r//(r,7i,72)+r?/(0,7i,72))5^ 



\2aJ ' ' ' ' ' (7.57) 

+ r/cos"2 (^^) (^drVi{T,'ji,'j2) + ^tan (^^j t^i(r, 71, 72)^ S'^ejki- 

Finally, we need the hatted connector for the parallel transport of the triads in the averaging 
process, which is defined as 

i>/(0, r; Cro) =gAO)VMT, 0; Co.)/(0) 

= cos^^ (^) (1 + VfiT, 7i, 72) - ^?/(0, 7i, 72)) |,^_gg^ 

+ 7] cos~^ j (^drVkir, 7i, 72) + ^ tan j 17^(1^, 7i, 72)^ efciii^*'- 

7.5.3 The Maximally Smooth Triad Field 



Since the metric (7.43) is proportional to unity, the best choice of reference field in geodesic 
coordinates is _ 

E^'i = cos' (^) (1 + r7/(r, 71, 72)) (7.59) 



The vector fields (7.20) are 



= (d - m^s...' - ,,^11,) , (7.60) 

For the differential equations that determine the maximally smooth triad field we need the 
value of these vector fields on the surface of constant geodesic distance r =, which is the 
boundary of the area R. With the function 

i^(r,7i,72)= / /(r',7i,72)dr' (7.61) 

the parameter representation of the surface dR can be expressed as 



«^(7i,72) = 2a tan (^) ^(0) - ^(0) + , 7' |^(0). 

^'''^ \2aJ dr^ ' ' cos^ (^) dr^' ' cos2 (^) dr^ ' 

(7.62) 
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The value of on the boundary dR is 



L' {z{r)) 



cos^ 



, sin i-^) 
1, 4oF(r, 71.72)— ^^^^ 



and the square root of the determinant of the metric on dR is 
\/9p {z{r)) = (l + T] 2/(r,7i,72) +r/ -F(r, 71, 72) tan 



(7.63) 



r \ 
, , cos" I — 
2aJ \2aJ 



(7.64) 



Thus, the values of the non vanishing vector fields ( |7.20[ ) on the boundary are 

(Q'^Y {z{r)) = -iQ^'Y {z{t)) = -iQ'Y (^W) = iQ^'Y {z{r)) 

^ ( iMi^ _ 1 _ ^m^Jl^M _ 2/(r) sin(^) \ dz^ 

l^acos3(^) ''2a2cos^(^) ^ cos^ (^) ^ a cos3 (^) J c^r ^ ^7.65) 



11/ \ 1 9f 



!;=z(r) 



and 



(gi2)2 (z(r)) = -iQ^'Y {z{r)) = -{Q'^Y {<r)) = (Q^'Y (^W) 
( i si^(l^) 1 F(r) sin^ (^) 2/(r) sin (I,) \ 

V acos3(i,) '2a2cos4(^) 'a^ cos^ (^) a cos3(^);ar 

11/ (9^3 \ 1 df 

-^^^ -^^(0^(0)) ^^^^ 



(0) 



(7.66) 



and 

'1 si^(^) F{r) 1 F{r) sin^i-J 2f{r) sin {j,) \ dz' 

cos^ (^) «^ cos^ (1^) « cos3 {^J ) dr^^ (7.67) 



yacos3(^) '2a2cos4(^) ' cos^ (|-) a cos^-J J Or {7 .67] 

-'i^ ("'"C'"') -^<^)|^<°') -'^4sl^L„. 

The normal unit vector to the surface of constant geodesic distance dR has the component; 

/ 2 f \ r if^\ F f \ 'iFcos^i^)\ . 

ni = cos — + r]/ cos U~ + ^— cot — - rj . , . sm7i COS72 

y \2aJ \2aJ a \2aJ a sm (^) ) ^^ ^g, 

1 f r \ f dF dF sin 72 2 s • - A 

+r]— cot — -— COS71 cos 72 - 7^ : h v cos 71 - f''sm7i sm72 , 

2a \2a/ \c'7i (772 sm7i / 
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and 



sm 7i sm 72 



1 f r \ f dF dF cos 72 1 o \ 

+77— cot — - — cos 71 sm 72 + : V cos 71 + sm 71 cos 72 , 

2a \2a/ \0'yi 072 sm7i / 



and 



1 -? f r \ , f r \ F fr\ 3Fcos^(:^) , 



1 ( "f \ ( dF 



d'ji 

For simplicity we introduce the three scalar fields 



dF 

sin 7i — v"^ sin 71 cos 'y2 + v^ sin 71 sin 72 I . 



(7.69) 



(7.70) 



A,{x\x',x') = p.,x^-7^^^. (7.71) 



We can then express the coupled differential equations for the maximally smooth triad 
field as 



A0 + cos^ AV^-sin^ • = - sm9{A2di9-Aid29) - cos9 cos^iAsdiO-AidsO) 

+ sm9sm'ilj{A3diilj-Aid3ilj) + cos6'sinV'(^3a26'-^2536') 

+ sin 9 cos ip{A^d2ip — A2d3^) 
A9 + sin 6* Vip-V(f) = cosijiAsdiip-Aidsip) - simp{A3d2iJ-A2d3ij) 

+ sin 9{A2di(l)-Aid2(t)) + cos 9 cos i/j{A2di(j)-Aid3(l)) 

— cos9sm^p{A^d2(t)—A2ds(t)), 

A^ + cos^ A0-sin^ ■ V0 = - cos ^(^3^1^- Ai^a^) - sin ^ sin ^'(^3910- Ai^g^) 

+ sin t/^ ( ^3^26' - A2936') - sin 6* cos (^3920 - ^2530) . 

(7.72) 

For the boundary conditions we introduce furthermore the quantities 



^i(T, 71,72) = ^ijk 



df 



dzp 

.97 



(0). 



(7.73) 
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Then the boundary conditions become 

r)rh dijj COS 9 f 1 dF 



dn dn cos^ (^) V 20^972 

sin 6 cos i}) f 1 dF 1 dF 

+ ^ ^^^co^72 + ;7^7^cot7iSin72 + /i2 r,7i,72) 

sin^sinV- f I dF 1 dF , , A 

+ '^^^ 1,2^ + 2^ + j ' 

89 sm9 f 1 dF 

dn = ^^^^ (,2^9^^°^^^ - 

COS 7/^ f \ dF 1 
+ — 7^p~ + cot 7i cos 72 + /ii(r, 71,72) ) , 
cos2 (^j \2a^ d7i 2a^ 



90 1 f i 9F 



(7.74) 



These coupled differential equations implicitly define the angles 0, 9 and and therefore 
the maximally smooth triad. Although the solution of these equations is very involved 
and therefore cannot present any specific example, the analysis shows that no principle 
difiiculties arise when the method is applied to more than two space dimensions. 



CHAPTER 8 



THE GAUSSIAN SHAPED PERTURBATION 



In this chapter we analyse numerically the results of chapter [6| in which we derived the 
equations governing the average of a two-sphere perturbed by a function /(x^,a;^), for the 
specific example of a Gaussian perturbation. The differential equation for the maximally 
smooth dyad field cannot be solved analytically and we make use of the numerical toolkit 
Gascoigne. After briefly introducing the method of finite element solutions we apply it 
to the problem. This chapter is based on work performed in collaboration with Thomas 
Richter. 



8.1 Finite Element Solutions 

Consider Poisson's equation, 

-A</) = (8.1) 
in the domain Q with a Neumann type boundary condition on dQ, 

^'^ V^ir,7)+Vi^,^rT^dJ r fis\^)ds') ^ g^, (8.2) 



which we will solve numerically. A detailed discussion of the use of the finite element 
approximation is given in [21], while we present a brief review of the essential points. 



The "weak formulation" of the boundary value problem is obtained by the multiplying 



(8.1 ) with a test function ^ and integrating over the domain f2 to give 

- I A0 ■ ^rfx = 0. (8.3) 

Integration by parts gives 



V0-Vec/a;= / ^-^f^o. (8.4) 



Q. Jdn 
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The weak formulation of the partial differential equation is given by G (where V is an 
appropriate function space) with 



V0 ■ Vidx = Qn- ^do, G V. (8.5) 
n Jan 

If is solution of this weak equation and is sufficiently regular - the Laplacian of has to 
exist - then we can reverse the integration by parts and see that is also the solution to 

- f A(j)-^dx = 0, with f ^■^do= [ Qn- ^do, G V. (8.6) 



One can then show that is a solution to equation (8.1) with the boundary condition (8.2) 
for every point of the domain VL. 



In the finite element method one solves the weak formulation (8.5) by introducing discrete 
subspaces Vh C V and searching for solutions Uh G Vh. This requires the discretization of 
the spatial domain Q with a regular triangulation Qh, a mesh of "finite elements", and the 
introduction of local basis functions of Vh on this mesh. The finite element triangulation 
Qh of the domain f2 is a set of open elements K eVL^ with the following properties: 

1. Every element K E VLh is a, quadrilateral and all interior angles a oi K are bound 
uniformly a < Uq < ti for all elements of VL^. 

2. The ratio of the diameter of the element dx and the radius of the largest inscribed 
circle pk is uniformly bounded in the grid dx/ Pk < Po- 

3. The closure of the elements covers the domain: 

VL = [jK. 

4. For every two elements and in VLh with ^ the intersection fl K"^ is 
empty and the intersection of the closures n K"^ is either empty, a common edge 
or a common point. 

On this finite element mesh Vth we define the space of piece-wise polynomial functions 

Vh = {v eV : v\K e span{l, x, y, xy, x'^,..., x^'y''}}. (8.7) 

Using the grid points Xj, the finite element basis {^j} of this function space is usually 
chosen as 

U^j) = (8.8) 

Thus, every basis function differs from zero only in the elements adjacent to the node 
Xi. In terms of the basis the function 0/^ G can be expressed as 

N 



J2 with 0^ G E. (8.9) 



i=l 



8.2 Approximation of curved boundaries 
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Figure 8.1: Approximation of a curved boundary with a linear and iso parametric (here 
quadratic) transformation to the computational element. 



The numerical solution of (8.5) in the finite element discretization is then the solution 

(ph e Vh to 

N „ N „ 

V 0' / ■ V^rfx = V 0* / . i^dx, Vj = 1, . . . , iV. (8.10) 

i=i i=i 
Since the basis functions have a small support, the integrals cover only small areas of the 
domain and every function couples only to a small number of other basis functions ^j. 
For a two dimensional computation with piece-wise quadratic finite elements every basis 
function couples usually to no more than 25 neighboring basis functions. Equation (8.10) 
gives a system of N linear equations with N unknowns. The resulting matrix is very sparse; 
every row has of the order of 25 non-zero entries. This system can be solved efficiently 
with a multi-grid algorithm 



8.2 Approximation of curved boundaries 

The discretized domain is not identical to the original continuous domain and the ap- 
proximation by quadrilaterals leads to a significant error when approximating the curved 
boundary. This must be carefully considered. To gain a better approximation of the 
boundary we use so called iso-parametric finite elements. Instead of defining the finite 



element space directly on the computational grid VL^^ as in (8.7), we define a reference 
function space on the unit square K = [0, 1]^ 

QJ- = {v e span{l, x, y, xy, a;^ . . . , x^'y''}}, (8.11) 

and use a nonlinear mapping onto the computational elements Tk : K ^ K E Qh- This 
mapping Tk is an element of the same reference space Q". The finite element space Vh is 
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Figure 8.2: The Gaussian shaped perturbation function. 



then given by 



VH = {veV:^v^ Q\v\K = Tk{v)]. (8.12) 
To improve the definition of the mapping on an element K at the curved boundary, 



we add further mesh points. Figure (8.1) shows an example in which this iso-parametric 



approach results in a much finer approximation of the boundary. 

This finite element method with iso-parametric elements is implemented in the software 
library Gascoigne 3D pQ developed by M. Braack, D. Meidner, R. Becker, T. Richter, B. 
Vexler at the universities of Heidelberg, Linz, Kiel and Pau, which we will employ in the 
remainder of this chapter. 



8.3 The Perturbed Two-Sphere Revisited 

We now apply the formalism outlined above to the analysis of the perturbed sphere pre- 
sented in chapter [6j As a concrete example, we choose a sphere of radius a = 10 and a 
perturbation function of Gaussian form. In general, the Gaussian function is of the form 



f{x^,x'^) = D exp 



a' 



(( 



,1\2 



X 



(8.13) 



To keep the perturbation small, we choose the amplitude D = 0.1 and standard deviation 



(7 = 1. The resulting function is shown in Figure 8.2 



We would expect a reasonable averaging process, when applied to the perturbed two-sphere, 
to significantly diminish and broaden the perturbation function., which should furthermore 
tend to disappear for increasingly large averaging radii. At the same time the radius of 
the background sphere might be slightly altered. 
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Figure 8.3: The Gaussian shaped perturbation function in the coordinate systems of the 
reference points (xqjXq) = (1,2) and (xqjXq) = (20,20). 



To take the average we must first transform the perturbation to the coordinate system of 



the reference point {xI,Xq) using (B.6). Doing so, we see that the transformed Gaussian 
function is 



f{x^, x^) = D exp 



1 



Dexp — 



16a^ {{xl)'' + (xg)^ + 2{xlx^ + xlx^) + {x^f + {x^f) 
(t2 16a4 - 8a2(xjP + xgP) + {{xlY + {xlY){{x^Y + (P)2) 



(8.14) 



From this we can see that the transformed function keeps the same amphtude but grows 



broader for increasing {x}^,x'^). In Figure 8.3 we show the transformed function for the 



reference points {x\,Xq) = (1,2) and (xo,Xq) = (20,20). This clearly shows this broaden- 
ing at distant reference points, which originates from the larger angle of rotation at such 
distances, and from the subsequent stereographic projection. 

In terms of the geodesic coordinates the function /(r, 7) is 

1 4a2((xJ)2 + (a;o)2+4atan(^)(xoCOS7+Xosin7)+4a2tan2(^)) 



/(r,7) = L)exp| 



(T^ 4a2 — 4atan(^) (xj COS7+X0 sin 7) + {{xiy + {xiy) tan' 



2/' 

2a' 



and the function /;,(r, 7), defined in (6.12) is 
/i(r,7) = 



(4a2 + (xj)2 + (xiy) (xj sin - xl cos 7) /(r, 7) 



,2 • 



0-2 cos2(^) (4a2-4atan(^)(x|^ cos7+x§ sin7) + ((x,^)2 + (xg)2) tan2(^))' 
This function defines the longitudinal deviation f('r, 7) of the perturbed geodesic via the 



differential equation (6.11). With f(T, 7) we can compute the connector (6.57) and the 



boundary of the geodesic region, dR, which in parameter representation is given by (6.42). 
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The differential equation for the scalar field x^), which defines the maximally smooth 



dyad field, is the Laplace equation (6.39) with the boundary condition (6.47) on dR. The 



numerical solution to this equation determined with Gascoigne for the reference point 
(xj, Xq) = (1,1) is shown in Figure 8.4 for the averaging radius r = 3 and in Figure 8.5 for 
the averaging radius r = 5. Since the perturbation function enters the differential equation 



only through the boundary condition. 



[X 



, X ) is essentially zero in the region far away 



from the perturbation and also just across the perturbation. Furthermore, the maximum 
value is larger for the smaller averaging radius. 



uOOO 

0.0198 




■3.55e-06 



■0.00659 



■D.D132 



■0,0198 



Figure 8.4: The result for (f){x^,x^) for the references point {xI,Xq) = (1, 1) and r = 3 
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Figure 8.5: The result for 0(x^,x^) for the references point 



^Xq, XqJ 



[1, 1) and r = 5. 



With the solution for </)(x , x ) for all reference points in a domain around the pertur- 
bation function /(x^. 



X 



we can construct the averaged metric from (6.96). The initial 
perturbation and the averaged perturbation function for the averaging radius r = 3 are 



shown in Figure 8.6 It is clear that the perturbation function is not much altered by 
the averaging process and that the effect of averaging is much smaller than expected. A 
reasonable averaging process should diminish the perturbation significantly and therefore 
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Figure 8.6: The perturbation before and after averaging over a region of radius r=3. 



we must conclude that, in its current form, the averaging process is of hmited use. We 
will therefore elaborate on potential problems in the process, and discuss possible improve- 
ments, briefly here and in greater depth in the next chapter. 

The problem might arise from the extended spread of the bundle of parallel transported 
dyads at each reference point, causing the averaged dyads to be shortened when compared 
to the maximally smooth dyad. A first ansatz to improve the averaging process is therefore 
to renormalize the averaged dyad at each reference point witht the variance ((0 — (0))^) = 
(0^) — (0)^, but this can be shown to be untenable. For the perturbed two-sphere we find 

2tt pr 



'■)-{^f = I-' / / /2(r)02(r,7)cos^(^)rf7rfr 

2n r ^2 (8-15) 



The solution for an averaging radius r = 3 and for all reference points is shown in Figure |8]7j 
Due to the symmetry of the perturbation function at the reference point {xq^xD = (0,0), 
the boundary condition vanishes and therefore 0(a;^,x^) = 0. The variance thus also van- 



ishes for this reference point which can be seen from the central dip in Figure 8.7 The 



variance can therefore not be used to change the length of the averaged dyad at this point. 

A more promising approach that we are currently working on is the renormalization of the 
averaged dyads with the renormalization factor 



Or 



1 + -"^^^^) > (8-16) 





AR 




\AR 


max 



where /z is a constant that must be defined according to the problem under consideration, R 
is the Ricci scalar and | Ai?|max is the maximum value of AR in the averaging domain. This 
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Figure 8.7: Variance for the averaging radius r = 3. 



renormalization factor appears only in the integrals in (6.85) and not in the normalization 



(6.77). For the perturbed two-sphere the factor becomes 



aren = l-w[ "A/ + {xd.{Af) + ydy{Af)) + —,A{Af) 



(8.17) 



and can completely be included in the scalar field i?(r, 7), defined in (6.84), to give 

\AR\ 



Sre„(r,7) = 5(r,7) 



|Ai?|, 



This approach must alter the length of the averaged dyads and results are expected in the 
near future. 



In this chapter we have applied the averaging process to the metric of a perturbed two- 
sphere for the specific example of a Gaussian perturbation. Solving the differential equation 
for the maximally smooth dyad field numerically with the toolkit Gascoigne has led to an 
averaged metric with a perturbation that is essentially unaltered by the averaging process. 
We elaborated on the reason for the small effect of the averaging process and presented an 
ansatz which might improve the process. Another approach which aims at the core of the 
problem will be presented in the next chapter. 



CHAPTER 9 



THE CURVATURE LAGRANGIAN 



In this chapter we argue that the maximally smooth Vielbein field employed in the av- 
eraging process thus far is too close to the specific Vielbein field that is invariant under 
the process. We suggest a different Vielbein field which will improve the averaging ef- 
fect by taking the local curvature of the manifold into account. This Vielbein field will be 
determined from a modification of the Lagrangian which determines the invariant Vielbein. 

9.1 Drawbacks of the Averaging Process 

In the previous section we applied the averaging process defined in chapter |4] to a two- 
sphere perturbed with a Gaussian perturbation. A reasonable averaging process would 
be expected to diminish this function and tend to remove it for larger averaging radii. 
Instead, our procedure left the perturbation basically invariant and therefore the process 
needs to be modified to overcome this problem. We presented in the previous chapter a 
renormalization factor that might address the shortening of the Vielbeins under parallel 
transport. In this chapter we instead consider a redefinition of the Lagrangian used to 
determine the Vielbeins in the first place. 

There is one particular Vielbein field that is invariant under the averaging process. This 
Vielbein field can be constructed by parallel transporting the reference Vielbein from the 
origin £^"^(0) to each point of the geodesic region R and therefore we will refer to it as 
the geodetic induced parallel field. Clearly, this Vielbein field is only determined up to 
global rotations. The averaging process parallel transports all Vielbeins back to the origin 
and the averaging returns the initial Vielbein E°'i[{)). The Lagrangian which defines the 
geodetic induced parallel field is 

Anv = {£tE\) {£tE\) Sabtn-', (9.1) 

where £t denotes the Lie derivative with respect to the tangent vector to the geodesic 
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bundle through the origin z'^. 

The Lagrangian we chose to determine the Vielbeins was that for the maximally smooth 



Vielbeins (4.43) 



Cms = {D,E\{x)) {D,E\{x)) g'^'' {x)gP^ {x)r^ab. (9.2) 
This Lagrangian leads to a differential equation whose solution (pMS is very similar to the 



solution 0inv for the differential equation defined by (9.1). The average of this Vielbein 
field is therefore only mildly altered. 

Knowing the Lagrangian that defines the invariant scalar field (f)\^^{x^ ^ x'^) allows us to 
modify this Lagrangian in order to find a scalar field that is definitely not invariant un- 
der the averaging process. The modification is chosen to be proportional to the Ricci 
scalar R. Since the modifier should be dimensionless and the dimension of the Ricci scalar 
is (length) we must balance this with a factor of the dimension (length)^. A reasonable 
choice is the length of the geodesic s^. We therefore define the curvature Lagrangian 

Cr = {£tE%) {£tE\) Sab {t'-t'' + g^-Rs^) . (9.3) 

For the fiat case with i? = the modification vanishes and the solution merges into the 
invariant scalar field. The concept of this Lagrangian is to construct the geodetic induced 
parallel field and to then rotate the Vielbeins according to the curvature at their location. 
The rotated Vielbeins are then parallel transported to a common point and averaged over 
to return the averaged Vielbein and therefore the averaged metric as before. 

To further explain our approach, we return to consider the perturbed two-sphere, starting 
with a detailed summary of the important quantities and equations. 

9.2 Summary of the Perturbed Two-Sphere 

In this section we review the important quantities and equations for the perturbed two- 
sphere, which were derived in chapter [6j We also compute the tangent vector that appears 
in the curvature Lagrangian in terms of the coordinates of the stereographic projection 
plane. 

The metric for the perturbed two-sphere after stereographic projection is 



(^7p),, = (l + 27^/(a;\a;^))(^-J 5., and [gpf = {I -2r^f{x\x^)) \^-^ (9.4) 
where I? = Aa? + {x^Y + (x^)^. The square root of the determinant is then 

^=(l + 2r//)(|]' (9.5) 



9.2 Summary of the Perturbed Two-Sphere 
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and a suitable choice for the reference dyad field is 



The geodesic bundle through the origin is 

cos 7 



z^(s, 7) = 2a tan ( — ) C0S7 — 
\2a/ 

z'^{s, 7) = 2a tan ^ 



V 



cos7-r/- 
2a/ cos"^ 



cos2(|^) Jo 
COS 7 



f{s,vds +1] sin7, 



cos 



2^^ -'O 



f{s,-f)ds +r/— — — sin7, 



cos^ 



where s is the arc length of the geodesies, 

s = T + rj / /(r', 7)(ir' and 



r = s-T] / f{s\-i)ds\ 



and f (r, 7) is the solution to the differential equation 

d^v{T) v{t) h{T) 



dr^ 



The function h{T,'j) is 
Of 



hir) 



{x\x2)=(z^{t),z2{t)) 



87^ 



a2 cos2(^)' 

5/ 



(a;i,z2) = (2i(-r),^2(r)) 



The tangent vectors to the geodesies (9.7) are defined as 
dz^ 



cos 7 



ds cos^ 



1 sin(2^) cos 7 

cos 7 / j[s , 7)as — r]- 



. 2a' 



a COS' 



COS^ 



■-2a' 



+ r] 



t>(s,7) sin( 



a cos 



2a > 



■ 2a' 



sin 7 + 



sin 7 (if 



(is 



cos 



2fJ_ 
■ 2a 



dz^ 



sin7 



ds cos^(^^ 



-V 



1 sin(2^) . /" sin7 
sin7 / f{s,vds -7]- 



a COS' 



■■2a' 







COS 



■ 2a' 



- V- 



ti(s,7) sin( 



2a' 



a COS' 



■■2a' 



COS 7 - ?7 — ^ 



cos 7 (if 



cos^l^^j ds 



In terms of geodesic coordinates s and 7 we thus find that 

4a2 



cos^ ' — 



2a' 



COS^ 



and therefore that the metric is 

i9p)ij 



(^l + 2r7/(s,7) + 2r/atan(^)^ 



f{s',l)ds' ] cos ( — j 6ij. 



(9.6) 



(9.7) 



(9.^ 



(9.9) 



(9.10) 



(9.11) 



(9.12) 



(9.13) 
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For constant arc length s = r, where r is the averaging radius, the geodesies (9.7) define 
the boundary dR of the geodesic region R. The square root of the determinant at the 
boundary dR is then 



n/^Im " ("^ ^ '^Vf{r, 7) + tan j ^ f{s', -f)ds'^ cos^ 



2a 



(9.14) 



and the reference dyad field becomes 
E 



IdR 



1 / r 
1 +rif{r,'y) + r]-tan — 
a \2a 



fis',l)ds' cos 



2a 



St- 



(9.15) 



The normal vector to the boundary dR is 



^^1(7) = (1 +^/)cos^ (^) cos7 + ?7^tan (^^j cos^ (^^) COS7 f f{s',-f)ds' 

J 



1 COS (2^) . d 

2asmQ^ rf7io 



/(.',7M.' + ,%^^sm7, 



2a sin(^) 



^2(7) = (1 +^/) cos^ (^) 

1 COS (|-) d 
2asm(^) djJo 



1 ( r 
sin 7 + 77- tan — 
a \2a 



cos 



(9.16) 



/ N . / "^(S' 7) cos (2^) 
/(s', 7)rfs' - ^ \ . ),{ cos 7, 



and so 



a/? 



1. 



(9.17) 



For the coordinate transformation from the geodesic coordinates to the coordinates of the 
stereographic projection plane, we need to invert (9.7). With 

F{s,^)= r f{s',^)ds' 



we find 



s = 2a arctan 



V/(X1)2 + (x2)2 



7 = arctan I ~f I + ^ 



X 



2a 



+ r]F{x^,x^) 



(9.18) 



(9.19) 



The tangent vectors (9.11) to the geodesies in terms of the coordinates x^ and are 



9.3 The Geodetic Induced Parallel Field 



119 



therefore 



X 



1\2 + (^2)2 



[X 



2\2 



x'^v{x^,x'^) 



16a4 + 



+ ^ 3^ — - + v 



1 2 

a; a; 



(9f 



16a4 + (x2)2) + ^15^4 ^3,1)2 + ^^,2)2) Q^2 ^ 



(X 



2\2 



dv 



t\x\x^) 



X 



-V- 



Furthermore, we find 



(9.20) 



16a4 ((xi)2 + (x2)2) dx^ ^ 16a4 {x^x^ + (x^)^) dx^ ' 



'rfs2 



= Tj^^ I X + X 



+ r] 



dx^ 



dv 
9x2 



16a4((xi)2 + (x 



2^2^ 



^^^2 +2x^x^ 



(9x1) 



9x19x2 



+ X 



2\2 



(9x2)2 



(9.21) 



We will need the tangent vectors (9.20) for the Lagrangian (9.1) and equation (9.21) for 



the proof that this Lagrangian defines the geodetic induced parallel field in the case of the 
perturbed two-sphere. This will be shown in the next section. 

9.3 The Geodetic Induced Parallel Field 

For the perturbed two-sphere we can compute the scalar field 0inv(a;^,x^) that defines the 



geodetic induced parallel field with the connector, which was defined in (6.56) 



V^^(r,0;Co.)=cos-2(^)(^ 



(l-r^/(r,7)+r//(0,7))5f 



+ ?7 fw'(r,7) + ^tan ^^(^,7) j e^^i 



(9.22) 



From this we find the hatted connector 
y,^(r,0;Co.) = (7.fe(r,7)V^(r,0;Co.)(7'^(0,7) 



cos 



{- 



(1 + r//(r, 7) - r//(0, l)) 5l + [ v'{t, 7) + - tan - v{r, 7) e,^ , 



2a 



(9.23) 



where the prime denotes the partial derivative with respect to r. This connector can be 
used to parallel transport the reference dyad from the origin, 



E^(0) = (l + r//(0))51^ 



(9.24) 
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to all points of the manifold. The local rotation of the resulting dyad field as compared to 
the reference dyad field then defines the scalar field ipinvix^yx"^), 



and the result for the scalar field is 



1 2 

X TT-r + X 



In geodesic coordinates this result is 

dv 



dv 



(9.25) 



(9.26) 



2a/ ds 



V5inv(s,7) = - -tan 
a 

The Lagrangian that describes the geodetic induced parallel field is 

-^^inv = {£tE"'i) i^tE^j) Sabt^'t^- 



(9.27) 



(9.28) 



With the definition of the Lie derivative (A. 13) this Lagrangian leads to the differential 
equation 



on the geodesic region R with the boundary condition 
on dR, where 



and 



= Pf = tHhdfE\Ef ^{Dkt')t\ 



(9.29) 

(9.30) 

(9.31) 
(9.32) 



With ( |9.6[ ) and ( |9.20[ ) we can see that 

( I dv dv 



Q + P = r]—+7] 



X 



+ 



Aa'^ \ dx^ dx 



+ ?7 



16a4((xi)2 + (x 



X 



1x2 d'^'V 1 2 



{dx 



l\2 



dx^dx'^ 



+ [X 



2\2 



d^v 
{dx'^y 



(9.33) 



An explicit calculation using (9.26) shows that 



t^djLfi, 



iP + Q) 



(9.34) 



and thus proves that 



(x^,x'^] 



^mvix\x ) is actually the solution to the differential 



equation (9.29) with the boundary condition (9.30), and hence is the geodetic induced 
parallel field needed in the curvature Lagrangian, which we now turn to. 



9-4 The Curvature Lagrangian 
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9.4 The Curvature Lagrangian 

For the perturbed two-sphere the curvature Lagrangian is 

Cr = {£tE\) {£tE'j) 6ab {tV + g'^Rs^) . (9.35) 
This Lagrangian yields the differential equation 

A (tV(l + 2Rs^){dj(t))) = -Di {P' + Q' + P}i + gy (9.36) 
with the boundary condition 

n,t¥(l + 2Rs^) {dj(f)) = -n, {P' + + + , (9.37) 

where 

P}, = Pnf = Rs''g^hdf&i&,{Dj,t')t\ (9.38) 

and 

= Q^f = Rs^&''edf{Di&j)&kth\ (9.39) 

If we decompose the solution into the scalar field 0inv and the rest (pR, 

(t){x\x') = (t);Ux\ x^) + Mx\x'), (9.40) 



some of the terms in the differential equation cancel due to (9.29) and we can add the 
remaining terms to the right-hand side. We then find 

A (tV(l + 2Rs^){dj(t)R)) = - A {P'r + Q'r + t¥2i?s2(9,0inv)) (9.41) 

with the boundary condition 

t={l + 2Rs^) {d,<PR) = -n, [PI, + + tH^2Rs\d,cp,^,)) . (9.42) 

The solutions to these equations then give us the scalar field (j)R{x^,x^) that takes the 
curvature into account and determines our new Vielbein. These equations are in principle 
solvable numerically and this research is ongoing, to be completed in the near future. 

9.5 Discussion and Conclusions 

In this chapter we presented the Lagrangian which defines the Vielbein field which is invari- 
ant under the averaging process and proved this invariance for the case of the perturbed 
two-sphere. Based on this Lagrangian we suggested a Lagrangian defining a Vielbein which 
differs from the invariant Vielbein provided that the Ricci scalar does not vanish. Since 
the modification term in the new Lagrangian includes the Ricci scalar, we expect the asso- 
ciated Vielbein to be a rotation of the invariant Vielbein about an angle which is larger in 
regions of higher curvature. The spread of the Vielbeins which are parallel transported to 
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a common point in the averaging process is therefore also larger in these regions and thus 
the average is sensitive to the local curvature. For dimensional reasons the modification 
term in the Lagrangian also includes a factor of s^. This factor vanishes at the origin and 
therefore we do not expect any issues with discontinuities. For more distant points the 
factor grows larger but its influence is counterbalanced by the averaging function in the 
process, which is required to weigh the contribution of Vielbeins from more distant points 
at a smaller rate. Future research will show if the curvature Lagrangian defines the Viel- 
bein fields that it was designed for and if the averaging process based on the new Vielbein 
field leads to the expected results. 



CHAPTER 10 



CONCLUSIONS 



In this thesis we have addressed two aspects of the averaging problem in general relativity. 
The first is the cosmological backreaction arising from the non-commutativity of volume 
averaging with the construction of the Einstein tensor and the second is the problem of 
constructing a generally covariant average metric. Both the average of a metric and the 
average of the Einstein tensor are of great importance in a variety of situations and cos- 
mology forms a useful testbed for both. 

It has been suggested in the literature that cosmological backreaction can account for the 
observed acceleration of the universe but there have been few quantitative studies. We 
review the standard cosmological model and some previous work in this area in chapter [2] 
In chapter [3] we considered a linearly-perturbed Friedmann-Lemaitre-Robertson- Walker 
universe. Previous work has tended to be in synchronous gauge which contains unphysical 
gauge modes. We chose instead to work in Newtonian gauge. Employing the Buchert 
formalism and the cmbeasy Boltzmann code we calculated the backreaction terms numer- 
ically for the WMAP concordance cosmology and for a toy Einstein de Sitter model. In 
both cases we showed that the backreaction of linear perturbations is too small to account 
for the observed acceleration of the universe and that the corrections act as a dark matter 
and not a dark energy. This result agrees with the naive expectations. Using the Halofit 
code we extended our analysis to mildly nonlinear scales. The backreaction from quasi- 
linear scales remains of the order of 10~^ for both the ACDM and EdS models and the 
corrections continue to act as a dark matter. However, although the backreaction effect is 
small it is nevertheless present and is expected to be significantly larger on smaller scales 
that we did not consider. On such scales vector and tensor perturbations become signifi- 
cant and for an unambiguous answer to whether backreaction effects significantly influence 
the evolution of the universe the Einstein equation will have to be averaged in a covariant 
manner. While our result is small it does not rule out a significant backreaction. It is also 
a prediction of the backreaction from viable cosmological models. 
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The other side of the averaging problem concerns the metric. In chapter |4] we presented 
a generally covariant averaging process which can be used to smooth metrics within the 
framework of general relativity. This process involves the decomposition of the metric into 
a specific choice of tetrads, namely the maximally smooth tetrad field. Within the geodesic 
region to be averaged across this tetrad field is determined by minimising an action whose 
Lagrangian describes the variation of the covariant derivative of the tetrads. The maxi- 
mally smooth tetrads are then parallel transported to a common point using a connector. 
This connector is a general relativistic reformulation of the Wegner- Wilson operator from 
QCD. At the common point the average can be unambiguously taken. The averaged metric 
is then obtained from the averaged tetrads. 

As a first example, we applied the process in chapter [5] to the two-dimensional spaces of 
constant curvature, the two-sphere, the two-plane and the two-hyperboloid. These spaces 
provide a convenient and simple tool to visualize the averaging process. Moreover, any 
reasonable averaging process should leave a space of constant curvature invariant and so 
they also act as a consistency test. Finally, the three-dimensional spaces of constant cur- 
vature form the basis of the FLRW models and so an understanding of averaging on such 
surfaces has an immediate cosmological application. Considering first the two-sphere we 
found that all constituents of the process, especailly the geodesic bundle through the ori- 
gin, can be expressed in an elegant way after the applications of stereographic projection. 
It was shown that the averaging indeed leaves the two-sphere unaffected as we required. 
The averaging process was then applied to the hyperbolic plane and the flat plane. Both 
were also invariant. 

In chapter [6] we extended the analysis to a slightly perturbed sphere. This serves as the 
first vital step towards a covariant average in three space dimensions as would be required 
to reconstruct the FLRW metric. The Lagrangian formulation was shown to yield a partial 
differential equation of Neumann type for the maximally smooth dyad field. This equation 
is not analytically tractable but can be evaluated numerically. 

In chapter [7] we showed that the process can in principle be applied to space of more than 
two space dimensions by considering the three-sphere. This space of constant curvature 
corresponds to a closed FLRW universe along a hypersurface of constant time. We showed 
that the three-sphere is also left invariant by the averaging procedure as should be expected 
but that the differential equations for the maximally smoothed triad gains a complicated 
structure that will be difficult to solve. In more convenient coordinates it will be solveable 
numerically. 

In chapter [8] we turned to solving the equations for the perturbed two-sphere. We speci- 
fied the perturbation function of the two-sphere to be a Gaussian function and solved the 
partial differential equation with the aid of the numerical toolkit Gascoigne. This led to 
the conclusion that the averaging process, as originally formulated, does not smooth the 
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perturbation function in the desired way. 

The reason for this was analysed in chapter |9| where we concluded that the Lagrangian 
generating the maximally smooth dyad field is similar to the Lagrangian that defines the 
dyad field which is completely invariant under the averaging process. This problem was 
overcome by adding a curvature-dependent term to the Lagrangian. By including the Ricci 
scalar at a general point in the Lagrangian we take into account the curvature of the man- 
ifold and directly link its fluctuations to the average. 

This averaging process can now be applied to different perturbation functions to study 
their interaction with each other and with the background sphere. Being applied to a 
three-dimensional sphere the outcome of the averaging would be of cosmological relevance 
for a closed FLRW model. Considering instead the fiat three-plane with perturbations 
would correspond to a fiat FLRW model. The immediate prospect is then of applying the 
formalism developed in chapter [7| along with the experience gained in chapters [8] and [9] to 
linearly-perturbed FLRW models and so test to some degree the cosmological principle. A 
more general study would require us to consider fluctuations of arbitrary size. Concerning 
the extension to four dimensions would involve an analysis of the boundary conditions on 
the congruence of light-like geodesies. 

The averaging problem in general relativity is one of the greatest unresolved problems in 
classical physics. In this thesis we have considered two aspects of it. We have calculated 
the corrections to the averaged Einstein equation compared to that for a smooth model 
that arise inevitably from the current concordance model and from a toy competitor model 
without dark energy. The impacts from the scales our approach considered are small and 
resemble a dark matter but indicate that backreaction remains a physical effect. We 
have then proposed a generally covariant averaging process for metrics based on parallel 
transport and much simpler than other existing procedures. We have applied this procedure 
to test metrics and analysed the perturbed two- and three-spheres and demonstrated its 
applicability. 



APPENDIX A I 

L FUNDAMENTAL CONCEPTS OF GENERAL RELATIVITY 



Throughout the thesis we use units in which c—1. 

In the general theory of relativity gravity is described as an effect of the curvature of 
spacetime, and the geometry of spacetime in turn is related to the energy and momentum 
content. 

Spacetime {A4, g) is described as a four-dimensional real smooth manifold with a symmetric 
nondegenerate metric g of Lorentzian signature (+, —,—,—). In a coordinate basis we can 
expand the metric in terms of its covariant components gfj_i, as 

Qf^u dxi" ^dx". (A.l) 

In abstract index notation the metric tensor is expressed as the line element 

ds^ = g^^dx^dx^ (A.2) 

The contravariant components of the metric are defined by 

g'^'g.p = S',- (A.3) 

If we change coordinate systems the contravariant components and the covariant com- 
ponents An of a vector A will transform according to 

^"-t?^' and (A.4) 

A geodesic Cx^x^ is the path of extremal length between two spacetime points Xq and Xi. 
If we choose an explicit parametrization for the geodesic 

= = ^'{r), ro<T< n, ^^(to) = x^, zf^in) = x>l}. (A.5) 
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and set the variation of the arc length s(C) 

s{C) = I dr^g,Mr))^{r)^{r) (A.6) 

with respect to r equal to zero, we find the general geodesic equation 

d I 1 , , , , dz^ , , \ 11 dq^u , , , , dz^ , , dz^ , , / . „v 

Tr \^M)^ir) j - = 0, (A.7) 

provided the variation vanishes at the end points and X\. Here, the arc length s as a 
function of r is given by 

sir) = £ dT'^g,Mr'))^{r')^{r'). (A.8) 

For ds"^ > wc can use the arc length s to parametrize the curve and therefore set r = s. 
Then the geodesic equation reduces to 

d^ d_z^dj^^ (A.9) 
ds^ ds ds 

with the Christoffel symbols 

The covariant derivative -D^ of a contravariant and a covariant vector are defined by 

D^A^^d^Af + r'I^.A'' and D^A, ^ d^A^ - F^^A,, (A.ll) 

respectively, with analogous definitions for higer order tensors. The covariant derivative of 
the metric vanishes, 

Dpg^u = 0. (A.12) 

The Lie derivative £a evaluates the change of a tensor field T'^^ along the fiow of vector 
field A''. The coordinate expression of the Lie derivative is 

£aT^. = A^iDpT^",) - TP^iDpA^') + T%{D,Ap). (A.13) 

The Riemann curvature tensor is given by 

i^^pA = 9pr(;, - a,r(;^ + r^^r^, - r^r^,. (a.m) 

The Ricci Tensor is the contraction of the Riemann curvature tensor, 

Rvp — R^upn- (A. 15) 
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Contracting once more yields the scalar curvature, sometimes also referred to as the Ricci 
scalar, 

R = g'^'^R,,. (A.16) 

The fundamental relation between the spacetime metric gi^,^ and the energy momentum 
tensor T^i, is given by Einstein's famous equation, which in standard notation is 

G^. = R^.u - Ig^^uR = SnGT^, + Ag^,. (A.17) 

Whether or not the cosmological term with the cosmological constant A needs to be in- 
cluded is still the subject of debate. 

The energy momentum conservation in general relativity is expressed as the vanishing of 
the divergence of the energy momentum tensor. 



(A.18) 



APPENDIX B 



.THE COORDINATE TRANSFORMATIONS FOR THE 

SPHERE 



Let the coordinates of the sphere be called x = (x, z) and the ones of the stereo- 
graphic projection plane at the south pole X5 = (0,0,— a) be {x^,x^). Let further- 
more the coordinates of the sphere after the rotation of the point Xq = {xo,yo,ZQ) — 
(a sin ^0 cos 00, sin ^0 sin 00, o cos ^0) to the south pole X5 be called x = (x,y,T) and 
the ones from the new stereographic projection plane {x^,x'^). Recall that in this case 
L2 = 4a2 + {x'f + 

The inverse stereographic projection from the coordinates (x^, x"^) to the coordinates of the 
sphere x = {x, y, z) is 

The rotation matrix is given by D — D^DyD'^ and we find 

(cos 00 —sin 00 0\ /cosh^o sinh^oX / cos0o sin0o 0^ 
sin 00 cos 00 1 — sin0o cos0o 

1/ ysinh^o cosh^o/ \ 1, 

(cosh ^0 cos^ 00 + sin^ 0o sin 0o cos 0o (cosh ^0 — 1) sinh 9q cos ipo 
sin 00 cos 00 (cosh ^o — 1) cosh Oq sin^ 0o + cos^ 0o sinh 9q sin 0o 
sinh 6*0 cos 0o sinh 9q sin 0o cosh 9q 



(B.2) 



After the rotation we apply the stereographic projection to the rotated sphere with the 
coordinates x = (x, y, z) and find the transformation to the coordinates (x^, x^) of the new 
projection plane to be 

J. = Jf^l^, i'^J^. (B.3) 

[a — z) [a — z) 
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The explicit transformation of the coordinates {x^,x'^) to the coordinates {x^,x^) and the 
first derivatives of this transformation are 



X — 



~2 
X 



(— cos 9o cos^ 00 + sin^ (j)o)8a^x^ — sin 0o cos 0o(cos 9o + l)8a'^x^ 
+ sin 9o cos 0o4ax^ + sin Oq sin (j)o4:ax'^ + cos ^o(-f'^ — 8a^) 
sin ^0 cos (/)o2a(L^ — 8a^) 
+ sin ^0 cos 0o4ax^ + sin sin 0o4ax^ + cos 9o{L'^ — 8a^) 
— sin 00 cos 0o(cos ^o + 'i-)8a'^x^ + (— cos 9o sin^ 0o + cos^ (l)o)Sa^x'^ 
+ sin ^0 cos 0o4ax^ + sin ^o sin (p^Aax'^ + cos ^o(-^^ — 8a^) 
sin ^0 cos (j)o2a{L'^ — 8a^) 
+ sin 9o cos 0o4aa;^ + sin 9o sin (poAax'^ + cos ^o(-f'^ ~ Sc^^) ' 



(B.4) 



dx^ _ d'P _ _8a2(4a2(cos0o-l)-4aa;2sin(/)osin^o-2((a;i)2-(x2)2)cos2 0ocos2(^o/2)) 
dx"^ (L^ + 4a sin^^o(a^^ cos 00 + sin 0o) + (-^^ — Sa^) cos 6*0)^ 

8a^ (sin^ 0o ( 1 + cos ^0) ( (a:^^ ) ^ - (a^^ ) ^ ) - 4x^ cos 0o {x^ sin 0o (cos ^0 + 1) + « sin 6*0) ) 
(L^ + 4a sin 9q[x^ cos 0o + sin 0o) + {L"^ — Sa^) cos 6*0)^ ' 

dx^ dx'^ 

32a^cos(^o/2)(a;^cos0o — a^^sin^o) 



dx'^ dx^ 

(x-^cos0oCOs(^o/2) + x'^ sin0oCOs(^o/2) + 2asin(^o/2)) 



X 



(L^ + 4a sin ^0(3:^^ cos 0o + x'^ sin 0o) + (L"^ — 8a^) cos ^0)^ 

(B.5) 

The inverse transformation from the coordinates {x^jx"^) to the coordinates {x^jx"^) and 
the first partial derivatives are 



x'^ 



x' 



(— COS 9o cos^ 00 + sin^ 0o)8a^x-^ — sin 0o cos 0o(cos 9o + l)8a^S^ 
— sin ^0 cos 0o4aP — sin ^0 sin 0o4ax^ + cos 9o{L'^ — 8a'^) 



sin ^^0 cos 0o2a(L^ — 8a^)) 
— sin ^^0 cos 0o4ax-'- — sin 9o sin 0o4a5^ + cos 9q{L'^ — 8c 

(B.6) 



— sin 00 cos 0o(cos 9o + l)8a^x-'^ + (— cos 9o sin^ 0o + cos^ 0o)8a^x^ 
— sin 9o cos 0o4aP — sin 9o sin 0o4aP + cos 9o{L'^ — Sa'^) 

sill (9() cos (5()2('/(L^ — 8«^)) 
— sin ^0 cos 0o4aP — sin ^0 sin 0o4aP + cos 9o{L'^ — 8a^) ' 
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B. The Coordinate Transformations for the Sphere 



dx^_dx^_ 8a^(4a^(cos Op - 1) +4a£^ sin 0o sin Oq -2{{x^f- {^f) cos^ 0o cos^(go/2)) 
dx^ 9P (L^ — 4a sin ^o(^''^ cos 0o + sin 0o) + (-^^ — So^) cos ^o)^ 

8a^(sin^ 0o(l+cos 9o){{x'^Y -{x^Y) -Ax^ cos 0o(^ sin 0o(cos 6*0 + 1) - a sin 6*0)) 
(L^ — 4a sin 9o{x^ cos 0o + sin 0o) + (-^^ — 8a^) cos ^0)^ 

5P = = 32a^cos(6'o/2)(x^cos0o - 2^ sin0o) 

(x^ cos0oCos(6'o/2) + 5;^ sin0o cos(^o/2) — 2a sin(^^o/2)) 
(L^ — 4a sin 6*0(2?-^ cos 4>q + 5?^ sin 0o) + (-^'^ — 8a^) cos 6'o)^ 

(B.7) 



APPENDIX C 



.THE COORDINATE TRANSFORMATIONS FOR THE 

HYPERBOLIC PLANE 



Let the coordinates of the hyperbohc plane be called x = {x, y, z) and the ones of the 
stereographic projection plane at the minimum point of the upper sheet = (0, 0, a) 
be Let furthermore the coordinates of the hyperbolic plane after the rotation 

of the point Xq = (xq, yo, -2^0) = (asinh^ocos^o, osinh^osin^o, Qcosh^o) to the minimum 
point xm be called x = (5;, i) and the ones from the new stereographic projection plane 
(x^, x^). Recall that in this case — Aa? — {x^Y — (x^)^. 

The inverse stereographic projection from (x^, x^) to x = (x, 7/, z) is 

4aV 4aV /Sa^ \ 

y^^^, z = a\^-j^-lj. (C.l) 

The rotation matrix of the hyperbolic plane is given hy D = DzDyD~^ and we find 

(cos 00 —sin 00 0\ /cosh^o sinh^oX / cos0o sin^o 0^ 
sin 00 cos 00 1 I I ~ 4>o cos 0o 

1/ \sinh^o cosh^o/ \ 1/ ^^^^ 

(cosh cos^ 00 + sin^ 00 sin 0o cos 0o (cosh ^^o — 1) sinh cos 0o\ 
sin 00 cos 00 (cosh ^^o — 1) cosh Oq sin^ 0o + cos^ 0o sinh 9o sin 0o 
sinh 9o cos 0o sinh Oq sin 0o cosh ^0 



After the rotation we apply the stereographic projection to the rotated hyperbolic plane 
with the coordinates x — {x,y,z) and find the transformation to the coordinates {x^,x^) 

3^ = ^, ^ = (C.3) 
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C. The Coordinate Transformations for the Hyperbolic Plane 



The explicit transformation of the coordinates {x^,x^) to the coordinates {x^,x^) and the 
first derivatives of this transformation are 



(cosh 6q cos^ 00 + sin^ (t)o)Sa'^x^ + sin (pQ cos 0o(cosh 6q — l)8a^a;^ 

— sinh 9q cos (j)Q'iax^ — sinh 6q sin (p^^ax^ + cosh ^0(80^ — L^) 

sinh Oq cos 0o2a(8a^ — L^) 

— sinh 6*0 cos 0o4aa;^ — sinh sin 4>o4ax'^ + cosh 6*0(80^ — L^) 

sin 00 cos 0o(cosh 9o — l)8a'^x^ + (cosh Oq sin^ 0o + cos^ 0o)8a^a;^ 

— sinh ^0 cos (poAax^ — sinh ^0 sin 0o4ax^ + cosh ^o(8a^ — L'^) 

sinh ^0 sin 0o2a(8a^ — L^) 

— sinh 6*0 cos (poiax^ — sinh 6*0 sin (poAax^ + cosh 6*0(80^ — L^) 



dx^ 8a^(4a^(l + cosh ^q) + 2x^x^(2 cosh 60 sin 0o cos 0o — sin(20o))) 
dx^ dx'^ {LP' — 4a sinh 9q{x^ cos 0o + x'^ sin 0o) + (8a^ — LP) cosh ^o)^ 

8a^(-4asinh6lo(x^cos0o + 2^^ sin 0o) + cos(20o)(l - cosh 6lo)((x2)^ - {x'^f)) 
(L^ — 4a sinh ^0(2^^ cos 0o + x"^ sin 0o) + (8a^ — -^^) cosh ^0)^ 

dap ^ ^ -32a^sinh(6'o/2)(x^cos0o -a;^sin0o) 

(2acosh(6'o/2) — sinh(^o/2)(2;^ cos 0o + sin0o)) 



X 



(L^ — 4a sinh 6*0(0;^ cos 0o + x"^ sin 0o) + (8a^ — L^) cosh 6*0)^ 

(C.5) 

The inverse transformation from the coordinates {x^,x^) to the coordinates {x^,x^) and 
the first partial derivatives are 



(cosh ^0 cos^ 00 + sin^ 0o)8a^5:^ + sin 0o cos 0o(cosh Oq — l)8a^5;^ 
+ sinh 9q cos (poiax^ + sinh 9o sin 0o4aP + cosh ^o(8a^ — L^) 

sinh ^0 cos 0o2a(8a^ — L^)) 
+ sinh 9o cos 0o4ax^ + sinh 9o sin (PqAox'^ + cosh ^o(8a2 — L^) 

sin 00 cos 00 (cosh 9o — l)8a^x-'^ + (cosh 9o sin^ 0o + cos^ 0o)8a^x^ 
+ sinh 9o cos 0o4aP + sinh 9o sin (f)oAax^ + cosh ^o(8a2 — L^) 

siuh^?() sin(5o2«(8('/' — L")) 
+ sinh ^0 cos 0o4aP + sinh ^0 sin 0o4aP + cosh ^o(8a2 — L^) 
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dx^ dx^ 8a^(4a^(l + cosh^^o) + 2a;^a;^(2 cosh^o sin0o cos 00 — sin(20o))) 
dx^ ~ ~ (^2 + 4a sinh O^ix^ cos 0o + sin 0o) + (Sa^ - Z2) cosh ^o)^ 

8a^(4asmh^o(5;^ cos0o + x'^ sin^o) + cos(20o)(l — cosh ^^o) ((5^^)^ — {x^Y)) 

I - — - ' — ' 5 

(L^ + 4a sinh ^o(^^ cos 0o + sin 0o) + (8a^ — L^) cosh ^o)^ (^q 
9P = - = 32a^ sinh(6'o/ 2) (x^ cos 0o - sin 0o) 

(2acosh(6'o/2) + sinh(6'o/2)(5'^ cos 0o + x^sin0o)) 
(L^ + 4a sinh 6q{x^ cos 0o + ^ sin 0o) + (Sa^ - L^) cosh 6*0)^ 
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